/ QRZ! Ham Radio 5 / QRZ Ham Radio Callsign Database - Volume 5.iso / files / satelite / orbits / spacetrk.for

< prev

next >

MacOS 8.1

Win98

DOS

browse contents

view JSON data

view as text

+--------+-------------------------+-------------------------+--------+--------+
|00000000| 43 75 74 20 66 69 6c 65 | 20 61 74 20 6c 69 6e 65 |Cut file| at line|
|00000010| 73 20 6f 66 20 74 68 65 | 20 66 6f 72 6d 20 22 2d |s of the| form "-|
|00000020| 2d 20 66 69 6c 65 6e 61 | 6d 65 2e 65 78 74 20 2d |- filena|me.ext -|
|00000030| 2d 22 20 61 6e 64 20 77 | 72 69 74 65 20 74 6f 20 |-" and w|rite to |
|00000040| 66 69 6c 65 0a 77 69 74 | 68 20 74 68 61 74 20 6e |file.wit|h that n|
|00000050| 61 6d 65 2e 0a 0a 2d 2d | 20 41 43 54 41 4e 2e 46 |ame...--| ACTAN.F|
|00000060| 4f 52 20 2d 2d 0a 5c 62 | 65 67 69 6e 7b 76 65 72 |OR --.\b|egin{ver|
|00000070| 62 61 74 69 6d 7d 0a 20 | 20 20 20 20 20 46 55 4e |batim}. |     FUN|
|00000080| 43 54 49 4f 4e 20 20 41 | 43 54 41 4e 28 53 49 4e |CTION  A|CTAN(SIN|
|00000090| 58 2c 43 4f 53 58 29 0a | 20 20 20 20 20 20 43 4f |X,COSX).|      CO|
|000000a0| 4d 4d 4f 4e 2f 43 32 2f | 44 45 32 52 41 2c 50 49 |MMON/C2/|DE2RA,PI|
|000000b0| 2c 50 49 4f 32 2c 54 57 | 4f 50 49 2c 58 33 50 49 |,PIO2,TW|OPI,X3PI|
|000000c0| 4f 32 0a 20 20 20 20 20 | 20 41 43 54 41 4e 3d 30 |O2.     | ACTAN=0|
|000000d0| 2e 0a 20 20 20 20 20 20 | 49 46 20 28 43 4f 53 58 |..      |IF (COSX|
|000000e0| 2e 45 51 2e 30 2e 20 20 | 29 20 47 4f 20 54 4f 20 |.EQ.0.  |) GO TO |
|000000f0| 35 0a 20 20 20 20 20 20 | 49 46 20 28 43 4f 53 58 |5.      |IF (COSX|
|00000100| 2e 47 54 2e 30 2e 20 20 | 29 20 47 4f 20 54 4f 20 |.GT.0.  |) GO TO |
|00000110| 31 0a 20 20 20 20 20 20 | 41 43 54 41 4e 3d 50 49 |1.      |ACTAN=PI|
|00000120| 0a 20 20 20 20 20 20 47 | 4f 20 54 4f 20 37 0a 20 |.      G|O TO 7. |
|00000130| 20 20 20 31 20 49 46 20 | 28 53 49 4e 58 2e 45 51 |   1 IF |(SINX.EQ|
|00000140| 2e 30 2e 20 20 29 20 47 | 4f 20 54 4f 20 38 0a 20 |.0.  ) G|O TO 8. |
|00000150| 20 20 20 20 20 49 46 20 | 28 53 49 4e 58 2e 47 54 |     IF |(SINX.GT|
|00000160| 2e 30 2e 20 20 29 20 47 | 4f 20 54 4f 20 37 0a 20 |.0.  ) G|O TO 7. |
|00000170| 20 20 20 20 20 41 43 54 | 41 4e 3d 54 57 4f 50 49 |     ACT|AN=TWOPI|
|00000180| 0a 20 20 20 20 20 20 47 | 4f 20 54 4f 20 37 0a 20 |.      G|O TO 7. |
|00000190| 20 20 20 35 20 49 46 20 | 28 53 49 4e 58 2e 45 51 |   5 IF |(SINX.EQ|
|000001a0| 2e 30 2e 20 20 29 20 47 | 4f 20 54 4f 20 38 0a 20 |.0.  ) G|O TO 8. |
|000001b0| 20 20 20 20 20 20 49 46 | 20 28 53 49 4e 58 2e 47 |      IF| (SINX.G|
|000001c0| 54 2e 30 2e 20 20 29 20 | 47 4f 20 54 4f 20 36 0a |T.0.  ) |GO TO 6.|
|000001d0| 20 20 20 20 20 20 41 43 | 54 41 4e 3d 58 33 50 49 |      AC|TAN=X3PI|
|000001e0| 4f 32 0a 20 20 20 20 20 | 20 47 4f 20 54 4f 20 38 |O2.     | GO TO 8|
|000001f0| 0a 20 20 20 20 36 20 41 | 43 54 41 4e 3d 50 49 4f |.    6 A|CTAN=PIO|
|00000200| 32 0a 20 20 20 20 20 20 | 47 4f 20 54 4f 20 38 0a |2.      |GO TO 8.|
|00000210| 20 20 20 20 37 20 54 45 | 4d 50 3d 53 49 4e 58 2f |    7 TE|MP=SINX/|
|00000220| 43 4f 53 58 0a 20 20 20 | 20 20 20 41 43 54 41 4e |COSX.   |   ACTAN|
|00000230| 3d 41 43 54 41 4e 2b 41 | 54 41 4e 28 54 45 4d 50 |=ACTAN+A|TAN(TEMP|
|00000240| 29 0a 20 20 20 20 38 20 | 52 45 54 55 52 4e 0a 20 |).    8 |RETURN. |
|00000250| 20 20 20 20 20 45 4e 44 | 0a 5c 65 6e 64 7b 76 65 |     END|.\end{ve|
|00000260| 72 62 61 74 69 6d 7d 0a | 2d 2d 20 44 45 45 50 2e |rbatim}.|-- DEEP.|
|00000270| 46 4f 52 20 2d 2d 0a 5c | 62 65 67 69 6e 7b 76 65 |FOR --.\|begin{ve|
|00000280| 72 62 61 74 69 6d 7d 0a | 2a 20 20 20 20 20 20 44 |rbatim}.|*      D|
|00000290| 45 45 50 20 53 50 41 43 | 45 20 20 20 20 20 20 20 |EEP SPAC|E       |
|000002a0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |        |        |
|000002b0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |        |        |
|000002c0| 20 20 33 31 20 4f 43 54 | 20 38 30 0a 20 20 20 20 |  31 OCT| 80.    |
|000002d0| 20 20 53 55 42 52 4f 55 | 54 49 4e 45 20 44 45 45 |  SUBROU|TINE DEE|
|000002e0| 50 0a 20 20 20 20 20 20 | 43 4f 4d 4d 4f 4e 2f 45 |P.      |COMMON/E|
|000002f0| 31 2f 58 4d 4f 2c 58 4e | 4f 44 45 4f 2c 4f 4d 45 |1/XMO,XN|ODEO,OME|
|00000300| 47 41 4f 2c 45 4f 2c 58 | 49 4e 43 4c 2c 58 4e 4f |GAO,EO,X|INCL,XNO|
|00000310| 2c 58 4e 44 54 32 4f 2c | 0a 20 20 20 20 20 31 20 |,XNDT2O,|.     1 |
|00000320| 20 20 20 20 20 20 20 20 | 20 20 58 4e 44 44 36 4f |        |  XNDD6O|
|00000330| 2c 42 53 54 41 52 2c 58 | 2c 59 2c 5a 2c 58 44 4f |,BSTAR,X|,Y,Z,XDO|
|00000340| 54 2c 59 44 4f 54 2c 5a | 44 4f 54 2c 45 50 4f 43 |T,YDOT,Z|DOT,EPOC|
|00000350| 48 2c 44 53 35 30 0a 20 | 20 20 20 20 20 43 4f 4d |H,DS50. |     COM|
|00000360| 4d 4f 4e 2f 43 31 2f 43 | 4b 32 2c 43 4b 34 2c 45 |MON/C1/C|K2,CK4,E|
|00000370| 36 41 2c 51 4f 4d 53 32 | 54 2c 53 2c 54 4f 54 48 |6A,QOMS2|T,S,TOTH|
|00000380| 52 44 2c 0a 20 20 20 20 | 20 31 20 20 20 20 20 20 |RD,.    | 1      |
|00000390| 20 20 20 20 20 58 4a 33 | 2c 58 4b 45 2c 58 4b 4d |     XJ3|,XKE,XKM|
|000003a0| 50 45 52 2c 58 4d 4e 50 | 44 41 2c 41 45 0a 20 20 |PER,XMNP|DA,AE.  |
|000003b0| 20 20 20 20 43 4f 4d 4d | 4f 4e 2f 43 32 2f 44 45 |    COMM|ON/C2/DE|
|000003c0| 32 52 41 2c 50 49 2c 50 | 49 4f 32 2c 54 57 4f 50 |2RA,PI,P|IO2,TWOP|
|000003d0| 49 2c 58 33 50 49 4f 32 | 0a 20 20 20 20 20 20 44 |I,X3PIO2|.      D|
|000003e0| 4f 55 42 4c 45 20 50 52 | 45 43 49 53 49 4f 4e 20 |OUBLE PR|ECISION |
|000003f0| 45 50 4f 43 48 2c 20 44 | 53 35 30 0a 20 20 20 20 |EPOCH, D|S50.    |
|00000400| 20 20 44 4f 55 42 4c 45 | 20 50 52 45 43 49 53 49 |  DOUBLE| PRECISI|
|00000410| 4f 4e 0a 20 20 20 20 20 | 2a 20 20 20 20 20 44 41 |ON.     |*     DA|
|00000420| 59 2c 50 52 45 45 50 2c | 58 4e 4f 44 43 45 2c 41 |Y,PREEP,|XNODCE,A|
|00000430| 54 49 4d 45 2c 44 45 4c | 54 2c 53 41 56 54 53 4e |TIME,DEL|T,SAVTSN|
|00000440| 2c 53 54 45 50 32 2c 53 | 54 45 50 4e 2c 53 54 45 |,STEP2,S|TEPN,STE|
|00000450| 50 50 0a 20 20 20 20 20 | 20 44 41 54 41 20 20 20 |PP.     | DATA   |
|00000460| 20 20 20 20 20 20 20 20 | 20 20 20 5a 4e 53 2c 20 |        |   ZNS, |
|00000470| 20 20 20 20 20 20 20 20 | 20 20 43 31 53 53 2c 20 |        |  C1SS, |
|00000480| 20 20 20 20 20 20 20 20 | 20 5a 45 53 2f 0a 20 20 |        | ZES/.  |
|00000490| 20 20 20 41 20 20 20 20 | 20 20 20 20 20 20 20 20 |   A    |        |
|000004a0| 20 20 20 20 20 20 31 2e | 31 39 34 35 39 45 2d 35 |      1.|19459E-5|
|000004b0| 2c 20 20 20 20 32 2e 39 | 38 36 34 37 39 37 45 2d |,    2.9|864797E-|
|000004c0| 36 2c 20 2e 30 31 36 37 | 35 2f 0a 20 20 20 20 20 |6, .0167|5/.     |
|000004d0| 20 44 41 54 41 20 20 20 | 20 20 20 20 20 20 20 20 | DATA   |        |
|000004e0| 20 20 20 5a 4e 4c 2c 20 | 20 20 20 20 20 20 20 20 |   ZNL, |        |
|000004f0| 20 20 43 31 4c 2c 20 20 | 20 20 20 20 20 20 20 20 |  C1L,  |        |
|00000500| 20 5a 45 4c 2f 0a 20 20 | 20 20 20 41 20 20 20 20 | ZEL/.  |   A    |
|00000510| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 31 2e |        |      1.|
|00000520| 35 38 33 35 32 31 38 45 | 2d 34 2c 20 20 34 2e 37 |5835218E|-4,  4.7|
|00000530| 39 36 38 30 36 35 45 2d | 37 2c 20 20 2e 30 35 34 |968065E-|7,  .054|
|00000540| 39 30 2f 0a 20 20 20 20 | 20 20 44 41 54 41 20 20 |90/.    |  DATA  |
|00000550| 20 20 20 20 20 20 20 20 | 20 20 20 20 5a 43 4f 53 |        |    ZCOS|
|00000560| 49 53 2c 20 20 20 20 20 | 20 20 20 5a 53 49 4e 49 |IS,     |   ZSINI|
|00000570| 53 2c 20 20 20 20 20 20 | 20 20 5a 53 49 4e 47 53 |S,      |  ZSINGS|
|00000580| 2f 0a 20 20 20 20 20 41 | 20 20 20 20 20 20 20 20 |/.     A|        |
|00000590| 20 20 20 20 20 20 20 20 | 20 20 2e 39 31 37 34 34 |        |  .91744|
|000005a0| 38 36 37 2c 20 20 20 20 | 20 2e 33 39 37 38 35 34 |867,    | .397854|
|000005b0| 31 36 2c 20 20 20 20 20 | 2d 2e 39 38 30 38 38 34 |16,     |-.980884|
|000005c0| 35 38 2f 0a 20 20 20 20 | 20 20 44 41 54 41 20 20 |58/.    |  DATA  |
|000005d0| 20 20 20 20 20 20 20 20 | 20 20 20 20 5a 43 4f 53 |        |    ZCOS|
|000005e0| 47 53 2c 20 20 20 20 20 | 20 20 20 5a 43 4f 53 48 |GS,     |   ZCOSH|
|000005f0| 53 2c 20 20 20 20 20 20 | 20 20 5a 53 49 4e 48 53 |S,      |  ZSINHS|
|00000600| 2f 0a 20 20 20 20 20 41 | 20 20 20 20 20 20 20 20 |/.     A|        |
|00000610| 20 20 20 20 20 20 20 20 | 20 20 2e 31 39 34 35 39 |        |  .19459|
|00000620| 30 35 2c 20 20 20 20 20 | 20 31 2e 30 2c 20 20 20 |05,     | 1.0,   |
|00000630| 20 20 20 20 20 20 20 20 | 30 2e 30 2f 0a 20 20 20 |        |0.0/.   |
|00000640| 20 20 20 44 41 54 41 20 | 51 32 32 2c 51 33 31 2c |   DATA |Q22,Q31,|
|00000650| 51 33 33 2f 31 2e 37 38 | 39 31 36 37 39 45 2d 36 |Q33/1.78|91679E-6|
|00000660| 2c 32 2e 31 34 36 30 37 | 34 38 45 2d 36 2c 32 2e |,2.14607|48E-6,2.|
|00000670| 32 31 32 33 30 31 35 45 | 2d 37 2f 0a 20 20 20 20 |2123015E|-7/.    |
|00000680| 20 20 44 41 54 41 20 47 | 32 32 2c 47 33 32 2f 35 |  DATA G|22,G32/5|
|00000690| 2e 37 36 38 36 33 39 36 | 2c 30 2e 39 35 32 34 30 |.7686396|,0.95240|
|000006a0| 38 39 38 2f 0a 20 20 20 | 20 20 20 44 41 54 41 20 |898/.   |   DATA |
|000006b0| 47 34 34 2c 47 35 32 2f | 31 2e 38 30 31 34 39 39 |G44,G52/|1.801499|
|000006c0| 38 2c 31 2e 30 35 30 38 | 33 33 30 2f 0a 20 20 20 |8,1.0508|330/.   |
|000006d0| 20 20 20 44 41 54 41 20 | 47 35 34 2f 34 2e 34 31 |   DATA |G54/4.41|
|000006e0| 30 38 38 39 38 2f 0a 20 | 20 20 20 20 20 44 41 54 |08898/. |     DAT|
|000006f0| 41 20 52 4f 4f 54 32 32 | 2c 52 4f 4f 54 33 32 2f |A ROOT22|,ROOT32/|
|00000700| 31 2e 37 38 39 31 36 37 | 39 45 2d 36 2c 33 2e 37 |1.789167|9E-6,3.7|
|00000710| 33 39 33 37 39 32 45 2d | 37 2f 0a 20 20 20 20 20 |393792E-|7/.     |
|00000720| 20 44 41 54 41 20 52 4f | 4f 54 34 34 2c 52 4f 4f | DATA RO|OT44,ROO|
|00000730| 54 35 32 2f 37 2e 33 36 | 33 36 39 35 33 45 2d 39 |T52/7.36|36953E-9|
|00000740| 2c 31 2e 31 34 32 38 36 | 33 39 45 2d 37 2f 0a 20 |,1.14286|39E-7/. |
|00000750| 20 20 20 20 20 44 41 54 | 41 20 52 4f 4f 54 35 34 |     DAT|A ROOT54|
|00000760| 2f 32 2e 31 37 36 35 38 | 30 33 45 2d 39 2f 0a 20 |/2.17658|03E-9/. |
|00000770| 20 20 20 20 20 44 41 54 | 41 20 54 48 44 54 2f 34 |     DAT|A THDT/4|
|00000780| 2e 33 37 35 32 36 39 31 | 45 2d 33 2f 0a 0a 2a 20 |.3752691|E-3/..* |
|00000790| 20 20 20 20 45 4e 54 52 | 41 4e 43 45 20 46 4f 52 |    ENTR|ANCE FOR|
|000007a0| 20 44 45 45 50 20 53 50 | 41 43 45 20 49 4e 49 54 | DEEP SP|ACE INIT|
|000007b0| 49 41 4c 49 5a 41 54 49 | 4f 4e 0a 0a 20 20 20 20 |IALIZATI|ON..    |
|000007c0| 20 20 45 4e 54 52 59 20 | 44 50 49 4e 49 54 28 45 |  ENTRY |DPINIT(E|
|000007d0| 51 53 51 2c 53 49 4e 49 | 51 2c 43 4f 53 49 51 2c |QSQ,SINI|Q,COSIQ,|
|000007e0| 52 54 45 51 53 51 2c 41 | 4f 2c 43 4f 53 51 32 2c |RTEQSQ,A|O,COSQ2,|
|000007f0| 53 49 4e 4f 4d 4f 2c 43 | 4f 53 4f 4d 4f 2c 0a 20 |SINOMO,C|OSOMO,. |
|00000800| 20 20 20 20 31 20 20 20 | 20 20 20 20 20 20 42 53 |    1   |      BS|
|00000810| 51 2c 58 4c 4c 44 4f 54 | 2c 4f 4d 47 44 54 2c 58 |Q,XLLDOT|,OMGDT,X|
|00000820| 4e 4f 44 4f 54 2c 58 4e | 4f 44 50 29 0a 20 20 20 |NODOT,XN|ODP).   |
|00000830| 20 20 20 54 48 47 52 3d | 54 48 45 54 41 47 28 45 |   THGR=|THETAG(E|
|00000840| 50 4f 43 48 29 0a 20 20 | 20 20 20 20 45 51 20 3d |POCH).  |    EQ =|
|00000850| 20 45 4f 0a 20 20 20 20 | 20 20 58 4e 51 20 3d 20 | EO.    |  XNQ = |
|00000860| 58 4e 4f 44 50 0a 20 20 | 20 20 20 20 41 51 4e 56 |XNODP.  |    AQNV|
|00000870| 20 3d 20 31 2e 2f 41 4f | 0a 20 20 20 20 20 20 58 | = 1./AO|.      X|
|00000880| 51 4e 43 4c 20 3d 20 58 | 49 4e 43 4c 0a 20 20 20 |QNCL = X|INCL.   |
|00000890| 20 20 20 58 4d 41 4f 3d | 58 4d 4f 0a 20 20 20 20 |   XMAO=|XMO.    |
|000008a0| 20 20 58 50 49 44 4f 54 | 3d 4f 4d 47 44 54 2b 58 |  XPIDOT|=OMGDT+X|
|000008b0| 4e 4f 44 4f 54 0a 20 20 | 20 20 20 20 53 49 4e 51 |NODOT.  |    SINQ|
|000008c0| 20 3d 20 53 49 4e 28 58 | 4e 4f 44 45 4f 29 0a 20 | = SIN(X|NODEO). |
|000008d0| 20 20 20 20 20 43 4f 53 | 51 20 3d 20 43 4f 53 28 |     COS|Q = COS(|
|000008e0| 58 4e 4f 44 45 4f 29 0a | 20 20 20 20 20 20 4f 4d |XNODEO).|      OM|
|000008f0| 45 47 41 51 20 3d 20 4f | 4d 45 47 41 4f 0a 0a 2a |EGAQ = O|MEGAO..*|
|00000900| 20 20 20 20 20 49 4e 49 | 54 49 41 4c 49 5a 45 20 |     INI|TIALIZE |
|00000910| 4c 55 4e 41 52 20 53 4f | 4c 41 52 20 54 45 52 4d |LUNAR SO|LAR TERM|
|00000920| 53 0a 0a 20 20 20 20 35 | 20 44 41 59 3d 44 53 35 |S..    5| DAY=DS5|
|00000930| 30 2b 31 38 32 36 31 2e | 35 44 30 0a 20 20 20 20 |0+18261.|5D0.    |
|00000940| 20 20 49 46 20 28 44 41 | 59 2e 45 51 2e 50 52 45 |  IF (DA|Y.EQ.PRE|
|00000950| 45 50 29 20 20 20 20 47 | 4f 20 54 4f 20 31 30 0a |EP)    G|O TO 10.|
|00000960| 20 20 20 20 20 20 50 52 | 45 45 50 20 3d 20 44 41 |      PR|EEP = DA|
|00000970| 59 0a 20 20 20 20 20 20 | 58 4e 4f 44 43 45 3d 34 |Y.      |XNODCE=4|
|00000980| 2e 35 32 33 36 30 32 30 | 2d 39 2e 32 34 32 32 30 |.5236020|-9.24220|
|00000990| 32 39 45 2d 34 2a 44 41 | 59 0a 20 20 20 20 20 20 |29E-4*DA|Y.      |
|000009a0| 53 54 45 4d 3d 44 53 49 | 4e 20 28 58 4e 4f 44 43 |STEM=DSI|N (XNODC|
|000009b0| 45 29 0a 20 20 20 20 20 | 20 43 54 45 4d 3d 44 43 |E).     | CTEM=DC|
|000009c0| 4f 53 20 28 58 4e 4f 44 | 43 45 29 0a 20 20 20 20 |OS (XNOD|CE).    |
|000009d0| 20 20 5a 43 4f 53 49 4c | 3d 2e 39 31 33 37 35 31 |  ZCOSIL|=.913751|
|000009e0| 36 34 2d 2e 30 33 35 36 | 38 30 39 36 2a 43 54 45 |64-.0356|8096*CTE|
|000009f0| 4d 0a 20 20 20 20 20 20 | 5a 53 49 4e 49 4c 3d 53 |M.      |ZSINIL=S|
|00000a00| 51 52 54 20 28 31 2e 2d | 5a 43 4f 53 49 4c 2a 5a |QRT (1.-|ZCOSIL*Z|
|00000a10| 43 4f 53 49 4c 29 0a 20 | 20 20 20 20 20 5a 53 49 |COSIL). |     ZSI|
|00000a20| 4e 48 4c 3d 20 2e 30 38 | 39 36 38 33 35 31 31 2a |NHL= .08|9683511*|
|00000a30| 53 54 45 4d 2f 5a 53 49 | 4e 49 4c 0a 20 20 20 20 |STEM/ZSI|NIL.    |
|00000a40| 20 20 5a 43 4f 53 48 4c | 3d 53 51 52 54 20 28 31 |  ZCOSHL|=SQRT (1|
|00000a50| 2e 2d 5a 53 49 4e 48 4c | 2a 5a 53 49 4e 48 4c 29 |.-ZSINHL|*ZSINHL)|
|00000a60| 0a 20 20 20 20 20 20 43 | 3d 34 2e 37 31 39 39 36 |.      C|=4.71996|
|00000a70| 37 32 2b 2e 32 32 39 39 | 37 31 35 30 2a 44 41 59 |72+.2299|7150*DAY|
|00000a80| 0a 20 20 20 20 20 20 47 | 41 4d 3d 35 2e 38 33 35 |.      G|AM=5.835|
|00000a90| 31 35 31 34 2b 2e 30 30 | 31 39 34 34 33 36 38 30 |1514+.00|19443680|
|00000aa0| 2a 44 41 59 0a 20 20 20 | 20 20 20 5a 4d 4f 4c 20 |*DAY.   |   ZMOL |
|00000ab0| 3d 20 46 4d 4f 44 32 50 | 28 43 2d 47 41 4d 29 0a |= FMOD2P|(C-GAM).|
|00000ac0| 20 20 20 20 20 20 5a 58 | 3d 20 2e 33 39 37 38 35 |      ZX|= .39785|
|00000ad0| 34 31 36 2a 53 54 45 4d | 2f 5a 53 49 4e 49 4c 0a |416*STEM|/ZSINIL.|
|00000ae0| 20 20 20 20 20 20 5a 59 | 3d 20 5a 43 4f 53 48 4c |      ZY|= ZCOSHL|
|00000af0| 2a 43 54 45 4d 2b 30 2e | 39 31 37 34 34 38 36 37 |*CTEM+0.|91744867|
|00000b00| 2a 5a 53 49 4e 48 4c 2a | 53 54 45 4d 0a 20 20 20 |*ZSINHL*|STEM.   |
|00000b10| 20 20 20 5a 58 3d 41 43 | 54 41 4e 28 5a 58 2c 5a |   ZX=AC|TAN(ZX,Z|
|00000b20| 59 29 0a 20 20 20 20 20 | 20 5a 58 3d 47 41 4d 2b |Y).     | ZX=GAM+|
|00000b30| 5a 58 2d 58 4e 4f 44 43 | 45 0a 20 20 20 20 20 20 |ZX-XNODC|E.      |
|00000b40| 5a 43 4f 53 47 4c 3d 43 | 4f 53 20 28 5a 58 29 0a |ZCOSGL=C|OS (ZX).|
|00000b50| 20 20 20 20 20 20 5a 53 | 49 4e 47 4c 3d 53 49 4e |      ZS|INGL=SIN|
|00000b60| 20 28 5a 58 29 0a 20 20 | 20 20 20 20 5a 4d 4f 53 | (ZX).  |    ZMOS|
|00000b70| 3d 36 2e 32 35 36 35 38 | 33 37 44 30 2b 2e 30 31 |=6.25658|37D0+.01|
|00000b80| 37 32 30 31 39 37 37 44 | 30 2a 44 41 59 0a 20 20 |7201977D|0*DAY.  |
|00000b90| 20 20 20 20 5a 4d 4f 53 | 3d 46 4d 4f 44 32 50 28 |    ZMOS|=FMOD2P(|
|00000ba0| 5a 4d 4f 53 29 0a 0a 2a | 20 20 20 20 20 44 4f 20 |ZMOS)..*|     DO |
|00000bb0| 53 4f 4c 41 52 20 54 45 | 52 4d 53 0a 0a 20 20 20 |SOLAR TE|RMS..   |
|00000bc0| 31 30 20 4c 53 20 3d 20 | 30 0a 20 20 20 20 20 20 |10 LS = |0.      |
|00000bd0| 53 41 56 54 53 4e 3d 31 | 2e 44 32 30 0a 20 20 20 |SAVTSN=1|.D20.   |
|00000be0| 20 20 20 5a 43 4f 53 47 | 3d 5a 43 4f 53 47 53 0a |   ZCOSG|=ZCOSGS.|
|00000bf0| 20 20 20 20 20 20 5a 53 | 49 4e 47 3d 5a 53 49 4e |      ZS|ING=ZSIN|
|00000c00| 47 53 0a 20 20 20 20 20 | 20 5a 43 4f 53 49 3d 5a |GS.     | ZCOSI=Z|
|00000c10| 43 4f 53 49 53 0a 20 20 | 20 20 20 20 5a 53 49 4e |COSIS.  |    ZSIN|
|00000c20| 49 3d 5a 53 49 4e 49 53 | 0a 20 20 20 20 20 20 5a |I=ZSINIS|.      Z|
|00000c30| 43 4f 53 48 3d 43 4f 53 | 51 0a 20 20 20 20 20 20 |COSH=COS|Q.      |
|00000c40| 5a 53 49 4e 48 3d 53 49 | 4e 51 0a 20 20 20 20 20 |ZSINH=SI|NQ.     |
|00000c50| 20 43 43 3d 43 31 53 53 | 0a 20 20 20 20 20 20 5a | CC=C1SS|.      Z|
|00000c60| 4e 3d 5a 4e 53 0a 20 20 | 20 20 20 20 5a 45 3d 5a |N=ZNS.  |    ZE=Z|
|00000c70| 45 53 0a 20 20 20 20 20 | 20 5a 4d 4f 3d 5a 4d 4f |ES.     | ZMO=ZMO|
|00000c80| 53 0a 20 20 20 20 20 20 | 58 4e 4f 49 3d 31 2e 2f |S.      |XNOI=1./|
|00000c90| 58 4e 51 0a 20 20 20 20 | 20 20 41 53 53 49 47 4e |XNQ.    |  ASSIGN|
|00000ca0| 20 33 30 20 54 4f 20 4c | 53 0a 20 20 20 32 30 20 | 30 TO L|S.   20 |
|00000cb0| 41 31 3d 5a 43 4f 53 47 | 2a 5a 43 4f 53 48 2b 5a |A1=ZCOSG|*ZCOSH+Z|
|00000cc0| 53 49 4e 47 2a 5a 43 4f | 53 49 2a 5a 53 49 4e 48 |SING*ZCO|SI*ZSINH|
|00000cd0| 0a 20 20 20 20 20 20 41 | 33 3d 2d 5a 53 49 4e 47 |.      A|3=-ZSING|
|00000ce0| 2a 5a 43 4f 53 48 2b 5a | 43 4f 53 47 2a 5a 43 4f |*ZCOSH+Z|COSG*ZCO|
|00000cf0| 53 49 2a 5a 53 49 4e 48 | 0a 20 20 20 20 20 20 41 |SI*ZSINH|.      A|
|00000d00| 37 3d 2d 5a 43 4f 53 47 | 2a 5a 53 49 4e 48 2b 5a |7=-ZCOSG|*ZSINH+Z|
|00000d10| 53 49 4e 47 2a 5a 43 4f | 53 49 2a 5a 43 4f 53 48 |SING*ZCO|SI*ZCOSH|
|00000d20| 0a 20 20 20 20 20 20 41 | 38 3d 5a 53 49 4e 47 2a |.      A|8=ZSING*|
|00000d30| 5a 53 49 4e 49 0a 20 20 | 20 20 20 20 41 39 3d 5a |ZSINI.  |    A9=Z|
|00000d40| 53 49 4e 47 2a 5a 53 49 | 4e 48 2b 5a 43 4f 53 47 |SING*ZSI|NH+ZCOSG|
|00000d50| 2a 5a 43 4f 53 49 2a 5a | 43 4f 53 48 0a 20 20 20 |*ZCOSI*Z|COSH.   |
|00000d60| 20 20 20 41 31 30 3d 5a | 43 4f 53 47 2a 5a 53 49 |   A10=Z|COSG*ZSI|
|00000d70| 4e 49 0a 20 20 20 20 20 | 20 41 32 3d 20 43 4f 53 |NI.     | A2= COS|
|00000d80| 49 51 2a 41 37 2b 20 53 | 49 4e 49 51 2a 41 38 0a |IQ*A7+ S|INIQ*A8.|
|00000d90| 20 20 20 20 20 20 41 34 | 3d 20 43 4f 53 49 51 2a |      A4|= COSIQ*|
|00000da0| 41 39 2b 20 53 49 4e 49 | 51 2a 41 31 30 0a 20 20 |A9+ SINI|Q*A10.  |
|00000db0| 20 20 20 20 41 35 3d 2d | 20 53 49 4e 49 51 2a 41 |    A5=-| SINIQ*A|
|00000dc0| 37 2b 20 43 4f 53 49 51 | 2a 41 38 0a 20 20 20 20 |7+ COSIQ|*A8.    |
|00000dd0| 20 20 41 36 3d 2d 20 53 | 49 4e 49 51 2a 41 39 2b |  A6=- S|INIQ*A9+|
|00000de0| 20 43 4f 53 49 51 2a 41 | 31 30 0a 43 0a 20 20 20 | COSIQ*A|10.C.   |
|00000df0| 20 20 20 58 31 3d 41 31 | 2a 43 4f 53 4f 4d 4f 2b |   X1=A1|*COSOMO+|
|00000e00| 41 32 2a 53 49 4e 4f 4d | 4f 0a 20 20 20 20 20 20 |A2*SINOM|O.      |
|00000e10| 58 32 3d 41 33 2a 43 4f | 53 4f 4d 4f 2b 41 34 2a |X2=A3*CO|SOMO+A4*|
|00000e20| 53 49 4e 4f 4d 4f 0a 20 | 20 20 20 20 20 58 33 3d |SINOMO. |     X3=|
|00000e30| 2d 41 31 2a 53 49 4e 4f | 4d 4f 2b 41 32 2a 43 4f |-A1*SINO|MO+A2*CO|
|00000e40| 53 4f 4d 4f 0a 20 20 20 | 20 20 20 58 34 3d 2d 41 |SOMO.   |   X4=-A|
|00000e50| 33 2a 53 49 4e 4f 4d 4f | 2b 41 34 2a 43 4f 53 4f |3*SINOMO|+A4*COSO|
|00000e60| 4d 4f 0a 20 20 20 20 20 | 20 58 35 3d 41 35 2a 53 |MO.     | X5=A5*S|
|00000e70| 49 4e 4f 4d 4f 0a 20 20 | 20 20 20 20 58 36 3d 41 |INOMO.  |    X6=A|
|00000e80| 36 2a 53 49 4e 4f 4d 4f | 0a 20 20 20 20 20 20 58 |6*SINOMO|.      X|
|00000e90| 37 3d 41 35 2a 43 4f 53 | 4f 4d 4f 0a 20 20 20 20 |7=A5*COS|OMO.    |
|00000ea0| 20 20 58 38 3d 41 36 2a | 43 4f 53 4f 4d 4f 0a 43 |  X8=A6*|COSOMO.C|
|00000eb0| 0a 20 20 20 20 20 20 5a | 33 31 3d 31 32 2e 2a 58 |.      Z|31=12.*X|
|00000ec0| 31 2a 58 31 2d 33 2e 2a | 58 33 2a 58 33 0a 20 20 |1*X1-3.*|X3*X3.  |
|00000ed0| 20 20 20 20 5a 33 32 3d | 32 34 2e 2a 58 31 2a 58 |    Z32=|24.*X1*X|
|00000ee0| 32 2d 36 2e 2a 58 33 2a | 58 34 0a 20 20 20 20 20 |2-6.*X3*|X4.     |
|00000ef0| 20 5a 33 33 3d 31 32 2e | 2a 58 32 2a 58 32 2d 33 | Z33=12.|*X2*X2-3|
|00000f00| 2e 2a 58 34 2a 58 34 0a | 20 20 20 20 20 20 5a 31 |.*X4*X4.|      Z1|
|00000f10| 3d 33 2e 2a 28 41 31 2a | 41 31 2b 41 32 2a 41 32 |=3.*(A1*|A1+A2*A2|
|00000f20| 29 2b 5a 33 31 2a 45 51 | 53 51 0a 20 20 20 20 20 |)+Z31*EQ|SQ.     |
|00000f30| 20 5a 32 3d 36 2e 2a 28 | 41 31 2a 41 33 2b 41 32 | Z2=6.*(|A1*A3+A2|
|00000f40| 2a 41 34 29 2b 5a 33 32 | 2a 45 51 53 51 0a 20 20 |*A4)+Z32|*EQSQ.  |
|00000f50| 20 20 20 20 5a 33 3d 33 | 2e 2a 28 41 33 2a 41 33 |    Z3=3|.*(A3*A3|
|00000f60| 2b 41 34 2a 41 34 29 2b | 5a 33 33 2a 45 51 53 51 |+A4*A4)+|Z33*EQSQ|
|00000f70| 0a 20 20 20 20 20 20 5a | 31 31 3d 2d 36 2e 2a 41 |.      Z|11=-6.*A|
|00000f80| 31 2a 41 35 2b 45 51 53 | 51 20 2a 28 2d 32 34 2e |1*A5+EQS|Q *(-24.|
|00000f90| 2a 58 31 2a 58 37 2d 36 | 2e 2a 58 33 2a 58 35 29 |*X1*X7-6|.*X3*X5)|
|00000fa0| 0a 20 20 20 20 20 20 5a | 31 32 3d 2d 36 2e 2a 28 |.      Z|12=-6.*(|
|00000fb0| 41 31 2a 41 36 2b 41 33 | 2a 41 35 29 2b 45 51 53 |A1*A6+A3|*A5)+EQS|
|00000fc0| 51 20 2a 28 2d 32 34 2e | 2a 28 58 32 2a 58 37 2b |Q *(-24.|*(X2*X7+|
|00000fd0| 58 31 2a 58 38 29 2d 36 | 2e 2a 28 58 33 2a 58 36 |X1*X8)-6|.*(X3*X6|
|00000fe0| 2b 58 34 2a 58 35 29 29 | 0a 20 20 20 20 20 20 5a |+X4*X5))|.      Z|
|00000ff0| 31 33 3d 2d 36 2e 2a 41 | 33 2a 41 36 2b 45 51 53 |13=-6.*A|3*A6+EQS|
|00001000| 51 20 2a 28 2d 32 34 2e | 2a 58 32 2a 58 38 2d 36 |Q *(-24.|*X2*X8-6|
|00001010| 2e 2a 58 34 2a 58 36 29 | 0a 20 20 20 20 20 20 5a |.*X4*X6)|.      Z|
|00001020| 32 31 3d 36 2e 2a 41 32 | 2a 41 35 2b 45 51 53 51 |21=6.*A2|*A5+EQSQ|
|00001030| 20 2a 28 32 34 2e 2a 58 | 31 2a 58 35 2d 36 2e 2a | *(24.*X|1*X5-6.*|
|00001040| 58 33 2a 58 37 29 0a 20 | 20 20 20 20 20 5a 32 32 |X3*X7). |     Z22|
|00001050| 3d 36 2e 2a 28 41 34 2a | 41 35 2b 41 32 2a 41 36 |=6.*(A4*|A5+A2*A6|
|00001060| 29 2b 45 51 53 51 20 2a | 28 32 34 2e 2a 28 58 32 |)+EQSQ *|(24.*(X2|
|00001070| 2a 58 35 2b 58 31 2a 58 | 36 29 2d 36 2e 2a 28 58 |*X5+X1*X|6)-6.*(X|
|00001080| 34 2a 58 37 2b 58 33 2a | 58 38 29 29 0a 20 20 20 |4*X7+X3*|X8)).   |
|00001090| 20 20 20 5a 32 33 3d 36 | 2e 2a 41 34 2a 41 36 2b |   Z23=6|.*A4*A6+|
|000010a0| 45 51 53 51 20 2a 28 32 | 34 2e 2a 58 32 2a 58 36 |EQSQ *(2|4.*X2*X6|
|000010b0| 2d 36 2e 2a 58 34 2a 58 | 38 29 0a 20 20 20 20 20 |-6.*X4*X|8).     |
|000010c0| 20 5a 31 3d 5a 31 2b 5a | 31 2b 42 53 51 2a 5a 33 | Z1=Z1+Z|1+BSQ*Z3|
|000010d0| 31 0a 20 20 20 20 20 20 | 5a 32 3d 5a 32 2b 5a 32 |1.      |Z2=Z2+Z2|
|000010e0| 2b 42 53 51 2a 5a 33 32 | 0a 20 20 20 20 20 20 5a |+BSQ*Z32|.      Z|
|000010f0| 33 3d 5a 33 2b 5a 33 2b | 42 53 51 2a 5a 33 33 0a |3=Z3+Z3+|BSQ*Z33.|
|00001100| 20 20 20 20 20 20 53 33 | 3d 43 43 2a 58 4e 4f 49 |      S3|=CC*XNOI|
|00001110| 0a 20 20 20 20 20 20 53 | 32 3d 2d 2e 35 2a 53 33 |.      S|2=-.5*S3|
|00001120| 2f 52 54 45 51 53 51 0a | 20 20 20 20 20 20 53 34 |/RTEQSQ.|      S4|
|00001130| 3d 53 33 2a 52 54 45 51 | 53 51 0a 20 20 20 20 20 |=S3*RTEQ|SQ.     |
|00001140| 20 53 31 3d 2d 31 35 2e | 2a 45 51 2a 53 34 0a 20 | S1=-15.|*EQ*S4. |
|00001150| 20 20 20 20 20 53 35 3d | 58 31 2a 58 33 2b 58 32 |     S5=|X1*X3+X2|
|00001160| 2a 58 34 0a 20 20 20 20 | 20 20 53 36 3d 58 32 2a |*X4.    |  S6=X2*|
|00001170| 58 33 2b 58 31 2a 58 34 | 0a 20 20 20 20 20 20 53 |X3+X1*X4|.      S|
|00001180| 37 3d 58 32 2a 58 34 2d | 58 31 2a 58 33 0a 20 20 |7=X2*X4-|X1*X3.  |
|00001190| 20 20 20 20 53 45 3d 53 | 31 2a 5a 4e 2a 53 35 0a |    SE=S|1*ZN*S5.|
|000011a0| 20 20 20 20 20 20 53 49 | 3d 53 32 2a 5a 4e 2a 28 |      SI|=S2*ZN*(|
|000011b0| 5a 31 31 2b 5a 31 33 29 | 0a 20 20 20 20 20 20 53 |Z11+Z13)|.      S|
|000011c0| 4c 3d 2d 5a 4e 2a 53 33 | 2a 28 5a 31 2b 5a 33 2d |L=-ZN*S3|*(Z1+Z3-|
|000011d0| 31 34 2e 2d 36 2e 2a 45 | 51 53 51 29 0a 20 20 20 |14.-6.*E|QSQ).   |
|000011e0| 20 20 20 53 47 48 3d 53 | 34 2a 5a 4e 2a 28 5a 33 |   SGH=S|4*ZN*(Z3|
|000011f0| 31 2b 5a 33 33 2d 36 2e | 29 0a 20 20 20 20 20 20 |1+Z33-6.|).      |
|00001200| 53 48 3d 2d 5a 4e 2a 53 | 32 2a 28 5a 32 31 2b 5a |SH=-ZN*S|2*(Z21+Z|
|00001210| 32 33 29 0a 20 20 20 20 | 20 20 49 46 28 58 51 4e |23).    |  IF(XQN|
|00001220| 43 4c 2e 4c 54 2e 35 2e | 32 33 35 39 38 37 37 45 |CL.LT.5.|2359877E|
|00001230| 2d 32 29 20 53 48 3d 30 | 2e 30 0a 20 20 20 20 20 |-2) SH=0|.0.     |
|00001240| 20 45 45 32 3d 32 2e 2a | 53 31 2a 53 36 0a 20 20 | EE2=2.*|S1*S6.  |
|00001250| 20 20 20 20 45 33 3d 32 | 2e 2a 53 31 2a 53 37 0a |    E3=2|.*S1*S7.|
|00001260| 20 20 20 20 20 20 58 49 | 32 3d 32 2e 2a 53 32 2a |      XI|2=2.*S2*|
|00001270| 5a 31 32 0a 20 20 20 20 | 20 20 58 49 33 3d 32 2e |Z12.    |  XI3=2.|
|00001280| 2a 53 32 2a 28 5a 31 33 | 2d 5a 31 31 29 0a 20 20 |*S2*(Z13|-Z11).  |
|00001290| 20 20 20 20 58 4c 32 3d | 2d 32 2e 2a 53 33 2a 5a |    XL2=|-2.*S3*Z|
|000012a0| 32 0a 20 20 20 20 20 20 | 58 4c 33 3d 2d 32 2e 2a |2.      |XL3=-2.*|
|000012b0| 53 33 2a 28 5a 33 2d 5a | 31 29 0a 20 20 20 20 20 |S3*(Z3-Z|1).     |
|000012c0| 20 58 4c 34 3d 2d 32 2e | 2a 53 33 2a 28 2d 32 31 | XL4=-2.|*S3*(-21|
|000012d0| 2e 2d 39 2e 2a 45 51 53 | 51 29 2a 5a 45 0a 20 20 |.-9.*EQS|Q)*ZE.  |
|000012e0| 20 20 20 20 58 47 48 32 | 3d 32 2e 2a 53 34 2a 5a |    XGH2|=2.*S4*Z|
|000012f0| 33 32 0a 20 20 20 20 20 | 20 58 47 48 33 3d 32 2e |32.     | XGH3=2.|
|00001300| 2a 53 34 2a 28 5a 33 33 | 2d 5a 33 31 29 0a 20 20 |*S4*(Z33|-Z31).  |
|00001310| 20 20 20 20 58 47 48 34 | 3d 2d 31 38 2e 2a 53 34 |    XGH4|=-18.*S4|
|00001320| 2a 5a 45 0a 20 20 20 20 | 20 20 58 48 32 3d 2d 32 |*ZE.    |  XH2=-2|
|00001330| 2e 2a 53 32 2a 5a 32 32 | 0a 20 20 20 20 20 20 58 |.*S2*Z22|.      X|
|00001340| 48 33 3d 2d 32 2e 2a 53 | 32 2a 28 5a 32 33 2d 5a |H3=-2.*S|2*(Z23-Z|
|00001350| 32 31 29 0a 20 20 20 20 | 20 20 47 4f 20 54 4f 20 |21).    |  GO TO |
|00001360| 4c 53 0a 0a 2a 20 20 20 | 20 20 44 4f 20 4c 55 4e |LS..*   |  DO LUN|
|00001370| 41 52 20 54 45 52 4d 53 | 0a 0a 20 20 20 33 30 20 |AR TERMS|..   30 |
|00001380| 53 53 45 20 3d 20 53 45 | 0a 20 20 20 20 20 20 53 |SSE = SE|.      S|
|00001390| 53 49 3d 53 49 0a 20 20 | 20 20 20 20 53 53 4c 3d |SI=SI.  |    SSL=|
|000013a0| 53 4c 0a 20 20 20 20 20 | 20 53 53 48 3d 53 48 2f |SL.     | SSH=SH/|
|000013b0| 53 49 4e 49 51 0a 20 20 | 20 20 20 20 53 53 47 3d |SINIQ.  |    SSG=|
|000013c0| 53 47 48 2d 43 4f 53 49 | 51 2a 53 53 48 0a 20 20 |SGH-COSI|Q*SSH.  |
|000013d0| 20 20 20 20 53 45 32 3d | 45 45 32 0a 20 20 20 20 |    SE2=|EE2.    |
|000013e0| 20 20 53 49 32 3d 58 49 | 32 0a 20 20 20 20 20 20 |  SI2=XI|2.      |
|000013f0| 53 4c 32 3d 58 4c 32 0a | 20 20 20 20 20 20 53 47 |SL2=XL2.|      SG|
|00001400| 48 32 3d 58 47 48 32 0a | 20 20 20 20 20 20 53 48 |H2=XGH2.|      SH|
|00001410| 32 3d 58 48 32 0a 20 20 | 20 20 20 20 53 45 33 3d |2=XH2.  |    SE3=|
|00001420| 45 33 0a 20 20 20 20 20 | 20 53 49 33 3d 58 49 33 |E3.     | SI3=XI3|
|00001430| 0a 20 20 20 20 20 20 53 | 4c 33 3d 58 4c 33 0a 20 |.      S|L3=XL3. |
|00001440| 20 20 20 20 20 53 47 48 | 33 3d 58 47 48 33 0a 20 |     SGH|3=XGH3. |
|00001450| 20 20 20 20 20 53 48 33 | 3d 58 48 33 0a 20 20 20 |     SH3|=XH3.   |
|00001460| 20 20 20 53 4c 34 3d 58 | 4c 34 0a 20 20 20 20 20 |   SL4=X|L4.     |
|00001470| 20 53 47 48 34 3d 58 47 | 48 34 0a 20 20 20 20 20 | SGH4=XG|H4.     |
|00001480| 20 4c 53 3d 31 0a 20 20 | 20 20 20 20 5a 43 4f 53 | LS=1.  |    ZCOS|
|00001490| 47 3d 5a 43 4f 53 47 4c | 0a 20 20 20 20 20 20 5a |G=ZCOSGL|.      Z|
|000014a0| 53 49 4e 47 3d 5a 53 49 | 4e 47 4c 0a 20 20 20 20 |SING=ZSI|NGL.    |
|000014b0| 20 20 5a 43 4f 53 49 3d | 5a 43 4f 53 49 4c 0a 20 |  ZCOSI=|ZCOSIL. |
|000014c0| 20 20 20 20 20 5a 53 49 | 4e 49 3d 5a 53 49 4e 49 |     ZSI|NI=ZSINI|
|000014d0| 4c 0a 20 20 20 20 20 20 | 5a 43 4f 53 48 3d 5a 43 |L.      |ZCOSH=ZC|
|000014e0| 4f 53 48 4c 2a 43 4f 53 | 51 2b 5a 53 49 4e 48 4c |OSHL*COS|Q+ZSINHL|
|000014f0| 2a 53 49 4e 51 0a 20 20 | 20 20 20 20 5a 53 49 4e |*SINQ.  |    ZSIN|
|00001500| 48 3d 53 49 4e 51 2a 5a | 43 4f 53 48 4c 2d 43 4f |H=SINQ*Z|COSHL-CO|
|00001510| 53 51 2a 5a 53 49 4e 48 | 4c 0a 20 20 20 20 20 20 |SQ*ZSINH|L.      |
|00001520| 5a 4e 3d 5a 4e 4c 0a 20 | 20 20 20 20 20 43 43 3d |ZN=ZNL. |     CC=|
|00001530| 43 31 4c 0a 20 20 20 20 | 20 20 5a 45 3d 5a 45 4c |C1L.    |  ZE=ZEL|
|00001540| 0a 20 20 20 20 20 20 5a | 4d 4f 3d 5a 4d 4f 4c 0a |.      Z|MO=ZMOL.|
|00001550| 20 20 20 20 20 20 41 53 | 53 49 47 4e 20 34 30 20 |      AS|SIGN 40 |
|00001560| 54 4f 20 4c 53 0a 20 20 | 20 20 20 20 47 4f 20 54 |TO LS.  |    GO T|
|00001570| 4f 20 32 30 0a 20 20 20 | 34 30 20 53 53 45 20 3d |O 20.   |40 SSE =|
|00001580| 20 53 53 45 2b 53 45 0a | 20 20 20 20 20 20 53 53 | SSE+SE.|      SS|
|00001590| 49 3d 53 53 49 2b 53 49 | 0a 20 20 20 20 20 20 53 |I=SSI+SI|.      S|
|000015a0| 53 4c 3d 53 53 4c 2b 53 | 4c 0a 20 20 20 20 20 20 |SL=SSL+S|L.      |
|000015b0| 53 53 47 3d 53 53 47 2b | 53 47 48 2d 43 4f 53 49 |SSG=SSG+|SGH-COSI|
|000015c0| 51 2f 53 49 4e 49 51 2a | 53 48 0a 20 20 20 20 20 |Q/SINIQ*|SH.     |
|000015d0| 20 53 53 48 3d 53 53 48 | 2b 53 48 2f 53 49 4e 49 | SSH=SSH|+SH/SINI|
|000015e0| 51 0a 0a 2a 20 20 20 20 | 20 47 45 4f 50 4f 54 45 |Q..*    | GEOPOTE|
|000015f0| 4e 54 49 41 4c 20 52 45 | 53 4f 4e 41 4e 43 45 20 |NTIAL RE|SONANCE |
|00001600| 49 4e 49 54 49 41 4c 49 | 5a 41 54 49 4f 4e 20 46 |INITIALI|ZATION F|
|00001610| 4f 52 20 31 32 20 48 4f | 55 52 20 4f 52 42 49 54 |OR 12 HO|UR ORBIT|
|00001620| 53 0a 0a 20 20 20 20 20 | 20 49 52 45 53 46 4c 3d |S..     | IRESFL=|
|00001630| 30 0a 20 20 20 20 20 20 | 49 53 59 4e 46 4c 3d 30 |0.      |ISYNFL=0|
|00001640| 0a 20 20 20 20 20 20 49 | 46 28 58 4e 51 2e 4c 54 |.      I|F(XNQ.LT|
|00001650| 2e 28 2e 30 30 35 32 33 | 35 39 38 37 37 29 2e 41 |.(.00523|59877).A|
|00001660| 4e 44 2e 58 4e 51 2e 47 | 54 2e 28 2e 30 30 33 34 |ND.XNQ.G|T.(.0034|
|00001670| 39 30 36 35 38 35 29 29 | 20 47 4f 20 54 4f 20 37 |906585))| GO TO 7|
|00001680| 30 0a 20 20 20 20 20 20 | 49 46 20 28 58 4e 51 2e |0.      |IF (XNQ.|
|00001690| 4c 54 2e 28 38 2e 32 36 | 45 2d 33 29 20 2e 4f 52 |LT.(8.26|E-3) .OR|
|000016a0| 2e 20 58 4e 51 2e 47 54 | 2e 28 39 2e 32 34 45 2d |. XNQ.GT|.(9.24E-|
|000016b0| 33 29 29 20 20 20 20 52 | 45 54 55 52 4e 0a 20 20 |3))    R|ETURN.  |
|000016c0| 20 20 20 20 49 46 20 28 | 45 51 2e 4c 54 2e 30 2e |    IF (|EQ.LT.0.|
|000016d0| 35 29 20 20 20 20 52 45 | 54 55 52 4e 0a 20 20 20 |5)    RE|TURN.   |
|000016e0| 20 20 20 49 52 45 53 46 | 4c 20 3d 31 0a 20 20 20 |   IRESF|L =1.   |
|000016f0| 20 20 20 45 4f 43 3d 45 | 51 2a 45 51 53 51 0a 20 |   EOC=E|Q*EQSQ. |
|00001700| 20 20 20 20 20 47 32 30 | 31 3d 2d 2e 33 30 36 2d |     G20|1=-.306-|
|00001710| 28 45 51 2d 2e 36 34 29 | 2a 2e 34 34 30 0a 20 20 |(EQ-.64)|*.440.  |
|00001720| 20 20 20 20 49 46 28 45 | 51 2e 47 54 2e 28 2e 36 |    IF(E|Q.GT.(.6|
|00001730| 35 29 29 20 47 4f 20 54 | 4f 20 34 35 0a 20 20 20 |5)) GO T|O 45.   |
|00001740| 20 20 20 47 32 31 31 3d | 33 2e 36 31 36 2d 31 33 |   G211=|3.616-13|
|00001750| 2e 32 34 37 2a 45 51 2b | 31 36 2e 32 39 30 2a 45 |.247*EQ+|16.290*E|
|00001760| 51 53 51 0a 20 20 20 20 | 20 20 47 33 31 30 3d 2d |QSQ.    |  G310=-|
|00001770| 31 39 2e 33 30 32 2b 31 | 31 37 2e 33 39 30 2a 45 |19.302+1|17.390*E|
|00001780| 51 2d 32 32 38 2e 34 31 | 39 2a 45 51 53 51 2b 31 |Q-228.41|9*EQSQ+1|
|00001790| 35 36 2e 35 39 31 2a 45 | 4f 43 0a 20 20 20 20 20 |56.591*E|OC.     |
|000017a0| 20 47 33 32 32 3d 2d 31 | 38 2e 39 30 36 38 2b 31 | G322=-1|8.9068+1|
|000017b0| 30 39 2e 37 39 32 37 2a | 45 51 2d 32 31 34 2e 36 |09.7927*|EQ-214.6|
|000017c0| 33 33 34 2a 45 51 53 51 | 2b 31 34 36 2e 35 38 31 |334*EQSQ|+146.581|
|000017d0| 36 2a 45 4f 43 0a 20 20 | 20 20 20 20 47 34 31 30 |6*EOC.  |    G410|
|000017e0| 3d 2d 34 31 2e 31 32 32 | 2b 32 34 32 2e 36 39 34 |=-41.122|+242.694|
|000017f0| 2a 45 51 2d 34 37 31 2e | 30 39 34 2a 45 51 53 51 |*EQ-471.|094*EQSQ|
|00001800| 2b 33 31 33 2e 39 35 33 | 2a 45 4f 43 0a 20 20 20 |+313.953|*EOC.   |
|00001810| 20 20 20 47 34 32 32 3d | 2d 31 34 36 2e 34 30 37 |   G422=|-146.407|
|00001820| 2b 38 34 31 2e 38 38 30 | 2a 45 51 2d 31 36 32 39 |+841.880|*EQ-1629|
|00001830| 2e 30 31 34 2a 45 51 53 | 51 2b 31 30 38 33 2e 34 |.014*EQS|Q+1083.4|
|00001840| 33 35 2a 45 4f 43 0a 20 | 20 20 20 20 20 47 35 32 |35*EOC. |     G52|
|00001850| 30 3d 2d 35 33 32 2e 31 | 31 34 2b 33 30 31 37 2e |0=-532.1|14+3017.|
|00001860| 39 37 37 2a 45 51 2d 35 | 37 34 30 2a 45 51 53 51 |977*EQ-5|740*EQSQ|
|00001870| 2b 33 37 30 38 2e 32 37 | 36 2a 45 4f 43 0a 20 20 |+3708.27|6*EOC.  |
|00001880| 20 20 20 20 47 4f 20 54 | 4f 20 35 35 0a 20 20 20 |    GO T|O 55.   |
|00001890| 34 35 20 47 32 31 31 3d | 2d 37 32 2e 30 39 39 2b |45 G211=|-72.099+|
|000018a0| 33 33 31 2e 38 31 39 2a | 45 51 2d 35 30 38 2e 37 |331.819*|EQ-508.7|
|000018b0| 33 38 2a 45 51 53 51 2b | 32 36 36 2e 37 32 34 2a |38*EQSQ+|266.724*|
|000018c0| 45 4f 43 0a 20 20 20 20 | 20 20 47 33 31 30 3d 2d |EOC.    |  G310=-|
|000018d0| 33 34 36 2e 38 34 34 2b | 31 35 38 32 2e 38 35 31 |346.844+|1582.851|
|000018e0| 2a 45 51 2d 32 34 31 35 | 2e 39 32 35 2a 45 51 53 |*EQ-2415|.925*EQS|
|000018f0| 51 2b 31 32 34 36 2e 31 | 31 33 2a 45 4f 43 0a 20 |Q+1246.1|13*EOC. |
|00001900| 20 20 20 20 20 47 33 32 | 32 3d 2d 33 34 32 2e 35 |     G32|2=-342.5|
|00001910| 38 35 2b 31 35 35 34 2e | 39 30 38 2a 45 51 2d 32 |85+1554.|908*EQ-2|
|00001920| 33 36 36 2e 38 39 39 2a | 45 51 53 51 2b 31 32 31 |366.899*|EQSQ+121|
|00001930| 35 2e 39 37 32 2a 45 4f | 43 0a 20 20 20 20 20 20 |5.972*EO|C.      |
|00001940| 47 34 31 30 3d 2d 31 30 | 35 32 2e 37 39 37 2b 34 |G410=-10|52.797+4|
|00001950| 37 35 38 2e 36 38 36 2a | 45 51 2d 37 31 39 33 2e |758.686*|EQ-7193.|
|00001960| 39 39 32 2a 45 51 53 51 | 2b 33 36 35 31 2e 39 35 |992*EQSQ|+3651.95|
|00001970| 37 2a 45 4f 43 0a 20 20 | 20 20 20 20 47 34 32 32 |7*EOC.  |    G422|
|00001980| 3d 2d 33 35 38 31 2e 36 | 39 2b 31 36 31 37 38 2e |=-3581.6|9+16178.|
|00001990| 31 31 2a 45 51 2d 32 34 | 34 36 32 2e 37 37 2a 45 |11*EQ-24|462.77*E|
|000019a0| 51 53 51 2b 31 32 34 32 | 32 2e 35 32 2a 45 4f 43 |QSQ+1242|2.52*EOC|
|000019b0| 0a 20 20 20 20 20 20 49 | 46 28 45 51 2e 47 54 2e |.      I|F(EQ.GT.|
|000019c0| 28 2e 37 31 35 29 29 20 | 47 4f 20 54 4f 20 35 30 |(.715)) |GO TO 50|
|000019d0| 0a 20 20 20 20 20 20 47 | 35 32 30 3d 31 34 36 34 |.      G|520=1464|
|000019e0| 2e 37 34 2d 34 36 36 34 | 2e 37 35 2a 45 51 2b 33 |.74-4664|.75*EQ+3|
|000019f0| 37 36 33 2e 36 34 2a 45 | 51 53 51 0a 20 20 20 20 |763.64*E|QSQ.    |
|00001a00| 20 20 47 4f 20 54 4f 20 | 35 35 0a 20 20 20 35 30 |  GO TO |55.   50|
|00001a10| 20 47 35 32 30 3d 2d 35 | 31 34 39 2e 36 36 2b 32 | G520=-5|149.66+2|
|00001a20| 39 39 33 36 2e 39 32 2a | 45 51 2d 35 34 30 38 37 |9936.92*|EQ-54087|
|00001a30| 2e 33 36 2a 45 51 53 51 | 2b 33 31 33 32 34 2e 35 |.36*EQSQ|+31324.5|
|00001a40| 36 2a 45 4f 43 0a 20 20 | 20 35 35 20 49 46 28 45 |6*EOC.  | 55 IF(E|
|00001a50| 51 2e 47 45 2e 28 2e 37 | 29 29 20 47 4f 20 54 4f |Q.GE.(.7|)) GO TO|
|00001a60| 20 36 30 0a 20 20 20 20 | 20 20 47 35 33 33 3d 2d | 60.    |  G533=-|
|00001a70| 39 31 39 2e 32 32 37 37 | 2b 34 39 38 38 2e 36 31 |919.2277|+4988.61|
|00001a80| 2a 45 51 2d 39 30 36 34 | 2e 37 37 2a 45 51 53 51 |*EQ-9064|.77*EQSQ|
|00001a90| 2b 35 35 34 32 2e 32 31 | 2a 45 4f 43 0a 20 20 20 |+5542.21|*EOC.   |
|00001aa0| 20 20 20 47 35 32 31 20 | 3d 20 2d 38 32 32 2e 37 |   G521 |= -822.7|
|00001ab0| 31 30 37 32 2b 34 35 36 | 38 2e 36 31 37 33 2a 45 |1072+456|8.6173*E|
|00001ac0| 51 2d 38 34 39 31 2e 34 | 31 34 36 2a 45 51 53 51 |Q-8491.4|146*EQSQ|
|00001ad0| 2b 35 33 33 37 2e 35 32 | 34 2a 45 4f 43 0a 20 20 |+5337.52|4*EOC.  |
|00001ae0| 20 20 20 20 47 35 33 32 | 20 3d 20 2d 38 35 33 2e |    G532| = -853.|
|00001af0| 36 36 36 2b 34 36 39 30 | 2e 32 35 2a 45 51 2d 38 |666+4690|.25*EQ-8|
|00001b00| 36 32 34 2e 37 37 2a 45 | 51 53 51 2b 35 33 34 31 |624.77*E|QSQ+5341|
|00001b10| 2e 34 2a 45 4f 43 0a 20 | 20 20 20 20 20 47 4f 20 |.4*EOC. |     GO |
|00001b20| 54 4f 20 36 35 0a 20 20 | 20 36 30 20 47 35 33 33 |TO 65.  | 60 G533|
|00001b30| 3d 2d 33 37 39 39 35 2e | 37 38 2b 31 36 31 36 31 |=-37995.|78+16161|
|00001b40| 36 2e 35 32 2a 45 51 2d | 32 32 39 38 33 38 2e 32 |6.52*EQ-|229838.2|
|00001b50| 2a 45 51 53 51 2b 31 30 | 39 33 37 37 2e 39 34 2a |*EQSQ+10|9377.94*|
|00001b60| 45 4f 43 0a 20 20 20 20 | 20 20 47 35 32 31 20 3d |EOC.    |  G521 =|
|00001b70| 20 2d 35 31 37 35 32 2e | 31 30 34 2b 32 31 38 39 | -51752.|104+2189|
|00001b80| 31 33 2e 39 35 2a 45 51 | 2d 33 30 39 34 36 38 2e |13.95*EQ|-309468.|
|00001b90| 31 36 2a 45 51 53 51 2b | 31 34 36 33 34 39 2e 34 |16*EQSQ+|146349.4|
|00001ba0| 32 2a 45 4f 43 0a 20 20 | 20 20 20 20 47 35 33 32 |2*EOC.  |    G532|
|00001bb0| 20 3d 20 2d 34 30 30 32 | 33 2e 38 38 2b 31 37 30 | = -4002|3.88+170|
|00001bc0| 34 37 30 2e 38 39 2a 45 | 51 2d 32 34 32 36 39 39 |470.89*E|Q-242699|
|00001bd0| 2e 34 38 2a 45 51 53 51 | 2b 31 31 35 36 30 35 2e |.48*EQSQ|+115605.|
|00001be0| 38 32 2a 45 4f 43 0a 20 | 20 20 36 35 20 53 49 4e |82*EOC. |  65 SIN|
|00001bf0| 49 32 3d 53 49 4e 49 51 | 2a 53 49 4e 49 51 0a 20 |I2=SINIQ|*SINIQ. |
|00001c00| 20 20 20 20 20 46 32 32 | 30 3d 2e 37 35 2a 28 31 |     F22|0=.75*(1|
|00001c10| 2e 2b 32 2e 2a 43 4f 53 | 49 51 2b 43 4f 53 51 32 |.+2.*COS|IQ+COSQ2|
|00001c20| 29 0a 20 20 20 20 20 20 | 46 32 32 31 3d 31 2e 35 |).      |F221=1.5|
|00001c30| 2a 53 49 4e 49 32 0a 20 | 20 20 20 20 20 46 33 32 |*SINI2. |     F32|
|00001c40| 31 3d 31 2e 38 37 35 2a | 53 49 4e 49 51 2a 28 31 |1=1.875*|SINIQ*(1|
|00001c50| 2e 2d 32 2e 2a 43 4f 53 | 49 51 2d 33 2e 2a 43 4f |.-2.*COS|IQ-3.*CO|
|00001c60| 53 51 32 29 0a 20 20 20 | 20 20 20 46 33 32 32 3d |SQ2).   |   F322=|
|00001c70| 2d 31 2e 38 37 35 2a 53 | 49 4e 49 51 2a 28 31 2e |-1.875*S|INIQ*(1.|
|00001c80| 2b 32 2e 2a 43 4f 53 49 | 51 2d 33 2e 2a 43 4f 53 |+2.*COSI|Q-3.*COS|
|00001c90| 51 32 29 0a 20 20 20 20 | 20 20 46 34 34 31 3d 33 |Q2).    |  F441=3|
|00001ca0| 35 2e 2a 53 49 4e 49 32 | 2a 46 32 32 30 0a 20 20 |5.*SINI2|*F220.  |
|00001cb0| 20 20 20 20 46 34 34 32 | 3d 33 39 2e 33 37 35 30 |    F442|=39.3750|
|00001cc0| 2a 53 49 4e 49 32 2a 53 | 49 4e 49 32 0a 20 20 20 |*SINI2*S|INI2.   |
|00001cd0| 20 20 20 46 35 32 32 3d | 39 2e 38 34 33 37 35 2a |   F522=|9.84375*|
|00001ce0| 53 49 4e 49 51 2a 28 53 | 49 4e 49 32 2a 28 31 2e |SINIQ*(S|INI2*(1.|
|00001cf0| 2d 32 2e 2a 43 4f 53 49 | 51 2d 35 2e 2a 43 4f 53 |-2.*COSI|Q-5.*COS|
|00001d00| 51 32 29 0a 20 20 20 20 | 20 31 20 20 20 20 20 2b |Q2).    | 1     +|
|00001d10| 2e 33 33 33 33 33 33 33 | 33 2a 28 2d 32 2e 2b 34 |.3333333|3*(-2.+4|
|00001d20| 2e 2a 43 4f 53 49 51 2b | 36 2e 2a 43 4f 53 51 32 |.*COSIQ+|6.*COSQ2|
|00001d30| 29 29 0a 20 20 20 20 20 | 20 46 35 32 33 20 3d 20 |)).     | F523 = |
|00001d40| 53 49 4e 49 51 2a 28 34 | 2e 39 32 31 38 37 35 31 |SINIQ*(4|.9218751|
|00001d50| 32 2a 53 49 4e 49 32 2a | 28 2d 32 2e 2d 34 2e 2a |2*SINI2*|(-2.-4.*|
|00001d60| 43 4f 53 49 51 2b 31 30 | 2e 2a 43 4f 53 51 32 29 |COSIQ+10|.*COSQ2)|
|00001d70| 0a 20 20 20 20 20 2a 20 | 20 20 20 20 20 2b 36 2e |.     * |     +6.|
|00001d80| 35 36 32 35 30 30 31 32 | 2a 28 31 2e 2b 32 2e 2a |56250012|*(1.+2.*|
|00001d90| 43 4f 53 49 51 2d 33 2e | 2a 43 4f 53 51 32 29 29 |COSIQ-3.|*COSQ2))|
|00001da0| 0a 20 20 20 20 20 20 46 | 35 34 32 20 3d 20 32 39 |.      F|542 = 29|
|00001db0| 2e 35 33 31 32 35 2a 53 | 49 4e 49 51 2a 28 32 2e |.53125*S|INIQ*(2.|
|00001dc0| 2d 38 2e 2a 43 4f 53 49 | 51 2b 43 4f 53 51 32 2a |-8.*COSI|Q+COSQ2*|
|00001dd0| 28 2d 31 32 2e 2b 38 2e | 2a 43 4f 53 49 51 0a 20 |(-12.+8.|*COSIQ. |
|00001de0| 20 20 20 20 2a 20 20 20 | 20 20 20 2b 31 30 2e 2a |    *   |   +10.*|
|00001df0| 43 4f 53 51 32 29 29 0a | 20 20 20 20 20 20 46 35 |COSQ2)).|      F5|
|00001e00| 34 33 3d 32 39 2e 35 33 | 31 32 35 2a 53 49 4e 49 |43=29.53|125*SINI|
|00001e10| 51 2a 28 2d 32 2e 2d 38 | 2e 2a 43 4f 53 49 51 2b |Q*(-2.-8|.*COSIQ+|
|00001e20| 43 4f 53 51 32 2a 28 31 | 32 2e 2b 38 2e 2a 43 4f |COSQ2*(1|2.+8.*CO|
|00001e30| 53 49 51 2d 31 30 2e 2a | 43 4f 53 51 32 29 29 0a |SIQ-10.*|COSQ2)).|
|00001e40| 20 20 20 20 20 20 58 4e | 4f 32 3d 58 4e 51 2a 58 |      XN|O2=XNQ*X|
|00001e50| 4e 51 0a 20 20 20 20 20 | 20 41 49 4e 56 32 3d 41 |NQ.     | AINV2=A|
|00001e60| 51 4e 56 2a 41 51 4e 56 | 0a 20 20 20 20 20 20 54 |QNV*AQNV|.      T|
|00001e70| 45 4d 50 31 20 3d 20 33 | 2e 2a 58 4e 4f 32 2a 41 |EMP1 = 3|.*XNO2*A|
|00001e80| 49 4e 56 32 0a 20 20 20 | 20 20 20 54 45 4d 50 20 |INV2.   |   TEMP |
|00001e90| 3d 20 54 45 4d 50 31 2a | 52 4f 4f 54 32 32 0a 20 |= TEMP1*|ROOT22. |
|00001ea0| 20 20 20 20 20 44 32 32 | 30 31 20 3d 20 54 45 4d |     D22|01 = TEM|
|00001eb0| 50 2a 46 32 32 30 2a 47 | 32 30 31 0a 20 20 20 20 |P*F220*G|201.    |
|00001ec0| 20 20 44 32 32 31 31 20 | 3d 20 54 45 4d 50 2a 46 |  D2211 |= TEMP*F|
|00001ed0| 32 32 31 2a 47 32 31 31 | 0a 20 20 20 20 20 20 54 |221*G211|.      T|
|00001ee0| 45 4d 50 31 20 3d 20 54 | 45 4d 50 31 2a 41 51 4e |EMP1 = T|EMP1*AQN|
|00001ef0| 56 0a 20 20 20 20 20 20 | 54 45 4d 50 20 3d 20 54 |V.      |TEMP = T|
|00001f00| 45 4d 50 31 2a 52 4f 4f | 54 33 32 0a 20 20 20 20 |EMP1*ROO|T32.    |
|00001f10| 20 20 44 33 32 31 30 20 | 3d 20 54 45 4d 50 2a 46 |  D3210 |= TEMP*F|
|00001f20| 33 32 31 2a 47 33 31 30 | 0a 20 20 20 20 20 20 44 |321*G310|.      D|
|00001f30| 33 32 32 32 20 3d 20 54 | 45 4d 50 2a 46 33 32 32 |3222 = T|EMP*F322|
|00001f40| 2a 47 33 32 32 0a 20 20 | 20 20 20 20 54 45 4d 50 |*G322.  |    TEMP|
|00001f50| 31 20 3d 20 54 45 4d 50 | 31 2a 41 51 4e 56 0a 20 |1 = TEMP|1*AQNV. |
|00001f60| 20 20 20 20 20 54 45 4d | 50 20 3d 20 32 2e 2a 54 |     TEM|P = 2.*T|
|00001f70| 45 4d 50 31 2a 52 4f 4f | 54 34 34 0a 20 20 20 20 |EMP1*ROO|T44.    |
|00001f80| 20 20 44 34 34 31 30 20 | 3d 20 54 45 4d 50 2a 46 |  D4410 |= TEMP*F|
|00001f90| 34 34 31 2a 47 34 31 30 | 0a 20 20 20 20 20 20 44 |441*G410|.      D|
|00001fa0| 34 34 32 32 20 3d 20 54 | 45 4d 50 2a 46 34 34 32 |4422 = T|EMP*F442|
|00001fb0| 2a 47 34 32 32 0a 20 20 | 20 20 20 20 54 45 4d 50 |*G422.  |    TEMP|
|00001fc0| 31 20 3d 20 54 45 4d 50 | 31 2a 41 51 4e 56 0a 20 |1 = TEMP|1*AQNV. |
|00001fd0| 20 20 20 20 20 54 45 4d | 50 20 3d 20 54 45 4d 50 |     TEM|P = TEMP|
|00001fe0| 31 2a 52 4f 4f 54 35 32 | 0a 20 20 20 20 20 20 44 |1*ROOT52|.      D|
|00001ff0| 35 32 32 30 20 3d 20 54 | 45 4d 50 2a 46 35 32 32 |5220 = T|EMP*F522|
|00002000| 2a 47 35 32 30 0a 20 20 | 20 20 20 20 44 35 32 33 |*G520.  |    D523|
|00002010| 32 20 3d 20 54 45 4d 50 | 2a 46 35 32 33 2a 47 35 |2 = TEMP|*F523*G5|
|00002020| 33 32 0a 20 20 20 20 20 | 20 54 45 4d 50 20 3d 20 |32.     | TEMP = |
|00002030| 32 2e 2a 54 45 4d 50 31 | 2a 52 4f 4f 54 35 34 0a |2.*TEMP1|*ROOT54.|
|00002040| 20 20 20 20 20 20 44 35 | 34 32 31 20 3d 20 54 45 |      D5|421 = TE|
|00002050| 4d 50 2a 46 35 34 32 2a | 47 35 32 31 0a 20 20 20 |MP*F542*|G521.   |
|00002060| 20 20 20 44 35 34 33 33 | 20 3d 20 54 45 4d 50 2a |   D5433| = TEMP*|
|00002070| 46 35 34 33 2a 47 35 33 | 33 0a 20 20 20 20 20 20 |F543*G53|3.      |
|00002080| 58 4c 41 4d 4f 20 3d 20 | 58 4d 41 4f 2b 58 4e 4f |XLAMO = |XMAO+XNO|
|00002090| 44 45 4f 2b 58 4e 4f 44 | 45 4f 2d 54 48 47 52 2d |DEO+XNOD|EO-THGR-|
|000020a0| 54 48 47 52 0a 20 20 20 | 20 20 20 42 46 41 43 54 |THGR.   |   BFACT|
|000020b0| 20 3d 20 58 4c 4c 44 4f | 54 2b 58 4e 4f 44 4f 54 | = XLLDO|T+XNODOT|
|000020c0| 2b 58 4e 4f 44 4f 54 2d | 54 48 44 54 2d 54 48 44 |+XNODOT-|THDT-THD|
|000020d0| 54 0a 20 20 20 20 20 20 | 42 46 41 43 54 3d 42 46 |T.      |BFACT=BF|
|000020e0| 41 43 54 2b 53 53 4c 2b | 53 53 48 2b 53 53 48 0a |ACT+SSL+|SSH+SSH.|
|000020f0| 20 20 20 20 20 20 47 4f | 20 54 4f 20 38 30 0a 0a |      GO| TO 80..|
|00002100| 2a 20 20 20 20 20 20 53 | 59 4e 43 48 52 4f 4e 4f |*      S|YNCHRONO|
|00002110| 55 53 20 52 45 53 4f 4e | 41 4e 43 45 20 54 45 52 |US RESON|ANCE TER|
|00002120| 4d 53 20 49 4e 49 54 49 | 41 4c 49 5a 41 54 49 4f |MS INITI|ALIZATIO|
|00002130| 4e 0a 0a 20 20 20 37 30 | 20 49 52 45 53 46 4c 3d |N..   70| IRESFL=|
|00002140| 31 0a 20 20 20 20 20 20 | 49 53 59 4e 46 4c 3d 31 |1.      |ISYNFL=1|
|00002150| 0a 20 20 20 20 20 20 47 | 32 30 30 3d 31 2e 30 2b |.      G|200=1.0+|
|00002160| 45 51 53 51 2a 28 2d 32 | 2e 35 2b 2e 38 31 32 35 |EQSQ*(-2|.5+.8125|
|00002170| 2a 45 51 53 51 29 0a 20 | 20 20 20 20 20 47 33 31 |*EQSQ). |     G31|
|00002180| 30 3d 31 2e 30 2b 32 2e | 30 2a 45 51 53 51 0a 20 |0=1.0+2.|0*EQSQ. |
|00002190| 20 20 20 20 20 47 33 30 | 30 3d 31 2e 30 2b 45 51 |     G30|0=1.0+EQ|
|000021a0| 53 51 2a 28 2d 36 2e 30 | 2b 36 2e 36 30 39 33 37 |SQ*(-6.0|+6.60937|
|000021b0| 2a 45 51 53 51 29 0a 20 | 20 20 20 20 20 46 32 32 |*EQSQ). |     F22|
|000021c0| 30 3d 2e 37 35 2a 28 31 | 2e 2b 43 4f 53 49 51 29 |0=.75*(1|.+COSIQ)|
|000021d0| 2a 28 31 2e 2b 43 4f 53 | 49 51 29 0a 20 20 20 20 |*(1.+COS|IQ).    |
|000021e0| 20 20 46 33 31 31 3d 2e | 39 33 37 35 2a 53 49 4e |  F311=.|9375*SIN|
|000021f0| 49 51 2a 53 49 4e 49 51 | 2a 28 31 2e 2b 33 2e 2a |IQ*SINIQ|*(1.+3.*|
|00002200| 43 4f 53 49 51 29 2d 2e | 37 35 2a 28 31 2e 2b 43 |COSIQ)-.|75*(1.+C|
|00002210| 4f 53 49 51 29 0a 20 20 | 20 20 20 20 46 33 33 30 |OSIQ).  |    F330|
|00002220| 3d 31 2e 2b 43 4f 53 49 | 51 0a 20 20 20 20 20 20 |=1.+COSI|Q.      |
|00002230| 46 33 33 30 3d 31 2e 38 | 37 35 2a 46 33 33 30 2a |F330=1.8|75*F330*|
|00002240| 46 33 33 30 2a 46 33 33 | 30 0a 20 20 20 20 20 20 |F330*F33|0.      |
|00002250| 44 45 4c 31 3d 33 2e 2a | 58 4e 51 2a 58 4e 51 2a |DEL1=3.*|XNQ*XNQ*|
|00002260| 41 51 4e 56 2a 41 51 4e | 56 0a 20 20 20 20 20 20 |AQNV*AQN|V.      |
|00002270| 44 45 4c 32 3d 32 2e 2a | 44 45 4c 31 2a 46 32 32 |DEL2=2.*|DEL1*F22|
|00002280| 30 2a 47 32 30 30 2a 51 | 32 32 0a 20 20 20 20 20 |0*G200*Q|22.     |
|00002290| 20 44 45 4c 33 3d 33 2e | 2a 44 45 4c 31 2a 46 33 | DEL3=3.|*DEL1*F3|
|000022a0| 33 30 2a 47 33 30 30 2a | 51 33 33 2a 41 51 4e 56 |30*G300*|Q33*AQNV|
|000022b0| 0a 20 20 20 20 20 20 44 | 45 4c 31 3d 44 45 4c 31 |.      D|EL1=DEL1|
|000022c0| 2a 46 33 31 31 2a 47 33 | 31 30 2a 51 33 31 2a 41 |*F311*G3|10*Q31*A|
|000022d0| 51 4e 56 0a 20 20 20 20 | 20 20 46 41 53 58 32 3d |QNV.    |  FASX2=|
|000022e0| 2e 31 33 31 33 30 39 30 | 38 0a 20 20 20 20 20 20 |.1313090|8.      |
|000022f0| 46 41 53 58 34 3d 32 2e | 38 38 34 33 31 39 38 0a |FASX4=2.|8843198.|
|00002300| 20 20 20 20 20 20 46 41 | 53 58 36 3d 2e 33 37 34 |      FA|SX6=.374|
|00002310| 34 38 30 38 37 0a 20 20 | 20 20 20 20 58 4c 41 4d |48087.  |    XLAM|
|00002320| 4f 3d 58 4d 41 4f 2b 58 | 4e 4f 44 45 4f 2b 4f 4d |O=XMAO+X|NODEO+OM|
|00002330| 45 47 41 4f 2d 54 48 47 | 52 0a 20 20 20 20 20 20 |EGAO-THG|R.      |
|00002340| 42 46 41 43 54 20 3d 20 | 58 4c 4c 44 4f 54 2b 58 |BFACT = |XLLDOT+X|
|00002350| 50 49 44 4f 54 2d 54 48 | 44 54 0a 20 20 20 20 20 |PIDOT-TH|DT.     |
|00002360| 20 42 46 41 43 54 3d 42 | 46 41 43 54 2b 53 53 4c | BFACT=B|FACT+SSL|
|00002370| 2b 53 53 47 2b 53 53 48 | 0a 20 20 20 38 30 20 58 |+SSG+SSH|.   80 X|
|00002380| 46 41 43 54 3d 42 46 41 | 43 54 2d 58 4e 51 0a 43 |FACT=BFA|CT-XNQ.C|
|00002390| 0a 43 20 20 20 20 20 49 | 4e 49 54 49 41 4c 49 5a |.C     I|NITIALIZ|
|000023a0| 45 20 49 4e 54 45 47 52 | 41 54 4f 52 0a 43 0a 20 |E INTEGR|ATOR.C. |
|000023b0| 20 20 20 20 20 58 4c 49 | 3d 58 4c 41 4d 4f 0a 20 |     XLI|=XLAMO. |
|000023c0| 20 20 20 20 20 58 4e 49 | 3d 58 4e 51 0a 20 20 20 |     XNI|=XNQ.   |
|000023d0| 20 20 20 41 54 49 4d 45 | 3d 30 2e 44 30 0a 20 20 |   ATIME|=0.D0.  |
|000023e0| 20 20 20 20 53 54 45 50 | 50 3d 37 32 30 2e 44 30 |    STEP|P=720.D0|
|000023f0| 0a 20 20 20 20 20 20 53 | 54 45 50 4e 3d 2d 37 32 |.      S|TEPN=-72|
|00002400| 30 2e 44 30 0a 20 20 20 | 20 20 20 53 54 45 50 32 |0.D0.   |   STEP2|
|00002410| 20 3d 20 32 35 39 32 30 | 30 2e 44 30 0a 20 20 20 | = 25920|0.D0.   |
|00002420| 20 20 20 52 45 54 55 52 | 4e 0a 0a 2a 20 20 20 20 |   RETUR|N..*    |
|00002430| 20 45 4e 54 52 41 4e 43 | 45 20 46 4f 52 20 44 45 | ENTRANC|E FOR DE|
|00002440| 45 50 20 53 50 41 43 45 | 20 53 45 43 55 4c 41 52 |EP SPACE| SECULAR|
|00002450| 20 45 46 46 45 43 54 53 | 0a 0a 20 20 20 20 20 20 | EFFECTS|..      |
|00002460| 45 4e 54 52 59 20 44 50 | 53 45 43 28 58 4c 4c 2c |ENTRY DP|SEC(XLL,|
|00002470| 4f 4d 47 41 53 4d 2c 58 | 4e 4f 44 45 53 2c 45 4d |OMGASM,X|NODES,EM|
|00002480| 2c 58 49 4e 43 2c 58 4e | 2c 54 29 0a 20 20 20 20 |,XINC,XN|,T).    |
|00002490| 20 20 58 4c 4c 3d 58 4c | 4c 2b 53 53 4c 2a 54 0a |  XLL=XL|L+SSL*T.|
|000024a0| 20 20 20 20 20 20 4f 4d | 47 41 53 4d 3d 4f 4d 47 |      OM|GASM=OMG|
|000024b0| 41 53 4d 2b 53 53 47 2a | 54 0a 20 20 20 20 20 20 |ASM+SSG*|T.      |
|000024c0| 58 4e 4f 44 45 53 3d 58 | 4e 4f 44 45 53 2b 53 53 |XNODES=X|NODES+SS|
|000024d0| 48 2a 54 0a 20 20 20 20 | 20 20 45 4d 3d 45 4f 2b |H*T.    |  EM=EO+|
|000024e0| 53 53 45 2a 54 0a 20 20 | 20 20 20 20 58 49 4e 43 |SSE*T.  |    XINC|
|000024f0| 3d 58 49 4e 43 4c 2b 53 | 53 49 2a 54 0a 20 20 20 |=XINCL+S|SI*T.   |
|00002500| 20 20 20 49 46 28 58 49 | 4e 43 20 2e 47 45 2e 20 |   IF(XI|NC .GE. |
|00002510| 30 2e 29 20 47 4f 20 54 | 4f 20 39 30 0a 20 20 20 |0.) GO T|O 90.   |
|00002520| 20 20 20 58 49 4e 43 20 | 3d 20 2d 58 49 4e 43 0a |   XINC |= -XINC.|
|00002530| 20 20 20 20 20 20 58 4e | 4f 44 45 53 20 3d 20 58 |      XN|ODES = X|
|00002540| 4e 4f 44 45 53 20 2b 20 | 50 49 0a 20 20 20 20 20 |NODES + |PI.     |
|00002550| 20 4f 4d 47 41 53 4d 20 | 3d 20 4f 4d 47 41 53 4d | OMGASM |= OMGASM|
|00002560| 20 2d 20 50 49 0a 20 20 | 20 39 30 20 49 46 28 49 | - PI.  | 90 IF(I|
|00002570| 52 45 53 46 4c 20 2e 45 | 51 2e 20 30 29 20 52 45 |RESFL .E|Q. 0) RE|
|00002580| 54 55 52 4e 0a 20 20 31 | 30 30 20 49 46 20 28 41 |TURN.  1|00 IF (A|
|00002590| 54 49 4d 45 2e 45 51 2e | 30 2e 44 30 29 20 20 20 |TIME.EQ.|0.D0)   |
|000025a0| 20 47 4f 20 54 4f 20 31 | 37 30 0a 20 20 20 20 20 | GO TO 1|70.     |
|000025b0| 20 49 46 28 54 2e 47 45 | 2e 28 30 2e 44 30 29 2e | IF(T.GE|.(0.D0).|
|000025c0| 41 4e 44 2e 41 54 49 4d | 45 2e 4c 54 2e 28 30 2e |AND.ATIM|E.LT.(0.|
|000025d0| 44 30 29 29 20 47 4f 20 | 54 4f 20 31 37 30 0a 20 |D0)) GO |TO 170. |
|000025e0| 20 20 20 20 20 49 46 28 | 54 2e 4c 54 2e 28 30 2e |     IF(|T.LT.(0.|
|000025f0| 44 30 29 2e 41 4e 44 2e | 41 54 49 4d 45 2e 47 45 |D0).AND.|ATIME.GE|
|00002600| 2e 28 30 2e 44 30 29 29 | 20 47 4f 20 54 4f 20 31 |.(0.D0))| GO TO 1|
|00002610| 37 30 0a 20 20 31 30 35 | 20 49 46 28 44 41 42 53 |70.  105| IF(DABS|
|00002620| 28 54 29 2e 47 45 2e 44 | 41 42 53 28 41 54 49 4d |(T).GE.D|ABS(ATIM|
|00002630| 45 29 29 20 47 4f 20 54 | 4f 20 31 32 30 0a 20 20 |E)) GO T|O 120.  |
|00002640| 20 20 20 20 44 45 4c 54 | 3d 53 54 45 50 50 0a 20 |    DELT|=STEPP. |
|00002650| 20 20 20 20 20 49 46 20 | 28 54 2e 47 45 2e 30 2e |     IF |(T.GE.0.|
|00002660| 44 30 29 20 20 20 20 44 | 45 4c 54 20 3d 20 53 54 |D0)    D|ELT = ST|
|00002670| 45 50 4e 0a 20 20 31 31 | 30 20 41 53 53 49 47 4e |EPN.  11|0 ASSIGN|
|00002680| 20 31 30 30 20 54 4f 20 | 49 52 45 54 0a 20 20 20 | 100 TO |IRET.   |
|00002690| 20 20 20 47 4f 20 54 4f | 20 31 36 30 0a 20 20 31 |   GO TO| 160.  1|
|000026a0| 32 30 20 44 45 4c 54 3d | 53 54 45 50 4e 0a 20 20 |20 DELT=|STEPN.  |
|000026b0| 20 20 20 20 49 46 20 28 | 54 2e 47 54 2e 30 2e 44 |    IF (|T.GT.0.D|
|000026c0| 30 29 20 20 20 20 44 45 | 4c 54 20 3d 20 53 54 45 |0)    DE|LT = STE|
|000026d0| 50 50 0a 20 20 31 32 35 | 20 49 46 20 28 44 41 42 |PP.  125| IF (DAB|
|000026e0| 53 28 54 2d 41 54 49 4d | 45 29 2e 4c 54 2e 53 54 |S(T-ATIM|E).LT.ST|
|000026f0| 45 50 50 29 20 20 20 20 | 47 4f 20 54 4f 20 31 33 |EPP)    |GO TO 13|
|00002700| 30 0a 20 20 20 20 20 20 | 41 53 53 49 47 4e 20 31 |0.      |ASSIGN 1|
|00002710| 32 35 20 54 4f 20 49 52 | 45 54 0a 20 20 20 20 20 |25 TO IR|ET.     |
|00002720| 20 47 4f 20 54 4f 20 31 | 36 30 0a 20 20 31 33 30 | GO TO 1|60.  130|
|00002730| 20 46 54 20 3d 20 54 2d | 41 54 49 4d 45 0a 20 20 | FT = T-|ATIME.  |
|00002740| 20 20 20 20 41 53 53 49 | 47 4e 20 31 34 30 20 54 |    ASSI|GN 140 T|
|00002750| 4f 20 49 52 45 54 4e 0a | 20 20 20 20 20 20 47 4f |O IRETN.|      GO|
|00002760| 20 54 4f 20 31 35 30 0a | 20 20 31 34 30 20 58 4e | TO 150.|  140 XN|
|00002770| 20 3d 20 58 4e 49 2b 58 | 4e 44 4f 54 2a 46 54 2b | = XNI+X|NDOT*FT+|
|00002780| 58 4e 44 44 54 2a 46 54 | 2a 46 54 2a 30 2e 35 0a |XNDDT*FT|*FT*0.5.|
|00002790| 20 20 20 20 20 20 58 4c | 20 3d 20 58 4c 49 2b 58 |      XL| = XLI+X|
|000027a0| 4c 44 4f 54 2a 46 54 2b | 58 4e 44 4f 54 2a 46 54 |LDOT*FT+|XNDOT*FT|
|000027b0| 2a 46 54 2a 30 2e 35 0a | 20 20 20 20 20 20 54 45 |*FT*0.5.|      TE|
|000027c0| 4d 50 20 3d 20 2d 58 4e | 4f 44 45 53 2b 54 48 47 |MP = -XN|ODES+THG|
|000027d0| 52 2b 54 2a 54 48 44 54 | 0a 20 20 20 20 20 20 58 |R+T*THDT|.      X|
|000027e0| 4c 4c 20 3d 20 58 4c 2d | 4f 4d 47 41 53 4d 2b 54 |LL = XL-|OMGASM+T|
|000027f0| 45 4d 50 0a 20 20 20 20 | 20 20 49 46 20 28 49 53 |EMP.    |  IF (IS|
|00002800| 59 4e 46 4c 2e 45 51 2e | 30 29 20 20 20 20 58 4c |YNFL.EQ.|0)    XL|
|00002810| 4c 20 3d 20 58 4c 2b 54 | 45 4d 50 2b 54 45 4d 50 |L = XL+T|EMP+TEMP|
|00002820| 0a 20 20 20 20 20 20 52 | 45 54 55 52 4e 0a 43 0a |.      R|ETURN.C.|
|00002830| 43 20 20 20 20 20 44 4f | 54 20 54 45 52 4d 53 20 |C     DO|T TERMS |
|00002840| 43 41 4c 43 55 4c 41 54 | 45 44 0a 43 0a 20 20 31 |CALCULAT|ED.C.  1|
|00002850| 35 30 20 49 46 20 28 49 | 53 59 4e 46 4c 2e 45 51 |50 IF (I|SYNFL.EQ|
|00002860| 2e 30 29 20 20 20 20 47 | 4f 20 54 4f 20 31 35 32 |.0)    G|O TO 152|
|00002870| 0a 20 20 20 20 20 20 58 | 4e 44 4f 54 3d 44 45 4c |.      X|NDOT=DEL|
|00002880| 31 2a 53 49 4e 20 28 58 | 4c 49 2d 46 41 53 58 32 |1*SIN (X|LI-FASX2|
|00002890| 29 2b 44 45 4c 32 2a 53 | 49 4e 20 28 32 2e 2a 28 |)+DEL2*S|IN (2.*(|
|000028a0| 58 4c 49 2d 46 41 53 58 | 34 29 29 0a 20 20 20 20 |XLI-FASX|4)).    |
|000028b0| 20 31 20 20 20 20 20 2b | 44 45 4c 33 2a 53 49 4e | 1     +|DEL3*SIN|
|000028c0| 20 28 33 2e 2a 28 58 4c | 49 2d 46 41 53 58 36 29 | (3.*(XL|I-FASX6)|
|000028d0| 29 0a 20 20 20 20 20 20 | 58 4e 44 44 54 20 3d 20 |).      |XNDDT = |
|000028e0| 44 45 4c 31 2a 43 4f 53 | 28 58 4c 49 2d 46 41 53 |DEL1*COS|(XLI-FAS|
|000028f0| 58 32 29 0a 20 20 20 20 | 20 2a 20 20 20 20 20 20 |X2).    | *      |
|00002900| 20 2b 32 2e 2a 44 45 4c | 32 2a 43 4f 53 28 32 2e | +2.*DEL|2*COS(2.|
|00002910| 2a 28 58 4c 49 2d 46 41 | 53 58 34 29 29 0a 20 20 |*(XLI-FA|SX4)).  |
|00002920| 20 20 20 2a 20 20 20 20 | 20 20 20 2b 33 2e 2a 44 |   *    |   +3.*D|
|00002930| 45 4c 33 2a 43 4f 53 28 | 33 2e 2a 28 58 4c 49 2d |EL3*COS(|3.*(XLI-|
|00002940| 46 41 53 58 36 29 29 0a | 20 20 20 20 20 20 47 4f |FASX6)).|      GO|
|00002950| 20 54 4f 20 31 35 34 0a | 20 20 31 35 32 20 58 4f | TO 154.|  152 XO|
|00002960| 4d 49 20 3d 20 4f 4d 45 | 47 41 51 2b 4f 4d 47 44 |MI = OME|GAQ+OMGD|
|00002970| 54 2a 41 54 49 4d 45 0a | 20 20 20 20 20 20 58 32 |T*ATIME.|      X2|
|00002980| 4f 4d 49 20 3d 20 58 4f | 4d 49 2b 58 4f 4d 49 0a |OMI = XO|MI+XOMI.|
|00002990| 20 20 20 20 20 20 58 32 | 4c 49 20 3d 20 58 4c 49 |      X2|LI = XLI|
|000029a0| 2b 58 4c 49 0a 20 20 20 | 20 20 20 58 4e 44 4f 54 |+XLI.   |   XNDOT|
|000029b0| 20 3d 20 44 32 32 30 31 | 2a 53 49 4e 28 58 32 4f | = D2201|*SIN(X2O|
|000029c0| 4d 49 2b 58 4c 49 2d 47 | 32 32 29 0a 20 20 20 20 |MI+XLI-G|22).    |
|000029d0| 20 2a 20 20 20 20 20 20 | 20 2b 44 32 32 31 31 2a | *      | +D2211*|
|000029e0| 53 49 4e 28 58 4c 49 2d | 47 32 32 29 0a 20 20 20 |SIN(XLI-|G22).   |
|000029f0| 20 20 2a 20 20 20 20 20 | 20 20 2b 44 33 32 31 30 |  *     |  +D3210|
|00002a00| 2a 53 49 4e 28 58 4f 4d | 49 2b 58 4c 49 2d 47 33 |*SIN(XOM|I+XLI-G3|
|00002a10| 32 29 0a 20 20 20 20 20 | 2a 20 20 20 20 20 20 20 |2).     |*       |
|00002a20| 2b 44 33 32 32 32 2a 53 | 49 4e 28 2d 58 4f 4d 49 |+D3222*S|IN(-XOMI|
|00002a30| 2b 58 4c 49 2d 47 33 32 | 29 0a 20 20 20 20 20 2a |+XLI-G32|).     *|
|00002a40| 20 20 20 20 20 20 20 2b | 44 34 34 31 30 2a 53 49 |       +|D4410*SI|
|00002a50| 4e 28 58 32 4f 4d 49 2b | 58 32 4c 49 2d 47 34 34 |N(X2OMI+|X2LI-G44|
|00002a60| 29 0a 20 20 20 20 20 2a | 20 20 20 20 20 20 20 2b |).     *|       +|
|00002a70| 44 34 34 32 32 2a 53 49 | 4e 28 58 32 4c 49 2d 47 |D4422*SI|N(X2LI-G|
|00002a80| 34 34 29 0a 20 20 20 20 | 20 2a 20 20 20 20 20 20 |44).    | *      |
|00002a90| 20 2b 44 35 32 32 30 2a | 53 49 4e 28 58 4f 4d 49 | +D5220*|SIN(XOMI|
|00002aa0| 2b 58 4c 49 2d 47 35 32 | 29 0a 20 20 20 20 20 2a |+XLI-G52|).     *|
|00002ab0| 20 20 20 20 20 20 20 2b | 44 35 32 33 32 2a 53 49 |       +|D5232*SI|
|00002ac0| 4e 28 2d 58 4f 4d 49 2b | 58 4c 49 2d 47 35 32 29 |N(-XOMI+|XLI-G52)|
|00002ad0| 0a 20 20 20 20 20 2a 20 | 20 20 20 20 20 20 2b 44 |.     * |      +D|
|00002ae0| 35 34 32 31 2a 53 49 4e | 28 58 4f 4d 49 2b 58 32 |5421*SIN|(XOMI+X2|
|00002af0| 4c 49 2d 47 35 34 29 0a | 20 20 20 20 20 2a 20 20 |LI-G54).|     *  |
|00002b00| 20 20 20 20 20 2b 44 35 | 34 33 33 2a 53 49 4e 28 |     +D5|433*SIN(|
|00002b10| 2d 58 4f 4d 49 2b 58 32 | 4c 49 2d 47 35 34 29 0a |-XOMI+X2|LI-G54).|
|00002b20| 20 20 20 20 20 20 58 4e | 44 44 54 20 3d 20 44 32 |      XN|DDT = D2|
|00002b30| 32 30 31 2a 43 4f 53 28 | 58 32 4f 4d 49 2b 58 4c |201*COS(|X2OMI+XL|
|00002b40| 49 2d 47 32 32 29 0a 20 | 20 20 20 20 2a 20 20 20 |I-G22). |    *   |
|00002b50| 20 20 20 20 2b 44 32 32 | 31 31 2a 43 4f 53 28 58 |    +D22|11*COS(X|
|00002b60| 4c 49 2d 47 32 32 29 0a | 20 20 20 20 20 2a 20 20 |LI-G22).|     *  |
|00002b70| 20 20 20 20 20 2b 44 33 | 32 31 30 2a 43 4f 53 28 |     +D3|210*COS(|
|00002b80| 58 4f 4d 49 2b 58 4c 49 | 2d 47 33 32 29 0a 20 20 |XOMI+XLI|-G32).  |
|00002b90| 20 20 20 2a 20 20 20 20 | 20 20 20 2b 44 33 32 32 |   *    |   +D322|
|00002ba0| 32 2a 43 4f 53 28 2d 58 | 4f 4d 49 2b 58 4c 49 2d |2*COS(-X|OMI+XLI-|
|00002bb0| 47 33 32 29 0a 20 20 20 | 20 20 2a 20 20 20 20 20 |G32).   |  *     |
|00002bc0| 20 20 2b 44 35 32 32 30 | 2a 43 4f 53 28 58 4f 4d |  +D5220|*COS(XOM|
|00002bd0| 49 2b 58 4c 49 2d 47 35 | 32 29 0a 20 20 20 20 20 |I+XLI-G5|2).     |
|00002be0| 2a 20 20 20 20 20 20 20 | 2b 44 35 32 33 32 2a 43 |*       |+D5232*C|
|00002bf0| 4f 53 28 2d 58 4f 4d 49 | 2b 58 4c 49 2d 47 35 32 |OS(-XOMI|+XLI-G52|
|00002c00| 29 0a 20 20 20 20 20 2a | 20 20 20 20 20 20 20 2b |).     *|       +|
|00002c10| 32 2e 2a 28 44 34 34 31 | 30 2a 43 4f 53 28 58 32 |2.*(D441|0*COS(X2|
|00002c20| 4f 4d 49 2b 58 32 4c 49 | 2d 47 34 34 29 0a 20 20 |OMI+X2LI|-G44).  |
|00002c30| 20 20 20 2a 20 20 20 20 | 20 20 20 2b 44 34 34 32 |   *    |   +D442|
|00002c40| 32 2a 43 4f 53 28 58 32 | 4c 49 2d 47 34 34 29 0a |2*COS(X2|LI-G44).|
|00002c50| 20 20 20 20 20 2a 20 20 | 20 20 20 20 20 2b 44 35 |     *  |     +D5|
|00002c60| 34 32 31 2a 43 4f 53 28 | 58 4f 4d 49 2b 58 32 4c |421*COS(|XOMI+X2L|
|00002c70| 49 2d 47 35 34 29 0a 20 | 20 20 20 20 2a 20 20 20 |I-G54). |    *   |
|00002c80| 20 20 20 20 2b 44 35 34 | 33 33 2a 43 4f 53 28 2d |    +D54|33*COS(-|
|00002c90| 58 4f 4d 49 2b 58 32 4c | 49 2d 47 35 34 29 29 0a |XOMI+X2L|I-G54)).|
|00002ca0| 20 20 31 35 34 20 58 4c | 44 4f 54 3d 58 4e 49 2b |  154 XL|DOT=XNI+|
|00002cb0| 58 46 41 43 54 0a 20 20 | 20 20 20 20 58 4e 44 44 |XFACT.  |    XNDD|
|00002cc0| 54 20 3d 20 58 4e 44 44 | 54 2a 58 4c 44 4f 54 0a |T = XNDD|T*XLDOT.|
|00002cd0| 20 20 20 20 20 20 47 4f | 20 54 4f 20 49 52 45 54 |      GO| TO IRET|
|00002ce0| 4e 0a 43 0a 43 20 20 20 | 20 20 49 4e 54 45 47 52 |N.C.C   |  INTEGR|
|00002cf0| 41 54 4f 52 0a 43 0a 20 | 20 31 36 30 20 41 53 53 |ATOR.C. | 160 ASS|
|00002d00| 49 47 4e 20 31 36 35 20 | 54 4f 20 49 52 45 54 4e |IGN 165 |TO IRETN|
|00002d10| 0a 20 20 20 20 20 20 47 | 4f 20 54 4f 20 31 35 30 |.      G|O TO 150|
|00002d20| 0a 20 20 31 36 35 20 58 | 4c 49 20 3d 20 58 4c 49 |.  165 X|LI = XLI|
|00002d30| 2b 58 4c 44 4f 54 2a 44 | 45 4c 54 2b 58 4e 44 4f |+XLDOT*D|ELT+XNDO|
|00002d40| 54 2a 53 54 45 50 32 0a | 20 20 20 20 20 20 58 4e |T*STEP2.|      XN|
|00002d50| 49 20 3d 20 58 4e 49 2b | 58 4e 44 4f 54 2a 44 45 |I = XNI+|XNDOT*DE|
|00002d60| 4c 54 2b 58 4e 44 44 54 | 2a 53 54 45 50 32 0a 20 |LT+XNDDT|*STEP2. |
|00002d70| 20 20 20 20 20 41 54 49 | 4d 45 3d 41 54 49 4d 45 |     ATI|ME=ATIME|
|00002d80| 2b 44 45 4c 54 0a 20 20 | 20 20 20 20 47 4f 20 54 |+DELT.  |    GO T|
|00002d90| 4f 20 49 52 45 54 0a 43 | 0a 43 20 20 20 20 20 45 |O IRET.C|.C     E|
|00002da0| 50 4f 43 48 20 52 45 53 | 54 41 52 54 0a 43 0a 20 |POCH RES|TART.C. |
|00002db0| 20 31 37 30 20 49 46 20 | 28 54 2e 47 45 2e 30 2e | 170 IF |(T.GE.0.|
|00002dc0| 44 30 29 20 20 20 20 47 | 4f 20 54 4f 20 31 37 35 |D0)    G|O TO 175|
|00002dd0| 0a 20 20 20 20 20 20 44 | 45 4c 54 3d 53 54 45 50 |.      D|ELT=STEP|
|00002de0| 4e 0a 20 20 20 20 20 20 | 47 4f 20 54 4f 20 31 38 |N.      |GO TO 18|
|00002df0| 30 0a 20 20 31 37 35 20 | 44 45 4c 54 20 3d 20 53 |0.  175 |DELT = S|
|00002e00| 54 45 50 50 0a 20 20 31 | 38 30 20 41 54 49 4d 45 |TEPP.  1|80 ATIME|
|00002e10| 20 3d 20 30 2e 44 30 0a | 20 20 20 20 20 20 58 4e | = 0.D0.|      XN|
|00002e20| 49 3d 58 4e 51 0a 20 20 | 20 20 20 20 58 4c 49 3d |I=XNQ.  |    XLI=|
|00002e30| 58 4c 41 4d 4f 0a 20 20 | 20 20 20 20 47 4f 20 54 |XLAMO.  |    GO T|
|00002e40| 4f 20 31 32 35 0a 43 0a | 43 20 20 20 20 20 45 4e |O 125.C.|C     EN|
|00002e50| 54 52 41 4e 43 45 53 20 | 46 4f 52 20 4c 55 4e 41 |TRANCES |FOR LUNA|
|00002e60| 52 2d 53 4f 4c 41 52 20 | 50 45 52 49 4f 44 49 43 |R-SOLAR |PERIODIC|
|00002e70| 53 0a 43 0a 43 0a 20 20 | 20 20 20 20 45 4e 54 52 |S.C.C.  |    ENTR|
|00002e80| 59 20 44 50 50 45 52 28 | 45 4d 2c 58 49 4e 43 2c |Y DPPER(|EM,XINC,|
|00002e90| 4f 4d 47 41 53 4d 2c 58 | 4e 4f 44 45 53 2c 58 4c |OMGASM,X|NODES,XL|
|00002ea0| 4c 29 0a 20 20 20 20 20 | 20 53 49 4e 49 53 20 3d |L).     | SINIS =|
|00002eb0| 20 53 49 4e 28 58 49 4e | 43 29 0a 20 20 20 20 20 | SIN(XIN|C).     |
|00002ec0| 20 43 4f 53 49 53 20 3d | 20 43 4f 53 28 58 49 4e | COSIS =| COS(XIN|
|00002ed0| 43 29 0a 20 20 20 20 20 | 20 49 46 20 28 44 41 42 |C).     | IF (DAB|
|00002ee0| 53 28 53 41 56 54 53 4e | 2d 54 29 2e 4c 54 2e 28 |S(SAVTSN|-T).LT.(|
|00002ef0| 33 30 2e 44 30 29 29 20 | 20 20 20 47 4f 20 54 4f |30.D0)) |   GO TO|
|00002f00| 20 32 31 30 0a 20 20 20 | 20 20 20 53 41 56 54 53 | 210.   |   SAVTS|
|00002f10| 4e 3d 54 0a 20 20 20 20 | 20 20 5a 4d 3d 5a 4d 4f |N=T.    |  ZM=ZMO|
|00002f20| 53 2b 5a 4e 53 2a 54 0a | 20 20 32 30 35 20 5a 46 |S+ZNS*T.|  205 ZF|
|00002f30| 3d 5a 4d 2b 32 2e 2a 5a | 45 53 2a 53 49 4e 20 28 |=ZM+2.*Z|ES*SIN (|
|00002f40| 5a 4d 29 0a 20 20 20 20 | 20 20 53 49 4e 5a 46 3d |ZM).    |  SINZF=|
|00002f50| 53 49 4e 20 28 5a 46 29 | 0a 20 20 20 20 20 20 46 |SIN (ZF)|.      F|
|00002f60| 32 3d 2e 35 2a 53 49 4e | 5a 46 2a 53 49 4e 5a 46 |2=.5*SIN|ZF*SINZF|
|00002f70| 2d 2e 32 35 0a 20 20 20 | 20 20 20 46 33 3d 2d 2e |-.25.   |   F3=-.|
|00002f80| 35 2a 53 49 4e 5a 46 2a | 43 4f 53 20 28 5a 46 29 |5*SINZF*|COS (ZF)|
|00002f90| 0a 20 20 20 20 20 20 53 | 45 53 3d 53 45 32 2a 46 |.      S|ES=SE2*F|
|00002fa0| 32 2b 53 45 33 2a 46 33 | 0a 20 20 20 20 20 20 53 |2+SE3*F3|.      S|
|00002fb0| 49 53 3d 53 49 32 2a 46 | 32 2b 53 49 33 2a 46 33 |IS=SI2*F|2+SI3*F3|
|00002fc0| 0a 20 20 20 20 20 20 53 | 4c 53 3d 53 4c 32 2a 46 |.      S|LS=SL2*F|
|00002fd0| 32 2b 53 4c 33 2a 46 33 | 2b 53 4c 34 2a 53 49 4e |2+SL3*F3|+SL4*SIN|
|00002fe0| 5a 46 0a 20 20 20 20 20 | 20 53 47 48 53 3d 53 47 |ZF.     | SGHS=SG|
|00002ff0| 48 32 2a 46 32 2b 53 47 | 48 33 2a 46 33 2b 53 47 |H2*F2+SG|H3*F3+SG|
|00003000| 48 34 2a 53 49 4e 5a 46 | 0a 20 20 20 20 20 20 53 |H4*SINZF|.      S|
|00003010| 48 53 3d 53 48 32 2a 46 | 32 2b 53 48 33 2a 46 33 |HS=SH2*F|2+SH3*F3|
|00003020| 0a 20 20 20 20 20 20 5a | 4d 3d 5a 4d 4f 4c 2b 5a |.      Z|M=ZMOL+Z|
|00003030| 4e 4c 2a 54 0a 20 20 20 | 20 20 20 5a 46 3d 5a 4d |NL*T.   |   ZF=ZM|
|00003040| 2b 32 2e 2a 5a 45 4c 2a | 53 49 4e 20 28 5a 4d 29 |+2.*ZEL*|SIN (ZM)|
|00003050| 0a 20 20 20 20 20 20 53 | 49 4e 5a 46 3d 53 49 4e |.      S|INZF=SIN|
|00003060| 20 28 5a 46 29 0a 20 20 | 20 20 20 20 46 32 3d 2e | (ZF).  |    F2=.|
|00003070| 35 2a 53 49 4e 5a 46 2a | 53 49 4e 5a 46 2d 2e 32 |5*SINZF*|SINZF-.2|
|00003080| 35 0a 20 20 20 20 20 20 | 46 33 3d 2d 2e 35 2a 53 |5.      |F3=-.5*S|
|00003090| 49 4e 5a 46 2a 43 4f 53 | 20 28 5a 46 29 0a 20 20 |INZF*COS| (ZF).  |
|000030a0| 20 20 20 20 53 45 4c 3d | 45 45 32 2a 46 32 2b 45 |    SEL=|EE2*F2+E|
|000030b0| 33 2a 46 33 0a 20 20 20 | 20 20 20 53 49 4c 3d 58 |3*F3.   |   SIL=X|
|000030c0| 49 32 2a 46 32 2b 58 49 | 33 2a 46 33 0a 20 20 20 |I2*F2+XI|3*F3.   |
|000030d0| 20 20 20 53 4c 4c 3d 58 | 4c 32 2a 46 32 2b 58 4c |   SLL=X|L2*F2+XL|
|000030e0| 33 2a 46 33 2b 58 4c 34 | 2a 53 49 4e 5a 46 0a 20 |3*F3+XL4|*SINZF. |
|000030f0| 20 20 20 20 20 53 47 48 | 4c 3d 58 47 48 32 2a 46 |     SGH|L=XGH2*F|
|00003100| 32 2b 58 47 48 33 2a 46 | 33 2b 58 47 48 34 2a 53 |2+XGH3*F|3+XGH4*S|
|00003110| 49 4e 5a 46 0a 20 20 20 | 20 20 20 53 48 4c 3d 58 |INZF.   |   SHL=X|
|00003120| 48 32 2a 46 32 2b 58 48 | 33 2a 46 33 0a 20 20 20 |H2*F2+XH|3*F3.   |
|00003130| 20 20 20 50 45 3d 53 45 | 53 2b 53 45 4c 0a 20 20 |   PE=SE|S+SEL.  |
|00003140| 20 20 20 20 50 49 4e 43 | 3d 53 49 53 2b 53 49 4c |    PINC|=SIS+SIL|
|00003150| 0a 20 20 20 20 20 20 50 | 4c 3d 53 4c 53 2b 53 4c |.      P|L=SLS+SL|
|00003160| 4c 0a 20 20 32 31 30 20 | 50 47 48 3d 53 47 48 53 |L.  210 |PGH=SGHS|
|00003170| 2b 53 47 48 4c 0a 20 20 | 20 20 20 20 50 48 3d 53 |+SGHL.  |    PH=S|
|00003180| 48 53 2b 53 48 4c 0a 20 | 20 20 20 20 20 58 49 4e |HS+SHL. |     XIN|
|00003190| 43 20 3d 20 58 49 4e 43 | 2b 50 49 4e 43 0a 20 20 |C = XINC|+PINC.  |
|000031a0| 20 20 20 20 45 4d 20 3d | 20 45 4d 2b 50 45 0a 20 |    EM =| EM+PE. |
|000031b0| 20 20 20 20 20 49 46 28 | 58 51 4e 43 4c 2e 4c 54 |     IF(|XQNCL.LT|
|000031c0| 2e 28 2e 32 29 29 20 47 | 4f 20 54 4f 20 32 32 30 |.(.2)) G|O TO 220|
|000031d0| 0a 20 20 20 20 20 20 47 | 4f 20 54 4f 20 32 31 38 |.      G|O TO 218|
|000031e0| 0a 43 0a 43 20 20 20 20 | 20 41 50 50 4c 59 20 50 |.C.C    | APPLY P|
|000031f0| 45 52 49 4f 44 49 43 53 | 20 44 49 52 45 43 54 4c |ERIODICS| DIRECTL|
|00003200| 59 0a 43 0a 20 20 32 31 | 38 20 50 48 3d 50 48 2f |Y.C.  21|8 PH=PH/|
|00003210| 53 49 4e 49 51 0a 20 20 | 20 20 20 20 50 47 48 3d |SINIQ.  |    PGH=|
|00003220| 50 47 48 2d 43 4f 53 49 | 51 2a 50 48 0a 20 20 20 |PGH-COSI|Q*PH.   |
|00003230| 20 20 20 4f 4d 47 41 53 | 4d 3d 4f 4d 47 41 53 4d |   OMGAS|M=OMGASM|
|00003240| 2b 50 47 48 0a 20 20 20 | 20 20 20 58 4e 4f 44 45 |+PGH.   |   XNODE|
|00003250| 53 3d 58 4e 4f 44 45 53 | 2b 50 48 0a 20 20 20 20 |S=XNODES|+PH.    |
|00003260| 20 20 58 4c 4c 20 3d 20 | 58 4c 4c 2b 50 4c 0a 20 |  XLL = |XLL+PL. |
|00003270| 20 20 20 20 20 47 4f 20 | 54 4f 20 32 33 30 0a 43 |     GO |TO 230.C|
|00003280| 0a 43 20 20 20 20 20 41 | 50 50 4c 59 20 50 45 52 |.C     A|PPLY PER|
|00003290| 49 4f 44 49 43 53 20 57 | 49 54 48 20 4c 59 44 44 |IODICS W|ITH LYDD|
|000032a0| 41 4e 45 20 4d 4f 44 49 | 46 49 43 41 54 49 4f 4e |ANE MODI|FICATION|
|000032b0| 0a 43 0a 20 20 32 32 30 | 20 53 49 4e 4f 4b 3d 53 |.C.  220| SINOK=S|
|000032c0| 49 4e 28 58 4e 4f 44 45 | 53 29 0a 20 20 20 20 20 |IN(XNODE|S).     |
|000032d0| 20 43 4f 53 4f 4b 3d 43 | 4f 53 28 58 4e 4f 44 45 | COSOK=C|OS(XNODE|
|000032e0| 53 29 0a 20 20 20 20 20 | 20 41 4c 46 44 50 3d 53 |S).     | ALFDP=S|
|000032f0| 49 4e 49 53 2a 53 49 4e | 4f 4b 0a 20 20 20 20 20 |INIS*SIN|OK.     |
|00003300| 20 42 45 54 44 50 3d 53 | 49 4e 49 53 2a 43 4f 53 | BETDP=S|INIS*COS|
|00003310| 4f 4b 0a 20 20 20 20 20 | 20 44 41 4c 46 3d 50 48 |OK.     | DALF=PH|
|00003320| 2a 43 4f 53 4f 4b 2b 50 | 49 4e 43 2a 43 4f 53 49 |*COSOK+P|INC*COSI|
|00003330| 53 2a 53 49 4e 4f 4b 0a | 20 20 20 20 20 20 44 42 |S*SINOK.|      DB|
|00003340| 45 54 3d 2d 50 48 2a 53 | 49 4e 4f 4b 2b 50 49 4e |ET=-PH*S|INOK+PIN|
|00003350| 43 2a 43 4f 53 49 53 2a | 43 4f 53 4f 4b 0a 20 20 |C*COSIS*|COSOK.  |
|00003360| 20 20 20 20 41 4c 46 44 | 50 3d 41 4c 46 44 50 2b |    ALFD|P=ALFDP+|
|00003370| 44 41 4c 46 0a 20 20 20 | 20 20 20 42 45 54 44 50 |DALF.   |   BETDP|
|00003380| 3d 42 45 54 44 50 2b 44 | 42 45 54 0a 20 20 20 20 |=BETDP+D|BET.    |
|00003390| 20 20 58 4c 53 20 3d 20 | 58 4c 4c 2b 4f 4d 47 41 |  XLS = |XLL+OMGA|
|000033a0| 53 4d 2b 43 4f 53 49 53 | 2a 58 4e 4f 44 45 53 0a |SM+COSIS|*XNODES.|
|000033b0| 20 20 20 20 20 20 44 4c | 53 3d 50 4c 2b 50 47 48 |      DL|S=PL+PGH|
|000033c0| 2d 50 49 4e 43 2a 58 4e | 4f 44 45 53 2a 53 49 4e |-PINC*XN|ODES*SIN|
|000033d0| 49 53 0a 20 20 20 20 20 | 20 58 4c 53 3d 58 4c 53 |IS.     | XLS=XLS|
|000033e0| 2b 44 4c 53 0a 20 20 20 | 20 20 20 58 4e 4f 44 45 |+DLS.   |   XNODE|
|000033f0| 53 3d 41 43 54 41 4e 28 | 41 4c 46 44 50 2c 42 45 |S=ACTAN(|ALFDP,BE|
|00003400| 54 44 50 29 0a 20 20 20 | 20 20 20 58 4c 4c 20 3d |TDP).   |   XLL =|
|00003410| 20 58 4c 4c 2b 50 4c 0a | 20 20 20 20 20 20 4f 4d | XLL+PL.|      OM|
|00003420| 47 41 53 4d 20 3d 20 58 | 4c 53 2d 58 4c 4c 2d 43 |GASM = X|LS-XLL-C|
|00003430| 4f 53 28 58 49 4e 43 29 | 2a 58 4e 4f 44 45 53 0a |OS(XINC)|*XNODES.|
|00003440| 20 20 32 33 30 20 43 4f | 4e 54 49 4e 55 45 0a 20 |  230 CO|NTINUE. |
|00003450| 20 20 20 20 20 52 45 54 | 55 52 4e 0a 20 20 20 20 |     RET|URN.    |
|00003460| 20 20 45 4e 44 0a 5c 65 | 6e 64 7b 76 65 72 62 61 |  END.\e|nd{verba|
|00003470| 74 69 6d 7d 0a 2d 2d 20 | 44 52 49 56 45 52 2e 46 |tim}.-- |DRIVER.F|
|00003480| 4f 52 20 2d 2d 0a 5c 62 | 65 67 69 6e 7b 76 65 72 |OR --.\b|egin{ver|
|00003490| 62 61 74 69 6d 7d 0a 2a | 20 20 20 20 20 20 44 52 |batim}.*|      DR|
|000034a0| 49 56 45 52 20 20 20 20 | 20 20 20 20 20 20 20 20 |IVER    |        |
|000034b0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |        |        |
|000034c0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |        |        |
|000034d0| 20 20 20 20 20 33 20 4e | 4f 56 20 38 30 0a 0a 2a |     3 N|OV 80..*|
|000034e0| 20 20 20 57 47 53 2d 37 | 32 20 50 48 59 53 49 43 |   WGS-7|2 PHYSIC|
|000034f0| 41 4c 20 41 4e 44 20 47 | 45 4f 50 4f 54 45 4e 54 |AL AND G|EOPOTENT|
|00003500| 49 41 4c 20 43 4f 4e 53 | 54 41 4e 54 53 0a 2a 20 |IAL CONS|TANTS.* |
|00003510| 20 20 20 20 20 20 20 20 | 20 43 4b 32 3d 20 2e 35 |        | CK2= .5|
|00003520| 2a 4a 32 2a 41 45 2a 2a | 32 20 20 20 20 20 43 4b |*J2*AE**|2     CK|
|00003530| 34 3d 2d 2e 33 37 35 2a | 4a 34 2a 41 45 2a 2a 34 |4=-.375*|J4*AE**4|
|00003540| 0a 0a 20 20 20 20 20 20 | 44 4f 55 42 4c 45 20 50 |..      |DOUBLE P|
|00003550| 52 45 43 49 53 49 4f 4e | 20 45 50 4f 43 48 2c 44 |RECISION| EPOCH,D|
|00003560| 53 35 30 0a 20 20 20 20 | 20 20 43 4f 4d 4d 4f 4e |S50.    |  COMMON|
|00003570| 2f 45 31 2f 58 4d 4f 2c | 58 4e 4f 44 45 4f 2c 4f |/E1/XMO,|XNODEO,O|
|00003580| 4d 45 47 41 4f 2c 45 4f | 2c 58 49 4e 43 4c 2c 58 |MEGAO,EO|,XINCL,X|
|00003590| 4e 4f 2c 58 4e 44 54 32 | 4f 2c 58 4e 44 44 36 4f |NO,XNDT2|O,XNDD6O|
|000035a0| 2c 42 53 54 41 52 2c 0a | 20 20 20 20 20 31 20 20 |,BSTAR,.|     1  |
|000035b0| 20 20 20 20 20 20 20 20 | 20 20 58 2c 59 2c 5a 2c |        |  X,Y,Z,|
|000035c0| 58 44 4f 54 2c 59 44 4f | 54 2c 5a 44 4f 54 2c 45 |XDOT,YDO|T,ZDOT,E|
|000035d0| 50 4f 43 48 2c 44 53 35 | 30 0a 20 20 20 20 20 20 |POCH,DS5|0.      |
|000035e0| 43 4f 4d 4d 4f 4e 2f 43 | 31 2f 43 4b 32 2c 43 4b |COMMON/C|1/CK2,CK|
|000035f0| 34 2c 45 36 41 2c 51 4f | 4d 53 32 54 2c 53 2c 54 |4,E6A,QO|MS2T,S,T|
|00003600| 4f 54 48 52 44 2c 0a 20 | 20 20 20 20 31 20 20 20 |OTHRD,. |    1   |
|00003610| 20 20 20 20 20 20 20 20 | 58 4a 33 2c 58 4b 45 2c |        |XJ3,XKE,|
|00003620| 58 4b 4d 50 45 52 2c 58 | 4d 4e 50 44 41 2c 41 45 |XKMPER,X|MNPDA,AE|
|00003630| 0a 20 20 20 20 20 20 43 | 4f 4d 4d 4f 4e 2f 43 32 |.      C|OMMON/C2|
|00003640| 2f 44 45 32 52 41 2c 50 | 49 2c 50 49 4f 32 2c 54 |/DE2RA,P|I,PIO2,T|
|00003650| 57 4f 50 49 2c 58 33 50 | 49 4f 32 0a 20 20 20 20 |WOPI,X3P|IO2.    |
|00003660| 20 20 44 41 54 41 20 49 | 48 47 2f 31 48 47 2f 0a |  DATA I|HG/1HG/.|
|00003670| 20 20 20 20 20 20 44 41 | 54 41 20 44 45 32 52 41 |      DA|TA DE2RA|
|00003680| 2c 45 36 41 2c 50 49 2c | 50 49 4f 32 2c 51 4f 2c |,E6A,PI,|PIO2,QO,|
|00003690| 53 4f 2c 54 4f 54 48 52 | 44 2c 54 57 4f 50 49 2c |SO,TOTHR|D,TWOPI,|
|000036a0| 58 33 50 49 4f 32 2c 58 | 4a 32 2c 58 4a 33 2c 0a |X3PIO2,X|J2,XJ3,.|
|000036b0| 20 20 20 20 20 31 20 20 | 20 20 20 20 20 20 20 20 |     1  |        |
|000036c0| 20 20 58 4a 34 2c 58 4b | 45 2c 58 4b 4d 50 45 52 |  XJ4,XK|E,XKMPER|
|000036d0| 2c 58 4d 4e 50 44 41 2c | 41 45 2f 2e 31 37 34 35 |,XMNPDA,|AE/.1745|
|000036e0| 33 32 39 32 35 45 2d 31 | 2c 31 2e 45 2d 36 2c 0a |32925E-1|,1.E-6,.|
|000036f0| 20 20 20 20 20 32 20 20 | 20 20 20 20 20 20 20 20 |     2  |        |
|00003700| 20 20 33 2e 31 34 31 35 | 39 32 36 35 2c 31 2e 35 |  3.1415|9265,1.5|
|00003710| 37 30 37 39 36 33 33 2c | 31 32 30 2e 30 2c 37 38 |7079633,|120.0,78|
|00003720| 2e 30 2c 2e 36 36 36 36 | 36 36 36 37 2c 0a 20 20 |.0,.6666|6667,.  |
|00003730| 20 20 20 34 20 20 20 20 | 20 20 20 20 20 20 20 20 |   4    |        |
|00003740| 36 2e 32 38 33 31 38 35 | 33 2c 34 2e 37 31 32 33 |6.283185|3,4.7123|
|00003750| 38 38 39 38 2c 31 2e 30 | 38 32 36 31 36 45 2d 33 |8898,1.0|82616E-3|
|00003760| 2c 2d 2e 32 35 33 38 38 | 31 45 2d 35 2c 0a 20 20 |,-.25388|1E-5,.  |
|00003770| 20 20 20 35 20 20 20 20 | 20 20 20 20 20 20 20 20 |   5    |        |
|00003780| 2d 31 2e 36 35 35 39 37 | 45 2d 36 2c 2e 37 34 33 |-1.65597|E-6,.743|
|00003790| 36 36 39 31 36 31 45 2d | 31 2c 36 33 37 38 2e 31 |669161E-|1,6378.1|
|000037a0| 33 35 2c 31 34 34 30 2e | 2c 31 2e 2f 0a 20 20 20 |35,1440.|,1./.   |
|000037b0| 20 20 20 44 49 4d 45 4e | 53 49 4f 4e 20 49 53 45 |   DIMEN|SION ISE|
|000037c0| 54 28 35 29 0a 20 20 20 | 20 20 20 43 48 41 52 41 |T(5).   |   CHARA|
|000037d0| 43 54 45 52 20 41 42 55 | 46 2a 38 30 28 32 29 0a |CTER ABU|F*80(2).|
|000037e0| 20 20 20 20 20 20 44 41 | 54 41 20 28 49 53 45 54 |      DA|TA (ISET|
|000037f0| 28 49 29 2c 49 3d 31 2c | 35 29 2f 33 48 53 47 50 |(I),I=1,|5)/3HSGP|
|00003800| 2c 34 48 53 47 50 34 2c | 34 48 53 44 50 34 2c 34 |,4HSGP4,|4HSDP4,4|
|00003810| 48 53 47 50 38 2c 34 48 | 53 44 50 38 2f 0a 0a 2a |HSGP8,4H|SDP8/..*|
|00003820| 20 20 20 20 20 53 45 4c | 45 43 54 20 45 50 48 45 |     SEL|ECT EPHE|
|00003830| 4d 45 52 49 53 20 54 59 | 50 45 20 41 4e 44 20 4f |MERIS TY|PE AND O|
|00003840| 55 54 50 55 54 20 54 49 | 4d 45 53 0a 0a 20 20 20 |UTPUT TI|MES..   |
|00003850| 20 20 20 43 4b 32 3d 2e | 35 2a 58 4a 32 2a 41 45 |   CK2=.|5*XJ2*AE|
|00003860| 2a 2a 32 0a 20 20 20 20 | 20 20 43 4b 34 3d 2d 2e |**2.    |  CK4=-.|
|00003870| 33 37 35 2a 58 4a 34 2a | 41 45 2a 2a 34 0a 20 20 |375*XJ4*|AE**4.  |
|00003880| 20 20 20 20 51 4f 4d 53 | 32 54 3d 28 28 51 4f 2d |    QOMS|2T=((QO-|
|00003890| 53 4f 29 2a 41 45 2f 58 | 4b 4d 50 45 52 29 2a 2a |SO)*AE/X|KMPER)**|
|000038a0| 34 0a 20 20 20 20 20 20 | 53 3d 41 45 2a 28 31 2e |4.      |S=AE*(1.|
|000038b0| 2b 53 4f 2f 58 4b 4d 50 | 45 52 29 0a 20 20 20 20 |+SO/XKMP|ER).    |
|000038c0| 32 20 52 45 41 44 20 28 | 35 2c 37 30 30 29 20 49 |2 READ (|5,700) I|
|000038d0| 45 50 54 2c 20 54 53 2c | 54 46 2c 44 45 4c 54 0a |EPT, TS,|TF,DELT.|
|000038e0| 20 20 20 20 20 20 49 46 | 28 49 45 50 54 2e 4c 45 |      IF|(IEPT.LE|
|000038f0| 2e 30 29 20 53 54 4f 50 | 0a 20 20 20 20 20 20 49 |.0) STOP|.      I|
|00003900| 44 45 45 50 3d 30 0a 0a | 2a 20 20 20 20 20 52 45 |DEEP=0..|*     RE|
|00003910| 41 44 20 49 4e 20 4d 45 | 41 4e 20 45 4c 45 4d 45 |AD IN ME|AN ELEME|
|00003920| 4e 54 53 20 46 52 4f 4d | 20 32 20 43 41 52 44 20 |NTS FROM| 2 CARD |
|00003930| 54 28 54 52 41 4e 53 29 | 20 4f 52 20 47 28 49 4e |T(TRANS)| OR G(IN|
|00003940| 54 45 52 4e 29 20 46 4f | 52 4d 41 54 0a 0a 20 20 |TERN) FO|RMAT..  |
|00003950| 20 20 20 20 52 45 41 44 | 20 28 35 2c 37 30 36 29 |    READ| (5,706)|
|00003960| 20 41 42 55 46 0a 20 20 | 20 20 20 20 44 45 43 4f | ABUF.  |    DECO|
|00003970| 44 45 28 41 42 55 46 28 | 31 29 2c 37 30 37 29 20 |DE(ABUF(|1),707) |
|00003980| 20 49 54 59 50 45 0a 20 | 20 20 20 20 20 49 46 28 | ITYPE. |     IF(|
|00003990| 49 54 59 50 45 2e 45 51 | 2e 49 48 47 29 20 47 4f |ITYPE.EQ|.IHG) GO|
|000039a0| 20 54 4f 20 35 0a 20 20 | 20 20 20 20 44 45 43 4f | TO 5.  |    DECO|
|000039b0| 44 45 20 28 41 42 55 46 | 2c 37 30 32 29 20 45 50 |DE (ABUF|,702) EP|
|000039c0| 4f 43 48 2c 58 4e 44 54 | 32 4f 2c 58 4e 44 44 36 |OCH,XNDT|2O,XNDD6|
|000039d0| 4f 2c 49 45 58 50 2c 42 | 53 54 41 52 2c 49 42 45 |O,IEXP,B|STAR,IBE|
|000039e0| 58 50 2c 58 49 4e 43 4c | 2c 0a 20 20 20 20 20 31 |XP,XINCL|,.     1|
|000039f0| 20 20 20 20 20 20 20 20 | 58 4e 4f 44 45 4f 2c 45 |        |XNODEO,E|
|00003a00| 4f 2c 4f 4d 45 47 41 4f | 2c 58 4d 4f 2c 58 4e 4f |O,OMEGAO|,XMO,XNO|
|00003a10| 0a 20 20 20 20 20 20 47 | 4f 20 54 4f 20 37 0a 20 |.      G|O TO 7. |
|00003a20| 20 20 20 35 20 44 45 43 | 4f 44 45 28 41 42 55 46 |   5 DEC|ODE(ABUF|
|00003a30| 2c 37 30 31 29 20 45 50 | 4f 43 48 2c 58 4d 4f 2c |,701) EP|OCH,XMO,|
|00003a40| 58 4e 4f 44 45 4f 2c 4f | 4d 45 47 41 4f 2c 45 4f |XNODEO,O|MEGAO,EO|
|00003a50| 2c 58 49 4e 43 4c 2c 58 | 4e 4f 2c 58 4e 44 54 32 |,XINCL,X|NO,XNDT2|
|00003a60| 4f 2c 0a 20 20 20 20 20 | 31 20 20 20 20 20 20 20 |O,.     |1       |
|00003a70| 20 58 4e 44 44 36 4f 2c | 49 45 58 50 2c 42 53 54 | XNDD6O,|IEXP,BST|
|00003a80| 41 52 2c 49 42 45 58 50 | 0a 20 20 20 37 20 20 49 |AR,IBEXP|.   7  I|
|00003a90| 46 28 58 4e 4f 2e 4c 45 | 2e 30 2e 29 20 53 54 4f |F(XNO.LE|.0.) STO|
|00003aa0| 50 0a 20 20 20 20 20 20 | 57 52 49 54 45 28 36 2c |P.      |WRITE(6,|
|00003ab0| 37 30 34 29 20 41 42 55 | 46 2c 49 53 45 54 28 49 |704) ABU|F,ISET(I|
|00003ac0| 45 50 54 29 0a 20 20 20 | 20 20 20 49 46 28 49 45 |EPT).   |   IF(IE|
|00003ad0| 50 54 2e 47 54 2e 35 29 | 20 47 4f 20 54 4f 20 39 |PT.GT.5)| GO TO 9|
|00003ae0| 30 30 0a 20 20 20 20 20 | 20 58 4e 44 44 36 4f 3d |00.     | XNDD6O=|
|00003af0| 58 4e 44 44 36 4f 2a 28 | 31 30 2e 2a 2a 49 45 58 |XNDD6O*(|10.**IEX|
|00003b00| 50 29 0a 20 20 20 20 20 | 20 58 4e 4f 44 45 4f 3d |P).     | XNODEO=|
|00003b10| 58 4e 4f 44 45 4f 2a 44 | 45 32 52 41 0a 20 20 20 |XNODEO*D|E2RA.   |
|00003b20| 20 20 20 4f 4d 45 47 41 | 4f 3d 4f 4d 45 47 41 4f |   OMEGA|O=OMEGAO|
|00003b30| 2a 44 45 32 52 41 0a 20 | 20 20 20 20 20 58 4d 4f |*DE2RA. |     XMO|
|00003b40| 3d 58 4d 4f 2a 44 45 32 | 52 41 0a 20 20 20 20 20 |=XMO*DE2|RA.     |
|00003b50| 20 58 49 4e 43 4c 3d 58 | 49 4e 43 4c 2a 44 45 32 | XINCL=X|INCL*DE2|
|00003b60| 52 41 0a 20 20 20 20 20 | 20 54 45 4d 50 3d 54 57 |RA.     | TEMP=TW|
|00003b70| 4f 50 49 2f 58 4d 4e 50 | 44 41 2f 58 4d 4e 50 44 |OPI/XMNP|DA/XMNPD|
|00003b80| 41 0a 20 20 20 20 20 20 | 58 4e 4f 3d 58 4e 4f 2a |A.      |XNO=XNO*|
|00003b90| 54 45 4d 50 2a 58 4d 4e | 50 44 41 0a 20 20 20 20 |TEMP*XMN|PDA.    |
|00003ba0| 20 20 58 4e 44 54 32 4f | 3d 58 4e 44 54 32 4f 2a |  XNDT2O|=XNDT2O*|
|00003bb0| 54 45 4d 50 0a 20 20 20 | 20 20 20 58 4e 44 44 36 |TEMP.   |   XNDD6|
|00003bc0| 4f 3d 58 4e 44 44 36 4f | 2a 54 45 4d 50 2f 58 4d |O=XNDD6O|*TEMP/XM|
|00003bd0| 4e 50 44 41 0a 0a 2a 20 | 20 20 20 20 49 4e 50 55 |NPDA..* |    INPU|
|00003be0| 54 20 43 48 45 43 4b 20 | 46 4f 52 20 50 45 52 49 |T CHECK |FOR PERI|
|00003bf0| 4f 44 20 56 53 20 45 50 | 48 45 4d 45 52 49 53 20 |OD VS EP|HEMERIS |
|00003c00| 53 45 4c 45 43 54 45 44 | 0a 2a 20 20 20 20 20 50 |SELECTED|.*     P|
|00003c10| 45 52 49 4f 44 20 47 45 | 20 32 32 35 20 4d 49 4e |ERIOD GE| 225 MIN|
|00003c20| 55 54 45 53 20 20 49 53 | 20 44 45 45 50 20 53 50 |UTES  IS| DEEP SP|
|00003c30| 41 43 45 0a 0a 20 20 20 | 20 20 20 41 31 3d 28 58 |ACE..   |   A1=(X|
|00003c40| 4b 45 2f 58 4e 4f 29 2a | 2a 54 4f 54 48 52 44 0a |KE/XNO)*|*TOTHRD.|
|00003c50| 20 20 20 20 20 20 54 45 | 4d 50 3d 31 2e 35 2a 43 |      TE|MP=1.5*C|
|00003c60| 4b 32 2a 28 33 2e 2a 43 | 4f 53 28 58 49 4e 43 4c |K2*(3.*C|OS(XINCL|
|00003c70| 29 2a 2a 32 2d 31 2e 29 | 2f 28 31 2e 2d 45 4f 2a |)**2-1.)|/(1.-EO*|
|00003c80| 45 4f 29 2a 2a 31 2e 35 | 0a 20 20 20 20 20 20 44 |EO)**1.5|.      D|
|00003c90| 45 4c 31 3d 54 45 4d 50 | 2f 28 41 31 2a 41 31 29 |EL1=TEMP|/(A1*A1)|
|00003ca0| 0a 20 20 20 20 20 20 41 | 4f 3d 41 31 2a 28 31 2e |.      A|O=A1*(1.|
|00003cb0| 2d 44 45 4c 31 2a 28 2e | 35 2a 54 4f 54 48 52 44 |-DEL1*(.|5*TOTHRD|
|00003cc0| 2b 44 45 4c 31 2a 28 31 | 2e 2b 31 33 34 2e 2f 38 |+DEL1*(1|.+134./8|
|00003cd0| 31 2e 2a 44 45 4c 31 29 | 29 29 0a 20 20 20 20 20 |1.*DEL1)|)).     |
|00003ce0| 20 44 45 4c 4f 3d 54 45 | 4d 50 2f 28 41 4f 2a 41 | DELO=TE|MP/(AO*A|
|00003cf0| 4f 29 0a 20 20 20 20 20 | 20 58 4e 4f 44 50 3d 58 |O).     | XNODP=X|
|00003d00| 4e 4f 2f 28 31 2e 2b 44 | 45 4c 4f 29 0a 20 20 20 |NO/(1.+D|ELO).   |
|00003d10| 20 20 20 49 46 28 28 54 | 57 4f 50 49 2f 58 4e 4f |   IF((T|WOPI/XNO|
|00003d20| 44 50 2f 58 4d 4e 50 44 | 41 29 20 2e 47 45 2e 20 |DP/XMNPD|A) .GE. |
|00003d30| 2e 31 35 36 32 35 29 20 | 49 44 45 45 50 3d 31 0a |.15625) |IDEEP=1.|
|00003d40| 0a 20 20 20 20 20 20 42 | 53 54 41 52 3d 42 53 54 |.      B|STAR=BST|
|00003d50| 41 52 2a 28 31 30 2e 2a | 2a 49 42 45 58 50 29 2f |AR*(10.*|*IBEXP)/|
|00003d60| 41 45 0a 20 20 20 20 20 | 20 54 53 49 4e 43 45 3d |AE.     | TSINCE=|
|00003d70| 54 53 0a 20 20 20 20 20 | 20 49 46 4c 41 47 3d 31 |TS.     | IFLAG=1|
|00003d80| 0a 20 20 20 20 20 20 49 | 46 28 49 44 45 45 50 20 |.      I|F(IDEEP |
|00003d90| 2e 45 51 2e 20 31 20 2e | 41 4e 44 2e 20 28 49 45 |.EQ. 1 .|AND. (IE|
|00003da0| 50 54 20 2e 45 51 2e 20 | 31 20 2e 4f 52 2e 20 49 |PT .EQ. |1 .OR. I|
|00003db0| 45 50 54 20 2e 45 51 2e | 20 32 0a 20 20 20 20 20 |EPT .EQ.| 2.     |
|00003dc0| 31 20 20 20 20 20 20 20 | 20 20 20 20 2e 4f 52 2e |1       |    .OR.|
|00003dd0| 20 49 45 50 54 20 2e 45 | 51 2e 20 34 29 29 20 47 | IEPT .E|Q. 4)) G|
|00003de0| 4f 20 54 4f 20 38 30 30 | 0a 20 20 20 20 39 20 49 |O TO 800|.    9 I|
|00003df0| 46 28 49 44 45 45 50 20 | 2e 45 51 2e 20 30 20 2e |F(IDEEP |.EQ. 0 .|
|00003e00| 41 4e 44 2e 20 28 49 45 | 50 54 20 2e 45 51 2e 20 |AND. (IE|PT .EQ. |
|00003e10| 33 20 2e 4f 52 2e 20 49 | 45 50 54 20 2e 45 51 2e |3 .OR. I|EPT .EQ.|
|00003e20| 20 35 29 29 0a 20 20 20 | 20 20 31 20 20 20 20 20 | 5)).   |  1     |
|00003e30| 20 20 20 20 20 20 47 4f | 20 54 4f 20 38 35 30 0a |      GO| TO 850.|
|00003e40| 20 20 20 31 30 20 47 4f | 20 54 4f 20 28 32 31 2c |   10 GO| TO (21,|
|00003e50| 32 32 2c 32 33 2c 32 34 | 2c 32 35 29 2c 20 49 45 |22,23,24|,25), IE|
|00003e60| 50 54 0a 20 20 20 32 31 | 20 43 41 4c 4c 20 53 47 |PT.   21| CALL SG|
|00003e70| 50 28 49 46 4c 41 47 2c | 54 53 49 4e 43 45 29 0a |P(IFLAG,|TSINCE).|
|00003e80| 20 20 20 20 20 20 47 4f | 20 54 4f 20 36 30 0a 20 |      GO| TO 60. |
|00003e90| 20 20 32 32 20 43 41 4c | 4c 20 53 47 50 34 28 49 |  22 CAL|L SGP4(I|
|00003ea0| 46 4c 41 47 2c 54 53 49 | 4e 43 45 29 0a 20 20 20 |FLAG,TSI|NCE).   |
|00003eb0| 20 20 20 47 4f 20 54 4f | 20 36 30 0a 20 20 20 32 |   GO TO| 60.   2|
|00003ec0| 33 20 43 41 4c 4c 20 53 | 44 50 34 28 49 46 4c 41 |3 CALL S|DP4(IFLA|
|00003ed0| 47 2c 54 53 49 4e 43 45 | 29 0a 20 20 20 20 20 20 |G,TSINCE|).      |
|00003ee0| 47 4f 20 54 4f 20 36 30 | 0a 20 20 20 32 34 20 43 |GO TO 60|.   24 C|
|00003ef0| 41 4c 4c 20 53 47 50 38 | 28 49 46 4c 41 47 2c 54 |ALL SGP8|(IFLAG,T|
|00003f00| 53 49 4e 43 45 29 0a 20 | 20 20 20 20 20 47 4f 20 |SINCE). |     GO |
|00003f10| 54 4f 20 36 30 0a 20 20 | 20 32 35 20 43 41 4c 4c |TO 60.  | 25 CALL|
|00003f20| 20 53 44 50 38 28 49 46 | 4c 41 47 2c 54 53 49 4e | SDP8(IF|LAG,TSIN|
|00003f30| 43 45 29 0a 20 20 20 36 | 30 20 58 3d 58 2a 58 4b |CE).   6|0 X=X*XK|
|00003f40| 4d 50 45 52 2f 41 45 0a | 20 20 20 20 20 20 59 3d |MPER/AE.|      Y=|
|00003f50| 59 2a 58 4b 4d 50 45 52 | 2f 41 45 0a 20 20 20 20 |Y*XKMPER|/AE.    |
|00003f60| 20 20 5a 3d 5a 2a 58 4b | 4d 50 45 52 2f 41 45 0a |  Z=Z*XK|MPER/AE.|
|00003f70| 20 20 20 20 20 20 58 44 | 4f 54 3d 58 44 4f 54 2a |      XD|OT=XDOT*|
|00003f80| 58 4b 4d 50 45 52 2f 41 | 45 2a 58 4d 4e 50 44 41 |XKMPER/A|E*XMNPDA|
|00003f90| 2f 38 36 34 30 30 2e 0a | 20 20 20 20 20 20 59 44 |/86400..|      YD|
|00003fa0| 4f 54 3d 59 44 4f 54 2a | 58 4b 4d 50 45 52 2f 41 |OT=YDOT*|XKMPER/A|
|00003fb0| 45 2a 58 4d 4e 50 44 41 | 2f 38 36 34 30 30 2e 0a |E*XMNPDA|/86400..|
|00003fc0| 20 20 20 20 20 20 5a 44 | 4f 54 3d 5a 44 4f 54 2a |      ZD|OT=ZDOT*|
|00003fd0| 58 4b 4d 50 45 52 2f 41 | 45 2a 58 4d 4e 50 44 41 |XKMPER/A|E*XMNPDA|
|00003fe0| 2f 38 36 34 30 30 2e 0a | 20 20 20 20 20 20 57 52 |/86400..|      WR|
|00003ff0| 49 54 45 28 36 2c 37 30 | 35 29 20 54 53 49 4e 43 |ITE(6,70|5) TSINC|
|00004000| 45 2c 58 2c 59 2c 5a 2c | 58 44 4f 54 2c 59 44 4f |E,X,Y,Z,|XDOT,YDO|
|00004010| 54 2c 5a 44 4f 54 0a 20 | 20 20 20 20 20 54 53 49 |T,ZDOT. |     TSI|
|00004020| 4e 43 45 3d 54 53 49 4e | 43 45 2b 44 45 4c 54 0a |NCE=TSIN|CE+DELT.|
|00004030| 20 20 20 20 20 20 49 46 | 28 41 42 53 28 54 53 49 |      IF|(ABS(TSI|
|00004040| 4e 43 45 29 20 2e 47 54 | 2e 20 41 42 53 28 54 46 |NCE) .GT|. ABS(TF|
|00004050| 29 29 20 47 4f 20 54 4f | 20 32 0a 20 20 20 20 20 |)) GO TO| 2.     |
|00004060| 20 47 4f 20 54 4f 20 31 | 30 0a 20 20 37 30 30 20 | GO TO 1|0.  700 |
|00004070| 46 4f 52 4d 41 54 28 49 | 31 2c 33 46 31 30 2e 30 |FORMAT(I|1,3F10.0|
|00004080| 29 0a 20 20 37 30 31 20 | 46 4f 52 4d 41 54 28 32 |).  701 |FORMAT(2|
|00004090| 39 58 2c 44 31 34 2e 38 | 2c 31 58 2c 33 46 38 2e |9X,D14.8|,1X,3F8.|
|000040a0| 34 2c 2f 2c 36 58 2c 46 | 38 2e 37 2c 46 38 2e 34 |4,/,6X,F|8.7,F8.4|
|000040b0| 2c 31 58 2c 32 46 31 31 | 2e 39 2c 31 58 2c 46 36 |,1X,2F11|.9,1X,F6|
|000040c0| 2e 35 2c 49 32 2c 0a 20 | 20 20 20 20 31 20 20 20 |.5,I2,. |    1   |
|000040d0| 20 20 20 20 34 58 2c 46 | 38 2e 37 2c 49 32 29 0a |    4X,F|8.7,I2).|
|000040e0| 20 20 37 30 32 20 46 4f | 52 4d 41 54 28 31 38 58 |  702 FO|RMAT(18X|
|000040f0| 2c 44 31 34 2e 38 2c 31 | 58 2c 46 31 30 2e 38 2c |,D14.8,1|X,F10.8,|
|00004100| 32 28 31 58 2c 46 36 2e | 35 2c 49 32 29 2c 2f 2c |2(1X,F6.|5,I2),/,|
|00004110| 37 58 2c 32 28 31 58 2c | 46 38 2e 34 29 2c 31 58 |7X,2(1X,|F8.4),1X|
|00004120| 2c 0a 20 20 20 20 20 31 | 20 20 20 20 20 20 20 46 |,.     1|       F|
|00004130| 37 2e 37 2c 32 28 31 58 | 2c 46 38 2e 34 29 2c 31 |7.7,2(1X|,F8.4),1|
|00004140| 58 2c 46 31 31 2e 38 29 | 0a 20 20 37 30 33 20 46 |X,F11.8)|.  703 F|
|00004150| 4f 52 4d 41 54 28 37 39 | 58 2c 41 31 29 0a 20 20 |ORMAT(79|X,A1).  |
|00004160| 37 30 34 20 46 4f 52 4d | 41 54 28 31 48 31 2c 41 |704 FORM|AT(1H1,A|
|00004170| 38 30 2c 2f 2c 31 58 2c | 41 38 30 2c 2f 2f 2c 31 |80,/,1X,|A80,//,1|
|00004180| 58 2c 41 34 2c 37 48 20 | 54 53 49 4e 43 45 2c 0a |X,A4,7H |TSINCE,.|
|00004190| 20 20 20 20 20 31 20 31 | 34 58 2c 31 48 58 2c 31 |     1 1|4X,1HX,1|
|000041a0| 36 58 2c 31 48 59 2c 31 | 36 58 2c 31 48 5a 2c 31 |6X,1HY,1|6X,1HZ,1|
|000041b0| 34 58 2c 0a 20 20 20 20 | 20 31 20 34 48 58 44 4f |4X,.    | 1 4HXDO|
|000041c0| 54 2c 31 33 58 2c 34 48 | 59 44 4f 54 2c 31 33 58 |T,13X,4H|YDOT,13X|
|000041d0| 2c 34 48 5a 44 4f 54 2c | 2f 2f 29 0a 20 20 37 30 |,4HZDOT,|//).  70|
|000041e0| 35 20 46 4f 52 4d 41 54 | 28 37 46 31 37 2e 38 29 |5 FORMAT|(7F17.8)|
|000041f0| 0a 20 20 37 30 36 20 46 | 4f 52 4d 41 54 28 41 38 |.  706 F|ORMAT(A8|
|00004200| 30 29 0a 20 20 37 30 37 | 20 46 4f 52 4d 41 54 28 |0).  707| FORMAT(|
|00004210| 37 39 58 2c 41 31 29 0a | 20 20 39 33 30 20 46 4f |79X,A1).|  930 FO|
|00004220| 52 4d 41 54 28 22 53 48 | 4f 55 4c 44 20 55 53 45 |RMAT("SH|OULD USE|
|00004230| 20 44 45 45 50 20 53 50 | 41 43 45 20 45 50 48 45 | DEEP SP|ACE EPHE|
|00004240| 4d 45 52 49 53 22 29 0a | 20 20 39 34 30 20 46 4f |MERIS").|  940 FO|
|00004250| 52 4d 41 54 28 22 53 48 | 4f 55 4c 44 20 55 53 45 |RMAT("SH|OULD USE|
|00004260| 20 4e 45 41 52 20 45 41 | 52 54 48 20 45 50 48 45 | NEAR EA|RTH EPHE|
|00004270| 4d 45 52 49 53 22 29 0a | 20 20 39 35 30 20 46 4f |MERIS").|  950 FO|
|00004280| 52 4d 41 54 28 22 45 50 | 48 45 4d 45 52 49 53 20 |RMAT("EP|HEMERIS |
|00004290| 4e 55 4d 42 45 52 22 2c | 49 32 2c 22 20 4e 4f 54 |NUMBER",|I2," NOT|
|000042a0| 20 4c 45 47 41 4c 2c 20 | 57 49 4c 4c 20 53 4b 49 | LEGAL, |WILL SKI|
|000042b0| 50 20 54 48 49 53 20 43 | 41 53 45 22 29 0a 20 20 |P THIS C|ASE").  |
|000042c0| 38 30 30 20 57 52 49 54 | 45 28 36 2c 39 33 30 29 |800 WRIT|E(6,930)|
|000042d0| 0a 20 20 20 20 20 20 47 | 4f 20 54 4f 20 39 0a 20 |.      G|O TO 9. |
|000042e0| 20 38 35 30 20 57 52 49 | 54 45 28 36 2c 39 34 30 | 850 WRI|TE(6,940|
|000042f0| 29 0a 20 20 20 20 20 20 | 47 4f 20 54 4f 20 31 30 |).      |GO TO 10|
|00004300| 0a 20 20 39 30 30 20 57 | 52 49 54 45 28 36 2c 39 |.  900 W|RITE(6,9|
|00004310| 35 30 29 20 49 45 50 54 | 0a 20 20 20 20 20 20 47 |50) IEPT|.      G|
|00004320| 4f 20 54 4f 20 32 0a 20 | 20 20 20 20 20 45 4e 44 |O TO 2. |     END|
|00004330| 0a 5c 65 6e 64 7b 76 65 | 72 62 61 74 69 6d 7d 0a |.\end{ve|rbatim}.|
|00004340| 2d 2d 20 46 4d 4f 44 32 | 50 2e 46 4f 52 20 2d 2d |-- FMOD2|P.FOR --|
|00004350| 0a 5c 62 65 67 69 6e 7b | 76 65 72 62 61 74 69 6d |.\begin{|verbatim|
|00004360| 7d 0a 20 20 20 20 20 20 | 46 55 4e 43 54 49 4f 4e |}.      |FUNCTION|
|00004370| 20 46 4d 4f 44 32 50 28 | 58 29 0a 20 20 20 20 20 | FMOD2P(|X).     |
|00004380| 20 43 4f 4d 4d 4f 4e 2f | 43 32 2f 44 45 32 52 41 | COMMON/|C2/DE2RA|
|00004390| 2c 50 49 2c 50 49 4f 32 | 2c 54 57 4f 50 49 2c 58 |,PI,PIO2|,TWOPI,X|
|000043a0| 33 50 49 4f 32 0a 20 20 | 20 20 20 20 46 4d 4f 44 |3PIO2.  |    FMOD|
|000043b0| 32 50 3d 58 0a 20 20 20 | 20 20 20 49 3d 46 4d 4f |2P=X.   |   I=FMO|
|000043c0| 44 32 50 2f 54 57 4f 50 | 49 0a 20 20 20 20 20 20 |D2P/TWOP|I.      |
|000043d0| 20 46 4d 4f 44 32 50 3d | 46 4d 4f 44 32 50 2d 49 | FMOD2P=|FMOD2P-I|
|000043e0| 2a 54 57 4f 50 49 0a 20 | 20 20 20 20 20 49 46 28 |*TWOPI. |     IF(|
|000043f0| 46 4d 4f 44 32 50 2e 4c | 54 2e 30 29 20 46 4d 4f |FMOD2P.L|T.0) FMO|
|00004400| 44 32 50 3d 46 4d 4f 44 | 32 50 2b 54 57 4f 50 49 |D2P=FMOD|2P+TWOPI|
|00004410| 0a 20 20 20 20 20 20 52 | 45 54 55 52 4e 0a 20 20 |.      R|ETURN.  |
|00004420| 20 20 20 20 45 4e 44 0a | 5c 65 6e 64 7b 76 65 72 |    END.|\end{ver|
|00004430| 62 61 74 69 6d 7d 0a 2d | 2d 20 53 47 50 2e 46 4f |batim}.-|- SGP.FO|
|00004440| 52 20 2d 2d 0a 5c 62 65 | 67 69 6e 7b 76 65 72 62 |R --.\be|gin{verb|
|00004450| 61 74 69 6d 7d 0a 2a 20 | 20 20 20 20 20 53 47 50 |atim}.* |     SGP|
|00004460| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |        |        |
|00004470| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |        |        |
|00004480| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 33 31 20 |        |     31 |
|00004490| 4f 43 54 20 38 30 0a 20 | 20 20 20 20 20 53 55 42 |OCT 80. |     SUB|
|000044a0| 52 4f 55 54 49 4e 45 20 | 53 47 50 28 49 46 4c 41 |ROUTINE |SGP(IFLA|
|000044b0| 47 2c 54 53 49 4e 43 45 | 29 0a 20 20 20 20 20 20 |G,TSINCE|).      |
|000044c0| 43 4f 4d 4d 4f 4e 2f 45 | 31 2f 58 4d 4f 2c 58 4e |COMMON/E|1/XMO,XN|
|000044d0| 4f 44 45 4f 2c 4f 4d 45 | 47 41 4f 2c 45 4f 2c 58 |ODEO,OME|GAO,EO,X|
|000044e0| 49 4e 43 4c 2c 58 4e 4f | 2c 58 4e 44 54 32 4f 2c |INCL,XNO|,XNDT2O,|
|000044f0| 58 4e 44 44 36 4f 2c 42 | 53 54 41 52 2c 0a 20 20 |XNDD6O,B|STAR,.  |
|00004500| 20 20 20 31 20 20 20 20 | 20 20 20 20 20 58 2c 59 |   1    |     X,Y|
|00004510| 2c 5a 2c 58 44 4f 54 2c | 59 44 4f 54 2c 5a 44 4f |,Z,XDOT,|YDOT,ZDO|
|00004520| 54 2c 45 50 4f 43 48 2c | 44 53 35 30 0a 20 20 20 |T,EPOCH,|DS50.   |
|00004530| 20 20 20 43 4f 4d 4d 4f | 4e 2f 43 31 2f 43 4b 32 |   COMMO|N/C1/CK2|
|00004540| 2c 43 4b 34 2c 45 36 41 | 2c 51 4f 4d 53 32 54 2c |,CK4,E6A|,QOMS2T,|
|00004550| 53 2c 54 4f 54 48 52 44 | 2c 0a 20 20 20 20 20 31 |S,TOTHRD|,.     1|
|00004560| 20 20 20 20 20 20 20 20 | 58 4a 33 2c 58 4b 45 2c |        |XJ3,XKE,|
|00004570| 58 4b 4d 50 45 52 2c 58 | 4d 4e 50 44 41 2c 41 45 |XKMPER,X|MNPDA,AE|
|00004580| 0a 20 20 20 20 20 20 44 | 4f 55 42 4c 45 20 50 52 |.      D|OUBLE PR|
|00004590| 45 43 49 53 49 4f 4e 20 | 45 50 4f 43 48 2c 20 44 |ECISION |EPOCH, D|
|000045a0| 53 35 30 0a 0a 20 20 20 | 20 20 20 49 46 28 49 46 |S50..   |   IF(IF|
|000045b0| 4c 41 47 2e 45 51 2e 30 | 29 20 47 4f 20 54 4f 20 |LAG.EQ.0|) GO TO |
|000045c0| 31 39 0a 0a 2a 20 20 20 | 20 20 20 49 4e 49 54 49 |19..*   |   INITI|
|000045d0| 41 4c 49 5a 41 54 49 4f | 4e 0a 0a 20 20 20 20 20 |ALIZATIO|N..     |
|000045e0| 20 43 31 3d 20 43 4b 32 | 2a 31 2e 35 0a 20 20 20 | C1= CK2|*1.5.   |
|000045f0| 20 20 20 43 32 3d 20 43 | 4b 32 2f 34 2e 30 0a 20 |   C2= C|K2/4.0. |
|00004600| 20 20 20 20 20 43 33 3d | 20 43 4b 32 2f 32 2e 30 |     C3=| CK2/2.0|
|00004610| 0a 20 20 20 20 20 20 43 | 34 3d 20 58 4a 33 2a 41 |.      C|4= XJ3*A|
|00004620| 45 2a 2a 33 2f 28 34 2e | 30 2a 43 4b 32 29 0a 20 |E**3/(4.|0*CK2). |
|00004630| 20 20 20 20 20 43 4f 53 | 49 4f 3d 43 4f 53 28 58 |     COS|IO=COS(X|
|00004640| 49 4e 43 4c 29 0a 20 20 | 20 20 20 20 53 49 4e 49 |INCL).  |    SINI|
|00004650| 4f 3d 53 49 4e 28 58 49 | 4e 43 4c 29 0a 20 20 20 |O=SIN(XI|NCL).   |
|00004660| 20 20 20 41 31 3d 28 58 | 4b 45 2f 58 4e 4f 29 2a |   A1=(X|KE/XNO)*|
|00004670| 2a 54 4f 54 48 52 44 0a | 20 20 20 20 20 20 44 31 |*TOTHRD.|      D1|
|00004680| 3d 20 20 20 20 20 43 31 | 2f 41 31 2f 41 31 2a 28 |=     C1|/A1/A1*(|
|00004690| 33 2e 2a 43 4f 53 49 4f | 2a 43 4f 53 49 4f 2d 31 |3.*COSIO|*COSIO-1|
|000046a0| 2e 29 2f 28 31 2e 2d 45 | 4f 2a 45 4f 29 2a 2a 31 |.)/(1.-E|O*EO)**1|
|000046b0| 2e 35 0a 20 20 20 20 20 | 20 41 4f 3d 41 31 2a 28 |.5.     | AO=A1*(|
|000046c0| 31 2e 2d 31 2e 2f 33 2e | 2a 44 31 2d 44 31 2a 44 |1.-1./3.|*D1-D1*D|
|000046d0| 31 2d 31 33 34 2e 2f 38 | 31 2e 2a 44 31 2a 44 31 |1-134./8|1.*D1*D1|
|000046e0| 2a 44 31 29 0a 20 20 20 | 20 20 20 50 4f 3d 41 4f |*D1).   |   PO=AO|
|000046f0| 2a 28 31 2e 2d 45 4f 2a | 45 4f 29 0a 20 20 20 20 |*(1.-EO*|EO).    |
|00004700| 20 20 51 4f 3d 41 4f 2a | 28 31 2e 2d 45 4f 29 0a |  QO=AO*|(1.-EO).|
|00004710| 20 20 20 20 20 20 58 4c | 4f 3d 58 4d 4f 2b 4f 4d |      XL|O=XMO+OM|
|00004720| 45 47 41 4f 2b 58 4e 4f | 44 45 4f 0a 20 20 20 20 |EGAO+XNO|DEO.    |
|00004730| 20 20 44 31 4f 3d 20 43 | 33 20 2a 53 49 4e 49 4f |  D1O= C|3 *SINIO|
|00004740| 2a 53 49 4e 49 4f 0a 20 | 20 20 20 20 20 44 32 4f |*SINIO. |     D2O|
|00004750| 3d 20 43 32 20 2a 28 37 | 2e 2a 43 4f 53 49 4f 2a |= C2 *(7|.*COSIO*|
|00004760| 43 4f 53 49 4f 2d 31 2e | 29 0a 20 20 20 20 20 20 |COSIO-1.|).      |
|00004770| 44 33 4f 3d 43 31 2a 43 | 4f 53 49 4f 0a 20 20 20 |D3O=C1*C|OSIO.   |
|00004780| 20 20 20 44 34 4f 3d 44 | 33 4f 2a 53 49 4e 49 4f |   D4O=D|3O*SINIO|
|00004790| 0a 20 20 20 20 20 20 50 | 4f 32 4e 4f 3d 58 4e 4f |.      P|O2NO=XNO|
|000047a0| 2f 28 50 4f 2a 50 4f 29 | 0a 20 20 20 20 20 20 4f |/(PO*PO)|.      O|
|000047b0| 4d 47 44 54 3d 43 31 2a | 50 4f 32 4e 4f 2a 28 35 |MGDT=C1*|PO2NO*(5|
|000047c0| 2e 2a 43 4f 53 49 4f 2a | 43 4f 53 49 4f 2d 31 2e |.*COSIO*|COSIO-1.|
|000047d0| 29 0a 20 20 20 20 20 20 | 58 4e 4f 44 4f 54 3d 2d |).      |XNODOT=-|
|000047e0| 32 2e 2a 44 33 4f 2a 50 | 4f 32 4e 4f 0a 20 20 20 |2.*D3O*P|O2NO.   |
|000047f0| 20 20 20 43 35 3d 2e 35 | 2a 43 34 2a 53 49 4e 49 |   C5=.5|*C4*SINI|
|00004800| 4f 2a 28 33 2e 2b 35 2e | 2a 43 4f 53 49 4f 29 2f |O*(3.+5.|*COSIO)/|
|00004810| 28 31 2e 2b 43 4f 53 49 | 4f 29 0a 20 20 20 20 20 |(1.+COSI|O).     |
|00004820| 20 43 36 3d 43 34 2a 53 | 49 4e 49 4f 0a 20 20 20 | C6=C4*S|INIO.   |
|00004830| 20 20 20 49 46 4c 41 47 | 3d 30 0a 0a 2a 20 20 20 |   IFLAG|=0..*   |
|00004840| 20 20 20 55 50 44 41 54 | 45 20 46 4f 52 20 53 45 |   UPDAT|E FOR SE|
|00004850| 43 55 4c 41 52 20 47 52 | 41 56 49 54 59 20 41 4e |CULAR GR|AVITY AN|
|00004860| 44 20 41 54 4d 4f 53 50 | 48 45 52 49 43 20 44 52 |D ATMOSP|HERIC DR|
|00004870| 41 47 0a 0a 20 20 20 31 | 39 20 41 3d 58 4e 4f 2b |AG..   1|9 A=XNO+|
|00004880| 28 32 2e 2a 58 4e 44 54 | 32 4f 2b 33 2e 2a 58 4e |(2.*XNDT|2O+3.*XN|
|00004890| 44 44 36 4f 2a 54 53 49 | 4e 43 45 29 2a 54 53 49 |DD6O*TSI|NCE)*TSI|
|000048a0| 4e 43 45 0a 20 20 20 20 | 20 20 41 3d 41 4f 2a 28 |NCE.    |  A=AO*(|
|000048b0| 58 4e 4f 2f 41 29 2a 2a | 54 4f 54 48 52 44 0a 20 |XNO/A)**|TOTHRD. |
|000048c0| 20 20 20 20 20 45 3d 45 | 36 41 0a 20 20 20 20 20 |     E=E|6A.     |
|000048d0| 20 49 46 28 41 2e 47 54 | 2e 51 4f 29 20 45 3d 31 | IF(A.GT|.QO) E=1|
|000048e0| 2e 2d 51 4f 2f 41 0a 20 | 20 20 20 20 20 50 3d 41 |.-QO/A. |     P=A|
|000048f0| 2a 28 31 2e 2d 45 2a 45 | 29 0a 20 20 20 20 20 20 |*(1.-E*E|).      |
|00004900| 58 4e 4f 44 45 53 3d 20 | 58 4e 4f 44 45 4f 2b 58 |XNODES= |XNODEO+X|
|00004910| 4e 4f 44 4f 54 2a 54 53 | 49 4e 43 45 0a 20 20 20 |NODOT*TS|INCE.   |
|00004920| 20 20 20 4f 4d 47 41 53 | 3d 20 4f 4d 45 47 41 4f |   OMGAS|= OMEGAO|
|00004930| 2b 4f 4d 47 44 54 2a 54 | 53 49 4e 43 45 0a 20 20 |+OMGDT*T|SINCE.  |
|00004940| 20 20 20 20 58 4c 53 3d | 46 4d 4f 44 32 50 28 58 |    XLS=|FMOD2P(X|
|00004950| 4c 4f 2b 28 58 4e 4f 2b | 4f 4d 47 44 54 2b 58 4e |LO+(XNO+|OMGDT+XN|
|00004960| 4f 44 4f 54 2b 28 58 4e | 44 54 32 4f 2b 58 4e 44 |ODOT+(XN|DT2O+XND|
|00004970| 44 36 4f 2a 54 53 49 4e | 43 45 29 2a 0a 20 20 20 |D6O*TSIN|CE)*.   |
|00004980| 20 20 31 20 54 53 49 4e | 43 45 29 2a 54 53 49 4e |  1 TSIN|CE)*TSIN|
|00004990| 43 45 29 0a 0a 2a 20 20 | 20 20 20 20 4c 4f 4e 47 |CE)..*  |    LONG|
|000049a0| 20 50 45 52 49 4f 44 20 | 50 45 52 49 4f 44 49 43 | PERIOD |PERIODIC|
|000049b0| 53 0a 0a 20 20 20 20 20 | 20 41 58 4e 53 4c 3d 45 |S..     | AXNSL=E|
|000049c0| 2a 43 4f 53 28 4f 4d 47 | 41 53 29 0a 20 20 20 20 |*COS(OMG|AS).    |
|000049d0| 20 20 41 59 4e 53 4c 3d | 45 2a 53 49 4e 28 4f 4d |  AYNSL=|E*SIN(OM|
|000049e0| 47 41 53 29 2d 43 36 2f | 50 0a 20 20 20 20 20 20 |GAS)-C6/|P.      |
|000049f0| 58 4c 3d 46 4d 4f 44 32 | 50 28 58 4c 53 2d 43 35 |XL=FMOD2|P(XLS-C5|
|00004a00| 2f 50 2a 41 58 4e 53 4c | 29 0a 0a 2a 20 20 20 20 |/P*AXNSL|)..*    |
|00004a10| 20 20 53 4f 4c 56 45 20 | 4b 45 50 4c 45 52 53 20 |  SOLVE |KEPLERS |
|00004a20| 45 51 55 41 54 49 4f 4e | 0a 0a 20 20 20 20 20 20 |EQUATION|..      |
|00004a30| 55 3d 46 4d 4f 44 32 50 | 28 58 4c 2d 58 4e 4f 44 |U=FMOD2P|(XL-XNOD|
|00004a40| 45 53 29 0a 20 20 20 20 | 20 20 49 54 45 4d 33 3d |ES).    |  ITEM3=|
|00004a50| 30 0a 20 20 20 20 20 20 | 45 4f 31 3d 55 0a 20 20 |0.      |EO1=U.  |
|00004a60| 20 20 20 20 54 45 4d 35 | 3d 31 2e 0a 20 20 20 32 |    TEM5|=1..   2|
|00004a70| 30 20 53 49 4e 45 4f 31 | 3d 53 49 4e 28 45 4f 31 |0 SINEO1|=SIN(EO1|
|00004a80| 29 0a 20 20 20 20 20 20 | 43 4f 53 45 4f 31 3d 43 |).      |COSEO1=C|
|00004a90| 4f 53 28 45 4f 31 29 0a | 20 20 20 20 20 20 49 46 |OS(EO1).|      IF|
|00004aa0| 28 41 42 53 28 54 45 4d | 35 29 2e 4c 54 2e 45 36 |(ABS(TEM|5).LT.E6|
|00004ab0| 41 29 20 47 4f 20 54 4f | 20 33 30 0a 20 20 20 20 |A) GO TO| 30.    |
|00004ac0| 20 20 49 46 28 49 54 45 | 4d 33 2e 47 45 2e 31 30 |  IF(ITE|M3.GE.10|
|00004ad0| 29 20 47 4f 20 54 4f 20 | 33 30 0a 20 20 20 20 20 |) GO TO |30.     |
|00004ae0| 20 49 54 45 4d 33 3d 49 | 54 45 4d 33 2b 31 0a 20 | ITEM3=I|TEM3+1. |
|00004af0| 20 20 20 20 20 54 45 4d | 35 3d 31 2e 2d 43 4f 53 |     TEM|5=1.-COS|
|00004b00| 45 4f 31 2a 41 58 4e 53 | 4c 2d 53 49 4e 45 4f 31 |EO1*AXNS|L-SINEO1|
|00004b10| 2a 41 59 4e 53 4c 0a 20 | 20 20 20 20 20 54 45 4d |*AYNSL. |     TEM|
|00004b20| 35 3d 28 55 2d 41 59 4e | 53 4c 2a 43 4f 53 45 4f |5=(U-AYN|SL*COSEO|
|00004b30| 31 2b 41 58 4e 53 4c 2a | 53 49 4e 45 4f 31 2d 45 |1+AXNSL*|SINEO1-E|
|00004b40| 4f 31 29 2f 54 45 4d 35 | 0a 20 20 20 20 20 20 54 |O1)/TEM5|.      T|
|00004b50| 45 4d 32 3d 41 42 53 28 | 54 45 4d 35 29 0a 20 20 |EM2=ABS(|TEM5).  |
|00004b60| 20 20 20 20 49 46 28 54 | 45 4d 32 2e 47 54 2e 31 |    IF(T|EM2.GT.1|
|00004b70| 2e 29 20 54 45 4d 35 3d | 54 45 4d 32 2f 54 45 4d |.) TEM5=|TEM2/TEM|
|00004b80| 35 0a 20 20 20 20 20 20 | 45 4f 31 3d 45 4f 31 2b |5.      |EO1=EO1+|
|00004b90| 54 45 4d 35 0a 20 20 20 | 20 20 20 47 4f 20 54 4f |TEM5.   |   GO TO|
|00004ba0| 20 32 30 0a 0a 2a 20 20 | 20 20 20 20 53 48 4f 52 | 20..*  |    SHOR|
|00004bb0| 54 20 50 45 52 49 4f 44 | 20 50 52 45 4c 49 4d 49 |T PERIOD| PRELIMI|
|00004bc0| 4e 41 52 59 20 51 55 41 | 4e 54 49 54 49 45 53 0a |NARY QUA|NTITIES.|
|00004bd0| 0a 20 20 20 33 30 20 45 | 43 4f 53 45 3d 41 58 4e |.   30 E|COSE=AXN|
|00004be0| 53 4c 2a 43 4f 53 45 4f | 31 2b 41 59 4e 53 4c 2a |SL*COSEO|1+AYNSL*|
|00004bf0| 53 49 4e 45 4f 31 0a 20 | 20 20 20 20 20 45 53 49 |SINEO1. |     ESI|
|00004c00| 4e 45 3d 41 58 4e 53 4c | 2a 53 49 4e 45 4f 31 2d |NE=AXNSL|*SINEO1-|
|00004c10| 41 59 4e 53 4c 2a 43 4f | 53 45 4f 31 0a 20 20 20 |AYNSL*CO|SEO1.   |
|00004c20| 20 20 20 45 4c 32 3d 41 | 58 4e 53 4c 2a 41 58 4e |   EL2=A|XNSL*AXN|
|00004c30| 53 4c 2b 41 59 4e 53 4c | 2a 41 59 4e 53 4c 0a 20 |SL+AYNSL|*AYNSL. |
|00004c40| 20 20 20 20 20 50 4c 3d | 41 2a 28 31 2e 2d 45 4c |     PL=|A*(1.-EL|
|00004c50| 32 29 0a 20 20 20 20 20 | 20 50 4c 32 3d 50 4c 2a |2).     | PL2=PL*|
|00004c60| 50 4c 0a 20 20 20 20 20 | 20 52 3d 41 2a 28 31 2e |PL.     | R=A*(1.|
|00004c70| 2d 45 43 4f 53 45 29 0a | 20 20 20 20 20 20 52 44 |-ECOSE).|      RD|
|00004c80| 4f 54 3d 58 4b 45 2a 53 | 51 52 54 28 41 29 2f 52 |OT=XKE*S|QRT(A)/R|
|00004c90| 2a 45 53 49 4e 45 0a 20 | 20 20 20 20 20 52 56 44 |*ESINE. |     RVD|
|00004ca0| 4f 54 3d 58 4b 45 2a 53 | 51 52 54 28 50 4c 29 2f |OT=XKE*S|QRT(PL)/|
|00004cb0| 52 0a 20 20 20 20 20 20 | 54 45 4d 50 3d 45 53 49 |R.      |TEMP=ESI|
|00004cc0| 4e 45 2f 28 31 2e 2b 53 | 51 52 54 28 31 2e 2d 45 |NE/(1.+S|QRT(1.-E|
|00004cd0| 4c 32 29 29 0a 20 20 20 | 20 20 20 53 49 4e 55 3d |L2)).   |   SINU=|
|00004ce0| 41 2f 52 2a 28 53 49 4e | 45 4f 31 2d 41 59 4e 53 |A/R*(SIN|EO1-AYNS|
|00004cf0| 4c 2d 41 58 4e 53 4c 2a | 54 45 4d 50 29 0a 20 20 |L-AXNSL*|TEMP).  |
|00004d00| 20 20 20 20 43 4f 53 55 | 3d 41 2f 52 2a 28 43 4f |    COSU|=A/R*(CO|
|00004d10| 53 45 4f 31 2d 41 58 4e | 53 4c 2b 41 59 4e 53 4c |SEO1-AXN|SL+AYNSL|
|00004d20| 2a 54 45 4d 50 29 0a 20 | 20 20 20 20 20 53 55 3d |*TEMP). |     SU=|
|00004d30| 41 43 54 41 4e 28 53 49 | 4e 55 2c 43 4f 53 55 29 |ACTAN(SI|NU,COSU)|
|00004d40| 0a 0a 2a 20 20 20 20 20 | 20 55 50 44 41 54 45 20 |..*     | UPDATE |
|00004d50| 46 4f 52 20 53 48 4f 52 | 54 20 50 45 52 49 4f 44 |FOR SHOR|T PERIOD|
|00004d60| 49 43 53 0a 0a 20 20 20 | 20 20 20 53 49 4e 32 55 |ICS..   |   SIN2U|
|00004d70| 3d 28 43 4f 53 55 2b 43 | 4f 53 55 29 2a 53 49 4e |=(COSU+C|OSU)*SIN|
|00004d80| 55 0a 20 20 20 20 20 20 | 43 4f 53 32 55 3d 31 2e |U.      |COS2U=1.|
|00004d90| 2d 32 2e 2a 53 49 4e 55 | 2a 53 49 4e 55 0a 20 20 |-2.*SINU|*SINU.  |
|00004da0| 20 20 20 20 52 4b 3d 52 | 2b 44 31 4f 2f 50 4c 2a |    RK=R|+D1O/PL*|
|00004db0| 43 4f 53 32 55 0a 20 20 | 20 20 20 20 55 4b 3d 53 |COS2U.  |    UK=S|
|00004dc0| 55 2d 44 32 4f 2f 50 4c | 32 2a 53 49 4e 32 55 0a |U-D2O/PL|2*SIN2U.|
|00004dd0| 20 20 20 20 20 20 58 4e | 4f 44 45 4b 3d 58 4e 4f |      XN|ODEK=XNO|
|00004de0| 44 45 53 2b 44 33 4f 2a | 53 49 4e 32 55 2f 50 4c |DES+D3O*|SIN2U/PL|
|00004df0| 32 0a 20 20 20 20 20 20 | 58 49 4e 43 4b 20 3d 58 |2.      |XINCK =X|
|00004e00| 49 4e 43 4c 2b 44 34 4f | 2f 50 4c 32 2a 43 4f 53 |INCL+D4O|/PL2*COS|
|00004e10| 32 55 0a 0a 2a 20 20 20 | 20 20 20 4f 52 49 45 4e |2U..*   |   ORIEN|
|00004e20| 54 41 54 49 4f 4e 20 56 | 45 43 54 4f 52 53 0a 0a |TATION V|ECTORS..|
|00004e30| 20 20 20 20 20 20 53 49 | 4e 55 4b 3d 53 49 4e 28 |      SI|NUK=SIN(|
|00004e40| 55 4b 29 0a 20 20 20 20 | 20 20 43 4f 53 55 4b 3d |UK).    |  COSUK=|
|00004e50| 43 4f 53 28 55 4b 29 0a | 20 20 20 20 20 20 53 49 |COS(UK).|      SI|
|00004e60| 4e 4e 4f 4b 3d 53 49 4e | 28 58 4e 4f 44 45 4b 29 |NNOK=SIN|(XNODEK)|
|00004e70| 0a 20 20 20 20 20 20 43 | 4f 53 4e 4f 4b 3d 43 4f |.      C|OSNOK=CO|
|00004e80| 53 28 58 4e 4f 44 45 4b | 29 0a 20 20 20 20 20 20 |S(XNODEK|).      |
|00004e90| 53 49 4e 49 4b 3d 53 49 | 4e 28 58 49 4e 43 4b 29 |SINIK=SI|N(XINCK)|
|00004ea0| 0a 20 20 20 20 20 20 43 | 4f 53 49 4b 3d 43 4f 53 |.      C|OSIK=COS|
|00004eb0| 28 58 49 4e 43 4b 29 0a | 20 20 20 20 20 20 58 4d |(XINCK).|      XM|
|00004ec0| 58 3d 2d 53 49 4e 4e 4f | 4b 2a 43 4f 53 49 4b 0a |X=-SINNO|K*COSIK.|
|00004ed0| 20 20 20 20 20 20 58 4d | 59 3d 43 4f 53 4e 4f 4b |      XM|Y=COSNOK|
|00004ee0| 2a 43 4f 53 49 4b 0a 20 | 20 20 20 20 20 55 58 3d |*COSIK. |     UX=|
|00004ef0| 58 4d 58 2a 53 49 4e 55 | 4b 2b 43 4f 53 4e 4f 4b |XMX*SINU|K+COSNOK|
|00004f00| 2a 43 4f 53 55 4b 0a 20 | 20 20 20 20 20 55 59 3d |*COSUK. |     UY=|
|00004f10| 58 4d 59 2a 53 49 4e 55 | 4b 2b 53 49 4e 4e 4f 4b |XMY*SINU|K+SINNOK|
|00004f20| 2a 43 4f 53 55 4b 0a 20 | 20 20 20 20 20 55 5a 3d |*COSUK. |     UZ=|
|00004f30| 53 49 4e 49 4b 2a 53 49 | 4e 55 4b 0a 20 20 20 20 |SINIK*SI|NUK.    |
|00004f40| 20 20 56 58 3d 58 4d 58 | 2a 43 4f 53 55 4b 2d 43 |  VX=XMX|*COSUK-C|
|00004f50| 4f 53 4e 4f 4b 2a 53 49 | 4e 55 4b 0a 20 20 20 20 |OSNOK*SI|NUK.    |
|00004f60| 20 20 56 59 3d 58 4d 59 | 2a 43 4f 53 55 4b 2d 53 |  VY=XMY|*COSUK-S|
|00004f70| 49 4e 4e 4f 4b 2a 53 49 | 4e 55 4b 0a 20 20 20 20 |INNOK*SI|NUK.    |
|00004f80| 20 20 56 5a 3d 53 49 4e | 49 4b 2a 43 4f 53 55 4b |  VZ=SIN|IK*COSUK|
|00004f90| 0a 0a 2a 20 20 20 20 20 | 20 50 4f 53 49 54 49 4f |..*     | POSITIO|
|00004fa0| 4e 20 41 4e 44 20 56 45 | 4c 4f 43 49 54 59 0a 0a |N AND VE|LOCITY..|
|00004fb0| 20 20 20 20 20 20 58 3d | 52 4b 2a 55 58 0a 20 20 |      X=|RK*UX.  |
|00004fc0| 20 20 20 20 59 3d 52 4b | 2a 55 59 0a 20 20 20 20 |    Y=RK|*UY.    |
|00004fd0| 20 20 5a 3d 52 4b 2a 55 | 5a 0a 20 20 20 20 20 20 |  Z=RK*U|Z.      |
|00004fe0| 58 44 4f 54 3d 52 44 4f | 54 2a 55 58 0a 20 20 20 |XDOT=RDO|T*UX.   |
|00004ff0| 20 20 20 59 44 4f 54 3d | 52 44 4f 54 2a 55 59 0a |   YDOT=|RDOT*UY.|
|00005000| 20 20 20 20 20 20 5a 44 | 4f 54 3d 52 44 4f 54 2a |      ZD|OT=RDOT*|
|00005010| 55 5a 0a 20 20 20 20 20 | 20 58 44 4f 54 3d 52 56 |UZ.     | XDOT=RV|
|00005020| 44 4f 54 2a 56 58 2b 58 | 44 4f 54 0a 20 20 20 20 |DOT*VX+X|DOT.    |
|00005030| 20 20 59 44 4f 54 3d 52 | 56 44 4f 54 2a 56 59 2b |  YDOT=R|VDOT*VY+|
|00005040| 59 44 4f 54 0a 20 20 20 | 20 20 20 5a 44 4f 54 3d |YDOT.   |   ZDOT=|
|00005050| 52 56 44 4f 54 2a 56 5a | 2b 5a 44 4f 54 0a 0a 20 |RVDOT*VZ|+ZDOT.. |
|00005060| 20 20 20 20 20 52 45 54 | 55 52 4e 0a 20 20 20 20 |     RET|URN.    |
|00005070| 20 20 45 4e 44 0a 5c 65 | 6e 64 7b 76 65 72 62 61 |  END.\e|nd{verba|
|00005080| 74 69 6d 7d 0a 2d 2d 20 | 53 47 50 34 2e 46 4f 52 |tim}.-- |SGP4.FOR|
|00005090| 20 2d 2d 0a 5c 62 65 67 | 69 6e 7b 76 65 72 62 61 | --.\beg|in{verba|
|000050a0| 74 69 6d 7d 0a 2a 20 20 | 20 20 20 20 53 47 50 34 |tim}.*  |    SGP4|
|000050b0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |        |        |
|000050c0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |        |        |
|000050d0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 33 |        |       3|
|000050e0| 20 4e 4f 56 20 38 30 0a | 20 20 20 20 20 20 53 55 | NOV 80.|      SU|
|000050f0| 42 52 4f 55 54 49 4e 45 | 20 53 47 50 34 28 49 46 |BROUTINE| SGP4(IF|
|00005100| 4c 41 47 2c 54 53 49 4e | 43 45 29 0a 20 20 20 20 |LAG,TSIN|CE).    |
|00005110| 20 20 43 4f 4d 4d 4f 4e | 2f 45 31 2f 58 4d 4f 2c |  COMMON|/E1/XMO,|
|00005120| 58 4e 4f 44 45 4f 2c 4f | 4d 45 47 41 4f 2c 45 4f |XNODEO,O|MEGAO,EO|
|00005130| 2c 58 49 4e 43 4c 2c 58 | 4e 4f 2c 58 4e 44 54 32 |,XINCL,X|NO,XNDT2|
|00005140| 4f 2c 0a 20 20 20 20 20 | 31 20 20 20 20 20 20 20 |O,.     |1       |
|00005150| 20 20 20 20 58 4e 44 44 | 36 4f 2c 42 53 54 41 52 |    XNDD|6O,BSTAR|
|00005160| 2c 58 2c 59 2c 5a 2c 58 | 44 4f 54 2c 59 44 4f 54 |,X,Y,Z,X|DOT,YDOT|
|00005170| 2c 5a 44 4f 54 2c 45 50 | 4f 43 48 2c 44 53 35 30 |,ZDOT,EP|OCH,DS50|
|00005180| 0a 20 20 20 20 20 20 43 | 4f 4d 4d 4f 4e 2f 43 31 |.      C|OMMON/C1|
|00005190| 2f 43 4b 32 2c 43 4b 34 | 2c 45 36 41 2c 51 4f 4d |/CK2,CK4|,E6A,QOM|
|000051a0| 53 32 54 2c 53 2c 54 4f | 54 48 52 44 2c 0a 20 20 |S2T,S,TO|THRD,.  |
|000051b0| 20 20 20 31 20 20 20 20 | 20 20 20 20 20 20 20 58 |   1    |       X|
|000051c0| 4a 33 2c 58 4b 45 2c 58 | 4b 4d 50 45 52 2c 58 4d |J3,XKE,X|KMPER,XM|
|000051d0| 4e 50 44 41 2c 41 45 0a | 20 20 20 20 20 20 44 4f |NPDA,AE.|      DO|
|000051e0| 55 42 4c 45 20 50 52 45 | 43 49 53 49 4f 4e 20 45 |UBLE PRE|CISION E|
|000051f0| 50 4f 43 48 2c 20 44 53 | 35 30 0a 0a 20 20 20 20 |POCH, DS|50..    |
|00005200| 20 20 49 46 20 28 49 46 | 4c 41 47 20 2e 45 51 2e |  IF (IF|LAG .EQ.|
|00005210| 20 30 29 20 47 4f 20 54 | 4f 20 31 30 30 0a 0a 2a | 0) GO T|O 100..*|
|00005220| 20 20 20 20 20 20 52 45 | 43 4f 56 45 52 20 4f 52 |      RE|COVER OR|
|00005230| 49 47 49 4e 41 4c 20 4d | 45 41 4e 20 4d 4f 54 49 |IGINAL M|EAN MOTI|
|00005240| 4f 4e 20 28 58 4e 4f 44 | 50 29 20 41 4e 44 20 53 |ON (XNOD|P) AND S|
|00005250| 45 4d 49 4d 41 4a 4f 52 | 20 41 58 49 53 20 28 41 |EMIMAJOR| AXIS (A|
|00005260| 4f 44 50 29 0a 2a 20 20 | 20 20 20 20 46 52 4f 4d |ODP).*  |    FROM|
|00005270| 20 49 4e 50 55 54 20 45 | 4c 45 4d 45 4e 54 53 0a | INPUT E|LEMENTS.|
|00005280| 0a 20 20 20 20 20 20 41 | 31 3d 28 58 4b 45 2f 58 |.      A|1=(XKE/X|
|00005290| 4e 4f 29 2a 2a 54 4f 54 | 48 52 44 0a 20 20 20 20 |NO)**TOT|HRD.    |
|000052a0| 20 20 43 4f 53 49 4f 3d | 43 4f 53 28 58 49 4e 43 |  COSIO=|COS(XINC|
|000052b0| 4c 29 0a 20 20 20 20 20 | 20 54 48 45 54 41 32 3d |L).     | THETA2=|
|000052c0| 43 4f 53 49 4f 2a 43 4f | 53 49 4f 0a 20 20 20 20 |COSIO*CO|SIO.    |
|000052d0| 20 20 58 33 54 48 4d 31 | 3d 33 2e 2a 54 48 45 54 |  X3THM1|=3.*THET|
|000052e0| 41 32 2d 31 2e 0a 20 20 | 20 20 20 20 45 4f 53 51 |A2-1..  |    EOSQ|
|000052f0| 3d 45 4f 2a 45 4f 0a 20 | 20 20 20 20 20 42 45 54 |=EO*EO. |     BET|
|00005300| 41 4f 32 3d 31 2e 2d 45 | 4f 53 51 0a 20 20 20 20 |AO2=1.-E|OSQ.    |
|00005310| 20 20 42 45 54 41 4f 3d | 53 51 52 54 28 42 45 54 |  BETAO=|SQRT(BET|
|00005320| 41 4f 32 29 0a 20 20 20 | 20 20 20 44 45 4c 31 3d |AO2).   |   DEL1=|
|00005330| 31 2e 35 2a 43 4b 32 2a | 58 33 54 48 4d 31 2f 28 |1.5*CK2*|X3THM1/(|
|00005340| 41 31 2a 41 31 2a 42 45 | 54 41 4f 2a 42 45 54 41 |A1*A1*BE|TAO*BETA|
|00005350| 4f 32 29 0a 20 20 20 20 | 20 20 41 4f 3d 41 31 2a |O2).    |  AO=A1*|
|00005360| 28 31 2e 2d 44 45 4c 31 | 2a 28 2e 35 2a 54 4f 54 |(1.-DEL1|*(.5*TOT|
|00005370| 48 52 44 2b 44 45 4c 31 | 2a 28 31 2e 2b 31 33 34 |HRD+DEL1|*(1.+134|
|00005380| 2e 2f 38 31 2e 2a 44 45 | 4c 31 29 29 29 0a 20 20 |./81.*DE|L1))).  |
|00005390| 20 20 20 20 44 45 4c 4f | 3d 31 2e 35 2a 43 4b 32 |    DELO|=1.5*CK2|
|000053a0| 2a 58 33 54 48 4d 31 2f | 28 41 4f 2a 41 4f 2a 42 |*X3THM1/|(AO*AO*B|
|000053b0| 45 54 41 4f 2a 42 45 54 | 41 4f 32 29 0a 20 20 20 |ETAO*BET|AO2).   |
|000053c0| 20 20 20 58 4e 4f 44 50 | 3d 58 4e 4f 2f 28 31 2e |   XNODP|=XNO/(1.|
|000053d0| 2b 44 45 4c 4f 29 0a 20 | 20 20 20 20 20 41 4f 44 |+DELO). |     AOD|
|000053e0| 50 3d 41 4f 2f 28 31 2e | 2d 44 45 4c 4f 29 0a 0a |P=AO/(1.|-DELO)..|
|000053f0| 2a 20 20 20 20 20 20 49 | 4e 49 54 49 41 4c 49 5a |*      I|NITIALIZ|
|00005400| 41 54 49 4f 4e 0a 0a 2a | 20 20 20 20 20 20 46 4f |ATION..*|      FO|
|00005410| 52 20 50 45 52 49 47 45 | 45 20 4c 45 53 53 20 54 |R PERIGE|E LESS T|
|00005420| 48 41 4e 20 32 32 30 20 | 4b 49 4c 4f 4d 45 54 45 |HAN 220 |KILOMETE|
|00005430| 52 53 2c 20 54 48 45 20 | 49 53 49 4d 50 20 46 4c |RS, THE |ISIMP FL|
|00005440| 41 47 20 49 53 20 53 45 | 54 20 41 4e 44 0a 2a 20 |AG IS SE|T AND.* |
|00005450| 20 20 20 20 20 54 48 45 | 20 45 51 55 41 54 49 4f |     THE| EQUATIO|
|00005460| 4e 53 20 41 52 45 20 54 | 52 55 4e 43 41 54 45 44 |NS ARE T|RUNCATED|
|00005470| 20 54 4f 20 4c 49 4e 45 | 41 52 20 56 41 52 49 41 | TO LINE|AR VARIA|
|00005480| 54 49 4f 4e 20 49 4e 20 | 53 51 52 54 20 41 20 41 |TION IN |SQRT A A|
|00005490| 4e 44 0a 2a 20 20 20 20 | 20 20 51 55 41 44 52 41 |ND.*    |  QUADRA|
|000054a0| 54 49 43 20 56 41 52 49 | 41 54 49 4f 4e 20 49 4e |TIC VARI|ATION IN|
|000054b0| 20 4d 45 41 4e 20 41 4e | 4f 4d 41 4c 59 2e 20 20 | MEAN AN|OMALY.  |
|000054c0| 41 4c 53 4f 2c 20 54 48 | 45 20 43 33 20 54 45 52 |ALSO, TH|E C3 TER|
|000054d0| 4d 2c 20 54 48 45 0a 2a | 20 20 20 20 20 20 44 45 |M, THE.*|      DE|
|000054e0| 4c 54 41 20 4f 4d 45 47 | 41 20 54 45 52 4d 2c 20 |LTA OMEG|A TERM, |
|000054f0| 41 4e 44 20 54 48 45 20 | 44 45 4c 54 41 20 4d 20 |AND THE |DELTA M |
|00005500| 54 45 52 4d 20 41 52 45 | 20 44 52 4f 50 50 45 44 |TERM ARE| DROPPED|
|00005510| 2e 0a 0a 20 20 20 20 20 | 20 49 53 49 4d 50 3d 30 |...     | ISIMP=0|
|00005520| 0a 20 20 20 20 20 20 49 | 46 28 28 41 4f 44 50 2a |.      I|F((AODP*|
|00005530| 28 31 2e 2d 45 4f 29 2f | 41 45 29 20 2e 4c 54 2e |(1.-EO)/|AE) .LT.|
|00005540| 20 28 32 32 30 2e 2f 58 | 4b 4d 50 45 52 2b 41 45 | (220./X|KMPER+AE|
|00005550| 29 29 20 49 53 49 4d 50 | 3d 31 0a 0a 2a 20 20 20 |)) ISIMP|=1..*   |
|00005560| 20 20 20 46 4f 52 20 50 | 45 52 49 47 45 45 20 42 |   FOR P|ERIGEE B|
|00005570| 45 4c 4f 57 20 31 35 36 | 20 4b 4d 2c 20 54 48 45 |ELOW 156| KM, THE|
|00005580| 20 56 41 4c 55 45 53 20 | 4f 46 0a 2a 20 20 20 20 | VALUES |OF.*    |
|00005590| 20 20 53 20 41 4e 44 20 | 51 4f 4d 53 32 54 20 41 |  S AND |QOMS2T A|
|000055a0| 52 45 20 41 4c 54 45 52 | 45 44 0a 0a 20 20 20 20 |RE ALTER|ED..    |
|000055b0| 20 20 53 34 3d 53 0a 20 | 20 20 20 20 20 51 4f 4d |  S4=S. |     QOM|
|000055c0| 53 32 34 3d 51 4f 4d 53 | 32 54 0a 20 20 20 20 20 |S24=QOMS|2T.     |
|000055d0| 20 50 45 52 49 47 45 3d | 28 41 4f 44 50 2a 28 31 | PERIGE=|(AODP*(1|
|000055e0| 2e 2d 45 4f 29 2d 41 45 | 29 2a 58 4b 4d 50 45 52 |.-EO)-AE|)*XKMPER|
|000055f0| 0a 20 20 20 20 20 20 49 | 46 28 50 45 52 49 47 45 |.      I|F(PERIGE|
|00005600| 20 2e 47 45 2e 20 31 35 | 36 2e 29 20 47 4f 20 54 | .GE. 15|6.) GO T|
|00005610| 4f 20 31 30 0a 20 20 20 | 20 20 20 53 34 3d 50 45 |O 10.   |   S4=PE|
|00005620| 52 49 47 45 2d 37 38 2e | 0a 20 20 20 20 20 20 49 |RIGE-78.|.      I|
|00005630| 46 28 50 45 52 49 47 45 | 20 2e 47 54 2e 20 39 38 |F(PERIGE| .GT. 98|
|00005640| 2e 29 20 47 4f 20 54 4f | 20 39 0a 20 20 20 20 20 |.) GO TO| 9.     |
|00005650| 20 53 34 3d 32 30 2e 0a | 20 20 20 20 39 20 51 4f | S4=20..|    9 QO|
|00005660| 4d 53 32 34 3d 28 28 31 | 32 30 2e 2d 53 34 29 2a |MS24=((1|20.-S4)*|
|00005670| 41 45 2f 58 4b 4d 50 45 | 52 29 2a 2a 34 0a 20 20 |AE/XKMPE|R)**4.  |
|00005680| 20 20 20 20 53 34 3d 53 | 34 2f 58 4b 4d 50 45 52 |    S4=S|4/XKMPER|
|00005690| 2b 41 45 0a 20 20 20 31 | 30 20 50 49 4e 56 53 51 |+AE.   1|0 PINVSQ|
|000056a0| 3d 31 2e 2f 28 41 4f 44 | 50 2a 41 4f 44 50 2a 42 |=1./(AOD|P*AODP*B|
|000056b0| 45 54 41 4f 32 2a 42 45 | 54 41 4f 32 29 0a 20 20 |ETAO2*BE|TAO2).  |
|000056c0| 20 20 20 20 54 53 49 3d | 31 2e 2f 28 41 4f 44 50 |    TSI=|1./(AODP|
|000056d0| 2d 53 34 29 0a 20 20 20 | 20 20 20 45 54 41 3d 41 |-S4).   |   ETA=A|
|000056e0| 4f 44 50 2a 45 4f 2a 54 | 53 49 0a 20 20 20 20 20 |ODP*EO*T|SI.     |
|000056f0| 20 45 54 41 53 51 3d 45 | 54 41 2a 45 54 41 0a 20 | ETASQ=E|TA*ETA. |
|00005700| 20 20 20 20 20 45 45 54 | 41 3d 45 4f 2a 45 54 41 |     EET|A=EO*ETA|
|00005710| 0a 20 20 20 20 20 20 50 | 53 49 53 51 3d 41 42 53 |.      P|SISQ=ABS|
|00005720| 28 31 2e 2d 45 54 41 53 | 51 29 0a 20 20 20 20 20 |(1.-ETAS|Q).     |
|00005730| 20 43 4f 45 46 3d 51 4f | 4d 53 32 34 2a 54 53 49 | COEF=QO|MS24*TSI|
|00005740| 2a 2a 34 0a 20 20 20 20 | 20 20 43 4f 45 46 31 3d |**4.    |  COEF1=|
|00005750| 43 4f 45 46 2f 50 53 49 | 53 51 2a 2a 33 2e 35 0a |COEF/PSI|SQ**3.5.|
|00005760| 20 20 20 20 20 20 43 32 | 3d 43 4f 45 46 31 2a 58 |      C2|=COEF1*X|
|00005770| 4e 4f 44 50 2a 28 41 4f | 44 50 2a 28 31 2e 2b 31 |NODP*(AO|DP*(1.+1|
|00005780| 2e 35 2a 45 54 41 53 51 | 2b 45 45 54 41 2a 28 34 |.5*ETASQ|+EETA*(4|
|00005790| 2e 2b 45 54 41 53 51 29 | 29 2b 2e 37 35 2a 0a 20 |.+ETASQ)|)+.75*. |
|000057a0| 20 20 20 20 31 20 20 20 | 20 20 20 20 20 20 43 4b |    1   |      CK|
|000057b0| 32 2a 54 53 49 2f 50 53 | 49 53 51 2a 58 33 54 48 |2*TSI/PS|ISQ*X3TH|
|000057c0| 4d 31 2a 28 38 2e 2b 33 | 2e 2a 45 54 41 53 51 2a |M1*(8.+3|.*ETASQ*|
|000057d0| 28 38 2e 2b 45 54 41 53 | 51 29 29 29 0a 20 20 20 |(8.+ETAS|Q))).   |
|000057e0| 20 20 20 43 31 3d 42 53 | 54 41 52 2a 43 32 0a 20 |   C1=BS|TAR*C2. |
|000057f0| 20 20 20 20 20 53 49 4e | 49 4f 3d 53 49 4e 28 58 |     SIN|IO=SIN(X|
|00005800| 49 4e 43 4c 29 0a 20 20 | 20 20 20 20 41 33 4f 56 |INCL).  |    A3OV|
|00005810| 4b 32 3d 2d 58 4a 33 2f | 43 4b 32 2a 41 45 2a 2a |K2=-XJ3/|CK2*AE**|
|00005820| 33 0a 20 20 20 20 20 20 | 43 33 3d 43 4f 45 46 2a |3.      |C3=COEF*|
|00005830| 54 53 49 2a 41 33 4f 56 | 4b 32 2a 58 4e 4f 44 50 |TSI*A3OV|K2*XNODP|
|00005840| 2a 41 45 2a 53 49 4e 49 | 4f 2f 45 4f 0a 20 20 20 |*AE*SINI|O/EO.   |
|00005850| 20 20 20 58 31 4d 54 48 | 32 3d 31 2e 2d 54 48 45 |   X1MTH|2=1.-THE|
|00005860| 54 41 32 0a 20 20 20 20 | 20 20 43 34 3d 32 2e 2a |TA2.    |  C4=2.*|
|00005870| 58 4e 4f 44 50 2a 43 4f | 45 46 31 2a 41 4f 44 50 |XNODP*CO|EF1*AODP|
|00005880| 2a 42 45 54 41 4f 32 2a | 28 45 54 41 2a 0a 20 20 |*BETAO2*|(ETA*.  |
|00005890| 20 20 20 31 20 20 20 20 | 20 20 20 20 20 28 32 2e |   1    |     (2.|
|000058a0| 2b 2e 35 2a 45 54 41 53 | 51 29 2b 45 4f 2a 28 2e |+.5*ETAS|Q)+EO*(.|
|000058b0| 35 2b 32 2e 2a 45 54 41 | 53 51 29 2d 32 2e 2a 43 |5+2.*ETA|SQ)-2.*C|
|000058c0| 4b 32 2a 54 53 49 2f 0a | 20 20 20 20 20 32 20 20 |K2*TSI/.|     2  |
|000058d0| 20 20 20 20 20 20 20 28 | 41 4f 44 50 2a 50 53 49 |       (|AODP*PSI|
|000058e0| 53 51 29 2a 28 2d 33 2e | 2a 58 33 54 48 4d 31 2a |SQ)*(-3.|*X3THM1*|
|000058f0| 28 31 2e 2d 32 2e 2a 45 | 45 54 41 2b 45 54 41 53 |(1.-2.*E|ETA+ETAS|
|00005900| 51 2a 0a 20 20 20 20 20 | 33 20 20 20 20 20 20 20 |Q*.     |3       |
|00005910| 20 20 28 31 2e 35 2d 2e | 35 2a 45 45 54 41 29 29 |  (1.5-.|5*EETA))|
|00005920| 2b 2e 37 35 2a 58 31 4d | 54 48 32 2a 28 32 2e 2a |+.75*X1M|TH2*(2.*|
|00005930| 45 54 41 53 51 2d 45 45 | 54 41 2a 0a 20 20 20 20 |ETASQ-EE|TA*.    |
|00005940| 20 34 20 20 20 20 20 20 | 20 20 20 28 31 2e 2b 45 | 4      |   (1.+E|
|00005950| 54 41 53 51 29 29 2a 43 | 4f 53 28 32 2e 2a 4f 4d |TASQ))*C|OS(2.*OM|
|00005960| 45 47 41 4f 29 29 29 0a | 20 20 20 20 20 20 43 35 |EGAO))).|      C5|
|00005970| 3d 32 2e 2a 43 4f 45 46 | 31 2a 41 4f 44 50 2a 42 |=2.*COEF|1*AODP*B|
|00005980| 45 54 41 4f 32 2a 28 31 | 2e 2b 32 2e 37 35 2a 28 |ETAO2*(1|.+2.75*(|
|00005990| 45 54 41 53 51 2b 45 45 | 54 41 29 2b 45 45 54 41 |ETASQ+EE|TA)+EETA|
|000059a0| 2a 45 54 41 53 51 29 0a | 20 20 20 20 20 20 54 48 |*ETASQ).|      TH|
|000059b0| 45 54 41 34 3d 54 48 45 | 54 41 32 2a 54 48 45 54 |ETA4=THE|TA2*THET|
|000059c0| 41 32 0a 20 20 20 20 20 | 20 54 45 4d 50 31 3d 33 |A2.     | TEMP1=3|
|000059d0| 2e 2a 43 4b 32 2a 50 49 | 4e 56 53 51 2a 58 4e 4f |.*CK2*PI|NVSQ*XNO|
|000059e0| 44 50 0a 20 20 20 20 20 | 20 54 45 4d 50 32 3d 54 |DP.     | TEMP2=T|
|000059f0| 45 4d 50 31 2a 43 4b 32 | 2a 50 49 4e 56 53 51 0a |EMP1*CK2|*PINVSQ.|
|00005a00| 20 20 20 20 20 20 54 45 | 4d 50 33 3d 31 2e 32 35 |      TE|MP3=1.25|
|00005a10| 2a 43 4b 34 2a 50 49 4e | 56 53 51 2a 50 49 4e 56 |*CK4*PIN|VSQ*PINV|
|00005a20| 53 51 2a 58 4e 4f 44 50 | 0a 20 20 20 20 20 20 58 |SQ*XNODP|.      X|
|00005a30| 4d 44 4f 54 3d 58 4e 4f | 44 50 2b 2e 35 2a 54 45 |MDOT=XNO|DP+.5*TE|
|00005a40| 4d 50 31 2a 42 45 54 41 | 4f 2a 58 33 54 48 4d 31 |MP1*BETA|O*X3THM1|
|00005a50| 2b 2e 30 36 32 35 2a 54 | 45 4d 50 32 2a 42 45 54 |+.0625*T|EMP2*BET|
|00005a60| 41 4f 2a 0a 20 20 20 20 | 20 31 20 20 20 20 20 20 |AO*.    | 1      |
|00005a70| 20 20 20 28 31 33 2e 2d | 37 38 2e 2a 54 48 45 54 |   (13.-|78.*THET|
|00005a80| 41 32 2b 31 33 37 2e 2a | 54 48 45 54 41 34 29 0a |A2+137.*|THETA4).|
|00005a90| 20 20 20 20 20 20 58 31 | 4d 35 54 48 3d 31 2e 2d |      X1|M5TH=1.-|
|00005aa0| 35 2e 2a 54 48 45 54 41 | 32 0a 20 20 20 20 20 20 |5.*THETA|2.      |
|00005ab0| 4f 4d 47 44 4f 54 3d 2d | 2e 35 2a 54 45 4d 50 31 |OMGDOT=-|.5*TEMP1|
|00005ac0| 2a 58 31 4d 35 54 48 2b | 2e 30 36 32 35 2a 54 45 |*X1M5TH+|.0625*TE|
|00005ad0| 4d 50 32 2a 28 37 2e 2d | 31 31 34 2e 2a 54 48 45 |MP2*(7.-|114.*THE|
|00005ae0| 54 41 32 2b 0a 20 20 20 | 20 20 31 20 20 20 20 20 |TA2+.   |  1     |
|00005af0| 20 20 20 20 33 39 35 2e | 2a 54 48 45 54 41 34 29 |    395.|*THETA4)|
|00005b00| 2b 54 45 4d 50 33 2a 28 | 33 2e 2d 33 36 2e 2a 54 |+TEMP3*(|3.-36.*T|
|00005b10| 48 45 54 41 32 2b 34 39 | 2e 2a 54 48 45 54 41 34 |HETA2+49|.*THETA4|
|00005b20| 29 0a 20 20 20 20 20 20 | 58 48 44 4f 54 31 3d 2d |).      |XHDOT1=-|
|00005b30| 54 45 4d 50 31 2a 43 4f | 53 49 4f 0a 20 20 20 20 |TEMP1*CO|SIO.    |
|00005b40| 20 20 58 4e 4f 44 4f 54 | 3d 58 48 44 4f 54 31 2b |  XNODOT|=XHDOT1+|
|00005b50| 28 2e 35 2a 54 45 4d 50 | 32 2a 28 34 2e 2d 31 39 |(.5*TEMP|2*(4.-19|
|00005b60| 2e 2a 54 48 45 54 41 32 | 29 2b 32 2e 2a 54 45 4d |.*THETA2|)+2.*TEM|
|00005b70| 50 33 2a 28 33 2e 2d 0a | 20 20 20 20 20 31 20 20 |P3*(3.-.|     1  |
|00005b80| 20 20 20 20 20 20 20 37 | 2e 2a 54 48 45 54 41 32 |       7|.*THETA2|
|00005b90| 29 29 2a 43 4f 53 49 4f | 0a 20 20 20 20 20 20 4f |))*COSIO|.      O|
|00005ba0| 4d 47 43 4f 46 3d 42 53 | 54 41 52 2a 43 33 2a 43 |MGCOF=BS|TAR*C3*C|
|00005bb0| 4f 53 28 4f 4d 45 47 41 | 4f 29 0a 20 20 20 20 20 |OS(OMEGA|O).     |
|00005bc0| 20 58 4d 43 4f 46 3d 2d | 54 4f 54 48 52 44 2a 43 | XMCOF=-|TOTHRD*C|
|00005bd0| 4f 45 46 2a 42 53 54 41 | 52 2a 41 45 2f 45 45 54 |OEF*BSTA|R*AE/EET|
|00005be0| 41 0a 20 20 20 20 20 20 | 58 4e 4f 44 43 46 3d 33 |A.      |XNODCF=3|
|00005bf0| 2e 35 2a 42 45 54 41 4f | 32 2a 58 48 44 4f 54 31 |.5*BETAO|2*XHDOT1|
|00005c00| 2a 43 31 0a 20 20 20 20 | 20 20 54 32 43 4f 46 3d |*C1.    |  T2COF=|
|00005c10| 31 2e 35 2a 43 31 0a 20 | 20 20 20 20 20 58 4c 43 |1.5*C1. |     XLC|
|00005c20| 4f 46 3d 2e 31 32 35 2a | 41 33 4f 56 4b 32 2a 53 |OF=.125*|A3OVK2*S|
|00005c30| 49 4e 49 4f 2a 28 33 2e | 2b 35 2e 2a 43 4f 53 49 |INIO*(3.|+5.*COSI|
|00005c40| 4f 29 2f 28 31 2e 2b 43 | 4f 53 49 4f 29 0a 20 20 |O)/(1.+C|OSIO).  |
|00005c50| 20 20 20 20 41 59 43 4f | 46 3d 2e 32 35 2a 41 33 |    AYCO|F=.25*A3|
|00005c60| 4f 56 4b 32 2a 53 49 4e | 49 4f 0a 20 20 20 20 20 |OVK2*SIN|IO.     |
|00005c70| 20 44 45 4c 4d 4f 3d 28 | 31 2e 2b 45 54 41 2a 43 | DELMO=(|1.+ETA*C|
|00005c80| 4f 53 28 58 4d 4f 29 29 | 2a 2a 33 0a 20 20 20 20 |OS(XMO))|**3.    |
|00005c90| 20 20 53 49 4e 4d 4f 3d | 53 49 4e 28 58 4d 4f 29 |  SINMO=|SIN(XMO)|
|00005ca0| 0a 20 20 20 20 20 20 58 | 37 54 48 4d 31 3d 37 2e |.      X|7THM1=7.|
|00005cb0| 2a 54 48 45 54 41 32 2d | 31 2e 0a 20 20 20 20 20 |*THETA2-|1..     |
|00005cc0| 20 49 46 28 49 53 49 4d | 50 20 2e 45 51 2e 20 31 | IF(ISIM|P .EQ. 1|
|00005cd0| 29 20 47 4f 20 54 4f 20 | 39 30 0a 20 20 20 20 20 |) GO TO |90.     |
|00005ce0| 20 43 31 53 51 3d 43 31 | 2a 43 31 0a 20 20 20 20 | C1SQ=C1|*C1.    |
|00005cf0| 20 20 44 32 3d 34 2e 2a | 41 4f 44 50 2a 54 53 49 |  D2=4.*|AODP*TSI|
|00005d00| 2a 43 31 53 51 0a 20 20 | 20 20 20 20 54 45 4d 50 |*C1SQ.  |    TEMP|
|00005d10| 3d 44 32 2a 54 53 49 2a | 43 31 2f 33 2e 0a 20 20 |=D2*TSI*|C1/3..  |
|00005d20| 20 20 20 20 44 33 3d 28 | 31 37 2e 2a 41 4f 44 50 |    D3=(|17.*AODP|
|00005d30| 2b 53 34 29 2a 54 45 4d | 50 0a 20 20 20 20 20 20 |+S4)*TEM|P.      |
|00005d40| 44 34 3d 2e 35 2a 54 45 | 4d 50 2a 41 4f 44 50 2a |D4=.5*TE|MP*AODP*|
|00005d50| 54 53 49 2a 28 32 32 31 | 2e 2a 41 4f 44 50 2b 33 |TSI*(221|.*AODP+3|
|00005d60| 31 2e 2a 53 34 29 2a 43 | 31 0a 20 20 20 20 20 20 |1.*S4)*C|1.      |
|00005d70| 54 33 43 4f 46 3d 44 32 | 2b 32 2e 2a 43 31 53 51 |T3COF=D2|+2.*C1SQ|
|00005d80| 0a 20 20 20 20 20 20 54 | 34 43 4f 46 3d 2e 32 35 |.      T|4COF=.25|
|00005d90| 2a 28 33 2e 2a 44 33 2b | 43 31 2a 28 31 32 2e 2a |*(3.*D3+|C1*(12.*|
|00005da0| 44 32 2b 31 30 2e 2a 43 | 31 53 51 29 29 0a 20 20 |D2+10.*C|1SQ)).  |
|00005db0| 20 20 20 20 54 35 43 4f | 46 3d 2e 32 2a 28 33 2e |    T5CO|F=.2*(3.|
|00005dc0| 2a 44 34 2b 31 32 2e 2a | 43 31 2a 44 33 2b 36 2e |*D4+12.*|C1*D3+6.|
|00005dd0| 2a 44 32 2a 44 32 2b 31 | 35 2e 2a 43 31 53 51 2a |*D2*D2+1|5.*C1SQ*|
|00005de0| 28 0a 20 20 20 20 20 31 | 20 20 20 20 20 20 20 20 |(.     1|        |
|00005df0| 20 32 2e 2a 44 32 2b 43 | 31 53 51 29 29 0a 20 20 | 2.*D2+C|1SQ)).  |
|00005e00| 20 39 30 20 49 46 4c 41 | 47 3d 30 0a 0a 2a 20 20 | 90 IFLA|G=0..*  |
|00005e10| 20 20 20 20 55 50 44 41 | 54 45 20 46 4f 52 20 53 |    UPDA|TE FOR S|
|00005e20| 45 43 55 4c 41 52 20 47 | 52 41 56 49 54 59 20 41 |ECULAR G|RAVITY A|
|00005e30| 4e 44 20 41 54 4d 4f 53 | 50 48 45 52 49 43 20 44 |ND ATMOS|PHERIC D|
|00005e40| 52 41 47 0a 0a 20 20 31 | 30 30 20 58 4d 44 46 3d |RAG..  1|00 XMDF=|
|00005e50| 58 4d 4f 2b 58 4d 44 4f | 54 2a 54 53 49 4e 43 45 |XMO+XMDO|T*TSINCE|
|00005e60| 0a 20 20 20 20 20 20 4f | 4d 47 41 44 46 3d 4f 4d |.      O|MGADF=OM|
|00005e70| 45 47 41 4f 2b 4f 4d 47 | 44 4f 54 2a 54 53 49 4e |EGAO+OMG|DOT*TSIN|
|00005e80| 43 45 0a 20 20 20 20 20 | 20 58 4e 4f 44 44 46 3d |CE.     | XNODDF=|
|00005e90| 58 4e 4f 44 45 4f 2b 58 | 4e 4f 44 4f 54 2a 54 53 |XNODEO+X|NODOT*TS|
|00005ea0| 49 4e 43 45 0a 20 20 20 | 20 20 20 4f 4d 45 47 41 |INCE.   |   OMEGA|
|00005eb0| 3d 4f 4d 47 41 44 46 0a | 20 20 20 20 20 20 58 4d |=OMGADF.|      XM|
|00005ec0| 50 3d 58 4d 44 46 0a 20 | 20 20 20 20 20 54 53 51 |P=XMDF. |     TSQ|
|00005ed0| 3d 54 53 49 4e 43 45 2a | 54 53 49 4e 43 45 0a 20 |=TSINCE*|TSINCE. |
|00005ee0| 20 20 20 20 20 58 4e 4f | 44 45 3d 58 4e 4f 44 44 |     XNO|DE=XNODD|
|00005ef0| 46 2b 58 4e 4f 44 43 46 | 2a 54 53 51 0a 20 20 20 |F+XNODCF|*TSQ.   |
|00005f00| 20 20 20 54 45 4d 50 41 | 3d 31 2e 2d 43 31 2a 54 |   TEMPA|=1.-C1*T|
|00005f10| 53 49 4e 43 45 0a 20 20 | 20 20 20 20 54 45 4d 50 |SINCE.  |    TEMP|
|00005f20| 45 3d 42 53 54 41 52 2a | 43 34 2a 54 53 49 4e 43 |E=BSTAR*|C4*TSINC|
|00005f30| 45 0a 20 20 20 20 20 20 | 54 45 4d 50 4c 3d 54 32 |E.      |TEMPL=T2|
|00005f40| 43 4f 46 2a 54 53 51 0a | 20 20 20 20 20 20 49 46 |COF*TSQ.|      IF|
|00005f50| 28 49 53 49 4d 50 20 2e | 45 51 2e 20 31 29 20 47 |(ISIMP .|EQ. 1) G|
|00005f60| 4f 20 54 4f 20 31 31 30 | 0a 20 20 20 20 20 20 44 |O TO 110|.      D|
|00005f70| 45 4c 4f 4d 47 3d 4f 4d | 47 43 4f 46 2a 54 53 49 |ELOMG=OM|GCOF*TSI|
|00005f80| 4e 43 45 0a 20 20 20 20 | 20 20 44 45 4c 4d 3d 58 |NCE.    |  DELM=X|
|00005f90| 4d 43 4f 46 2a 28 28 31 | 2e 2b 45 54 41 2a 43 4f |MCOF*((1|.+ETA*CO|
|00005fa0| 53 28 58 4d 44 46 29 29 | 2a 2a 33 2d 44 45 4c 4d |S(XMDF))|**3-DELM|
|00005fb0| 4f 29 0a 20 20 20 20 20 | 20 54 45 4d 50 3d 44 45 |O).     | TEMP=DE|
|00005fc0| 4c 4f 4d 47 2b 44 45 4c | 4d 0a 20 20 20 20 20 20 |LOMG+DEL|M.      |
|00005fd0| 58 4d 50 3d 58 4d 44 46 | 2b 54 45 4d 50 0a 20 20 |XMP=XMDF|+TEMP.  |
|00005fe0| 20 20 20 20 4f 4d 45 47 | 41 3d 4f 4d 47 41 44 46 |    OMEG|A=OMGADF|
|00005ff0| 2d 54 45 4d 50 0a 20 20 | 20 20 20 20 54 43 55 42 |-TEMP.  |    TCUB|
|00006000| 45 3d 54 53 51 2a 54 53 | 49 4e 43 45 0a 20 20 20 |E=TSQ*TS|INCE.   |
|00006010| 20 20 20 54 46 4f 55 52 | 3d 54 53 49 4e 43 45 2a |   TFOUR|=TSINCE*|
|00006020| 54 43 55 42 45 0a 20 20 | 20 20 20 20 54 45 4d 50 |TCUBE.  |    TEMP|
|00006030| 41 3d 54 45 4d 50 41 2d | 44 32 2a 54 53 51 2d 44 |A=TEMPA-|D2*TSQ-D|
|00006040| 33 2a 54 43 55 42 45 2d | 44 34 2a 54 46 4f 55 52 |3*TCUBE-|D4*TFOUR|
|00006050| 0a 20 20 20 20 20 20 54 | 45 4d 50 45 3d 54 45 4d |.      T|EMPE=TEM|
|00006060| 50 45 2b 42 53 54 41 52 | 2a 43 35 2a 28 53 49 4e |PE+BSTAR|*C5*(SIN|
|00006070| 28 58 4d 50 29 2d 53 49 | 4e 4d 4f 29 0a 20 20 20 |(XMP)-SI|NMO).   |
|00006080| 20 20 20 54 45 4d 50 4c | 3d 54 45 4d 50 4c 2b 54 |   TEMPL|=TEMPL+T|
|00006090| 33 43 4f 46 2a 54 43 55 | 42 45 2b 0a 20 20 20 20 |3COF*TCU|BE+.    |
|000060a0| 20 31 20 20 20 20 20 20 | 20 20 20 54 46 4f 55 52 | 1      |   TFOUR|
|000060b0| 2a 28 54 34 43 4f 46 2b | 54 53 49 4e 43 45 2a 54 |*(T4COF+|TSINCE*T|
|000060c0| 35 43 4f 46 29 0a 20 20 | 31 31 30 20 41 3d 41 4f |5COF).  |110 A=AO|
|000060d0| 44 50 2a 54 45 4d 50 41 | 2a 2a 32 0a 20 20 20 20 |DP*TEMPA|**2.    |
|000060e0| 20 20 45 3d 45 4f 2d 54 | 45 4d 50 45 0a 20 20 20 |  E=EO-T|EMPE.   |
|000060f0| 20 20 20 58 4c 3d 58 4d | 50 2b 4f 4d 45 47 41 2b |   XL=XM|P+OMEGA+|
|00006100| 58 4e 4f 44 45 2b 58 4e | 4f 44 50 2a 54 45 4d 50 |XNODE+XN|ODP*TEMP|
|00006110| 4c 0a 20 20 20 20 20 20 | 42 45 54 41 3d 53 51 52 |L.      |BETA=SQR|
|00006120| 54 28 31 2e 2d 45 2a 45 | 29 0a 20 20 20 20 20 20 |T(1.-E*E|).      |
|00006130| 58 4e 3d 58 4b 45 2f 41 | 2a 2a 31 2e 35 0a 0a 2a |XN=XKE/A|**1.5..*|
|00006140| 20 20 20 20 20 20 4c 4f | 4e 47 20 50 45 52 49 4f |      LO|NG PERIO|
|00006150| 44 20 50 45 52 49 4f 44 | 49 43 53 0a 0a 20 20 20 |D PERIOD|ICS..   |
|00006160| 20 20 20 41 58 4e 3d 45 | 2a 43 4f 53 28 4f 4d 45 |   AXN=E|*COS(OME|
|00006170| 47 41 29 0a 20 20 20 20 | 20 20 54 45 4d 50 3d 31 |GA).    |  TEMP=1|
|00006180| 2e 2f 28 41 2a 42 45 54 | 41 2a 42 45 54 41 29 0a |./(A*BET|A*BETA).|
|00006190| 20 20 20 20 20 20 58 4c | 4c 3d 54 45 4d 50 2a 58 |      XL|L=TEMP*X|
|000061a0| 4c 43 4f 46 2a 41 58 4e | 0a 20 20 20 20 20 20 41 |LCOF*AXN|.      A|
|000061b0| 59 4e 4c 3d 54 45 4d 50 | 2a 41 59 43 4f 46 0a 20 |YNL=TEMP|*AYCOF. |
|000061c0| 20 20 20 20 20 58 4c 54 | 3d 58 4c 2b 58 4c 4c 0a |     XLT|=XL+XLL.|
|000061d0| 20 20 20 20 20 20 41 59 | 4e 3d 45 2a 53 49 4e 28 |      AY|N=E*SIN(|
|000061e0| 4f 4d 45 47 41 29 2b 41 | 59 4e 4c 0a 0a 2a 20 20 |OMEGA)+A|YNL..*  |
|000061f0| 20 20 20 20 53 4f 4c 56 | 45 20 4b 45 50 4c 45 52 |    SOLV|E KEPLER|
|00006200| 53 20 45 51 55 41 54 49 | 4f 4e 0a 0a 20 20 20 20 |S EQUATI|ON..    |
|00006210| 20 20 43 41 50 55 3d 46 | 4d 4f 44 32 50 28 58 4c |  CAPU=F|MOD2P(XL|
|00006220| 54 2d 58 4e 4f 44 45 29 | 0a 20 20 20 20 20 20 54 |T-XNODE)|.      T|
|00006230| 45 4d 50 32 3d 43 41 50 | 55 0a 20 20 20 20 20 20 |EMP2=CAP|U.      |
|00006240| 44 4f 20 31 33 30 20 49 | 3d 31 2c 31 30 0a 20 20 |DO 130 I|=1,10.  |
|00006250| 20 20 20 20 53 49 4e 45 | 50 57 3d 53 49 4e 28 54 |    SINE|PW=SIN(T|
|00006260| 45 4d 50 32 29 0a 20 20 | 20 20 20 20 43 4f 53 45 |EMP2).  |    COSE|
|00006270| 50 57 3d 43 4f 53 28 54 | 45 4d 50 32 29 0a 20 20 |PW=COS(T|EMP2).  |
|00006280| 20 20 20 20 54 45 4d 50 | 33 3d 41 58 4e 2a 53 49 |    TEMP|3=AXN*SI|
|00006290| 4e 45 50 57 0a 20 20 20 | 20 20 20 54 45 4d 50 34 |NEPW.   |   TEMP4|
|000062a0| 3d 41 59 4e 2a 43 4f 53 | 45 50 57 0a 20 20 20 20 |=AYN*COS|EPW.    |
|000062b0| 20 20 54 45 4d 50 35 3d | 41 58 4e 2a 43 4f 53 45 |  TEMP5=|AXN*COSE|
|000062c0| 50 57 0a 20 20 20 20 20 | 20 54 45 4d 50 36 3d 41 |PW.     | TEMP6=A|
|000062d0| 59 4e 2a 53 49 4e 45 50 | 57 0a 20 20 20 20 20 20 |YN*SINEP|W.      |
|000062e0| 45 50 57 3d 28 43 41 50 | 55 2d 54 45 4d 50 34 2b |EPW=(CAP|U-TEMP4+|
|000062f0| 54 45 4d 50 33 2d 54 45 | 4d 50 32 29 2f 28 31 2e |TEMP3-TE|MP2)/(1.|
|00006300| 2d 54 45 4d 50 35 2d 54 | 45 4d 50 36 29 2b 54 45 |-TEMP5-T|EMP6)+TE|
|00006310| 4d 50 32 0a 20 20 20 20 | 20 20 49 46 28 41 42 53 |MP2.    |  IF(ABS|
|00006320| 28 45 50 57 2d 54 45 4d | 50 32 29 20 2e 4c 45 2e |(EPW-TEM|P2) .LE.|
|00006330| 20 45 36 41 29 20 47 4f | 20 54 4f 20 31 34 30 0a | E6A) GO| TO 140.|
|00006340| 20 20 31 33 30 20 54 45 | 4d 50 32 3d 45 50 57 0a |  130 TE|MP2=EPW.|
|00006350| 0a 2a 20 20 20 20 20 20 | 53 48 4f 52 54 20 50 45 |.*      |SHORT PE|
|00006360| 52 49 4f 44 20 50 52 45 | 4c 49 4d 49 4e 41 52 59 |RIOD PRE|LIMINARY|
|00006370| 20 51 55 41 4e 54 49 54 | 49 45 53 0a 0a 20 20 31 | QUANTIT|IES..  1|
|00006380| 34 30 20 45 43 4f 53 45 | 3d 54 45 4d 50 35 2b 54 |40 ECOSE|=TEMP5+T|
|00006390| 45 4d 50 36 0a 20 20 20 | 20 20 20 45 53 49 4e 45 |EMP6.   |   ESINE|
|000063a0| 3d 54 45 4d 50 33 2d 54 | 45 4d 50 34 0a 20 20 20 |=TEMP3-T|EMP4.   |
|000063b0| 20 20 20 45 4c 53 51 3d | 41 58 4e 2a 41 58 4e 2b |   ELSQ=|AXN*AXN+|
|000063c0| 41 59 4e 2a 41 59 4e 0a | 20 20 20 20 20 20 54 45 |AYN*AYN.|      TE|
|000063d0| 4d 50 3d 31 2e 2d 45 4c | 53 51 0a 20 20 20 20 20 |MP=1.-EL|SQ.     |
|000063e0| 20 50 4c 3d 41 2a 54 45 | 4d 50 0a 20 20 20 20 20 | PL=A*TE|MP.     |
|000063f0| 20 52 3d 41 2a 28 31 2e | 2d 45 43 4f 53 45 29 0a | R=A*(1.|-ECOSE).|
+--------+-------------------------+-------------------------+--------+--------+